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SUMMARY 


The  purpose  of  the  Coreless  Linear  Conduction  Pump  (CLCP)  was  to 
evaluate  the  feasibility  of  the  CLCP  as  a  means  of  building 
Thermoelectromagnetic  Pumps  (TEMPs)  that  are  not,  in  themselves, 
temperatue  limited  by  the  Curie  point  of  ferromagnetic  cores  or  do  not 
limit  the  temperature  capability  of  the  Nuclear  Space  Power  Systems 
(NSPSs)  in  which  they  are  used.  A  TEMP  is  a  single  component  in  which  a 
thermoelectric  generator  (TEG)  and  an  electromagnetic  pump  (EMP)  are 
highly  integrated . 

The  CLCP  concept  employs  coreless  linear  magnets  rather  than 
circular  magnets  with  or  without  ferromagnetic  cores.  Also  in  the 
concept,  the  counter-current,  parallel  busbars  which  lie  on  either  side 
of  the  pump  duct  and  that  constitute  a  linear  magnet  carries  the 
integrated  field  and  armature  currents  to  provide  the  magnetic  field  at 
the  entrances  to  the  magnet  and  the  pumping  section.  As  the  current 
proceeds  down  one  busbar  a  portion  of  it  becomes  armature  current  by 
crossing  the  pump  duct  to  the  other  busbar  while  the  remaining  field 
current  continues  down  the  busbar.  The  partitioning  of  the  currents 
happens  all  along  the  busbars  and  duct  until  it  is  all  consumed. 

The  optimum  partitioning  factor  of  the  total  current  to  provide 
armature  currents  was  found  to  be  ~l/2  at  each  point  along  the  pump  duct 
and  Linear  Magnet  (LM)  length.  This  leads  to  geometric  degradation  of 
pumping  capability  from  entrcince  to  exit  so  that  the  upstream  end  of  the 
pump  produces  nearly  all  of  the  pumping  power.  Truncation  of  the 
unproductive  end  leads  to  higher  mass  efficiency,  but  leaves  the 
situation  nearly  the  same  as  occurs  without  the  integrated  current 
concept  and  without  the  complications  and  difficulties  of  detailed 
control  of  the  currents.  The  integrated  current  concept  was  abandoned 
in  favor  of  separate  field  and  armature  circuits. 

However,  coreless  magnets,  including  the  Helmholtz  and  linear 
magnets,  require  substantially  greater  electric  power  for  their 
operation  which  dominates  the  pump's  power  requirement.  This,  in  turn, 
requires  a  TEG  with  a  specific  mass  four  to  five  times  greater  than  that 
of  the  coreless  pump.  In  fact,  there  is  no  way  that  such  TEMPs  can 
provide  enough  heat  exchanger  surface  to  provide  the  TEG  with  the  heat 
rate  needed  for  the  TEG  to  meet  the  power  demand.  Thus,  coreless  TEMP 
were  not  found  to  be  feasible.  A  pumping  subsystem,  that  is  not  a  TEMP, 
in  which  separate  components  consisting  of  a  coreless  EMP,  heat 
exchanger,  and  TEG  may  provide  the  functions  of  a  TEMP  and  may  be 
feasible  but  too  massive  for  NSPSs. 

The  substitution  of  ZrBz  for  Mo  as  the  electrical  conductors  in  the 
TEMPs  at  1200K  clearly  indicated  the  superiority  of  ZrBz  in  this  appli¬ 
cation.  The  HM-TEMP,  using  ZrBz  as  a  conductor,  was  23%  lower  in  mass 
at  the  same  pumping  capacity  than  the  same  type  TEMP  optimized  for  the 
use  of  Mo  conductors.  However,  the  improvement  was  not  sufficient  to 
make  coreless  TEMPs  feasible.  The  CLCP  evaluated  here  is  associated 
with  a  64  kWe  NSPS.  For  larger  NSPS  of  interest  to  SDI,  coreless  TEMPs 
move  still  further  away  from  feasibility. 
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INTRODUCTION 


The  purpose  of  this  program  is  to  evaluate  the  Coreless  Linear 
Conduction  Pump  (CLCP)  concept  as  a  potential  component  of  a  nuclear 
space  power  system  that  could  result  in  significant  improvement  in 
survivability  and  mass  reduction  of  such  systems.  The  CLCP  is  an 
electromagnetic  pump  which,  when  powered  by  a  thermoeletric  power 
supply,  can  circulate  the  liquid  metal  coolant  during  normal  operation 
and/or  after  reactor  shutdown  for  as  long  as  there  is  decay  heat  in  the 
nuclear  reactor  to  be  removed.  The  significant  factor  in  realization  of 
the  potential  advantages  is  that  such  a  pump  would  eliminate  the  current 
temperature  limits  of  the  Curie  point  of  the  ferromagnetic  cores  of  the 
elec tromagnetic  pumps  previously  used,  or  being  developed  for  use,  in 
the  next  or  future  generations  of  nuclear  space  power  systems.  In  this 
project,  no  detailed  consideration  is  given  to  the  thermoelectric  power 
source  that  would  be  ultimately  integrated  with  the  pump  to  form  a  TEMP, 
if  the  CLCP  is  a  viable  concept. 

BACKGROUND 

In  1967  the  USA  launched  into  orbit  its  first  and,  so  far,  only 
space  power  system,  SNAPSHOT,  employing  a  nuclear  reactor  as  the  primary 
energy  source.  The  experimental  system  was  designed  to  produce  500 
W(e).  It  operated  for  only  42  days  before  it  was  inadvertently  shut 
down.  A  thermoelectric  electromagnetic  pump  (or  a  thermoelectromagnetic 
pump  (TEMP))  was  employed  as  the  main  pump  to  circulate  the  coolant  both 
during  normal  operation  and  after  reactor  shutdown.  The  magnetic  field 
in  the  TEMP  was  produced  by  a  permanent  magnet.  It  is  not  known  how 
long  the  TEMP  continued  to  circulate  the  coolant  after  shutdown. 

The  USA,  as  part  of  its  Strategic  Defense  Initiative,  is  currently 
developing  a  new  and  larger  (10  —  1000  kW(e))  Nulear  Space  Power  System 
(NSPS)  in  the  SP— 100  Program.  Its  main  and  after — shutdown  coolant 
circulation  pumps  will  be  TEMPs  using  electromagnets  to  produce  the 
magnetic  field. 

In  a  previous  study<^^,  TEMPs  using  permanent  magnets,  electro^ 
magnets,  and  Helmholtz  coils  (coreless)  to  produce  the  magnetic  fields 
were  evaluated  for  SP-lOO-type  NSPSs  over  a  range  of  powers  from  100 
kW(e)  to  10  MW(b).  The  coreless  Helmholtz  coils  did  not  perform  very 
well  in  producing  the  magnetic  field.  They  were  excessively  large  and 
massive.  Their  circular  nature  did  not  fit  well  to  the  linear  pumping 
section  where  the  magnetic  flux  was  needed.  From  this  observation  arose 
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the  concept  of  the  CLCP  ae  described  below. 

THE  CORELEBS  LINEAR  PUMP  CONCEPT 

The  basic  concept  of  the  coreless  pump  is  shown  in  Figure  1. 
geometry  of  the  coreless  pump  has  some  similarities,  at  least  at 
first  glance,  with  that  of  rail  guns.  The  two  major  current 

carrying  busbars  are  similar  to  the  rails  of  the  guns.  Current 
crosses  from  one  busbar  or  rail  to  the  other  as  it  moves  away  from 
the  power  source.  However,  in  the  rail  gun,  the  cross  conductor 
(plasma)  moves  down  the  rail  with  time  in  contrast  to  the  coreless 
pump  where  the  pump  cell  remains  in  place,  and  liquid  thrust  out  of 
the  cell  is  replaced  by  liquid  flowing  in.  Current  in  the  linear 
coreless  conduction  pump  crosses  all  along  the  busbars  from  one  to 
the  other. 

The  coreless  pump  concept  can  be  described  as  incorporating  a 
rectangular  pump  duct  of  height,  H,  width,  W,  and  length,  LT 
through  which  a  conductive  liquid  metal  (such  as  the  alkali  metals) 
flows  and  in  which  the  pumping  action  is  applied.  Busbars,  on 
opposite  sides  of  the  duct  width,  run  along  the  length  of  the  duct. 
The  busbars  are  attached  along  the  entire  length  of  the  duct  so  as 
to  supply  electric  current  for  passage  across  the  width  of  the  duct 
and  normal  to  the  flow  and  pumping  direction.  The  busbars  are 
connected  from  one  end  of  the  duct  to  the  power  supply  terminals, 
thereby  requiring  the  current  to  flow  in  opposite  directions  in 
each  busbar.  With  this  orientation,  the  magnetic  fluxes  induced  by 
the  currents  in  the  busbars  are  additive  between  the  busbars,  i.e., 
in  the  duct,  and  are  normal  to  both  the  current  and  the  flow. 
Thus,  the  conditions  required  for  pumping  are  created.  Outside  the 
busbars,  the  fluxes  induced  by  the  currents  in  the  busbars  are 
subtractive  and  diminish  rapidly  away  from  the  busbars. 

As  the  current  enters  the  positive  busbar,  it  flows  across  the 
duct  to  the  negative  busbar.  Assuming  constant  duct  and  busbar 
dimensions,  this  diminishes  the  density  of  the  current  continuing 
along  the  busbar  before  crossing  the  duct  further  along.  Thus, 
both  the  flux  and  the  pumping  current  fall  off  along  the  duct 
moving  away  from  the  power  source.  By  tailoring  the  thicknesses 
and  dimensions  from  end  to  end  of  the  duct,  it  is  possible  to 
control  the  flux  and  current  density  crossing  between  busbars  over 
the  length  of  the  pump.  The  width  and  height  of  the  duct  can  be 
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varied  along  its  length.  Also  the  thicknesses,  i.e.,  the  width,  of 
the  duct  walls  and  busbars  can  be  varied  over  the  length  ot  the 
pump.  Other  adjustments  are  possible  by  selection  of  the  materials 
of  the  wall  and  the  busbars. 


Figure  1.  Coreless  Linear  Conduction  Pump  Concept 
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The  second  major  objective  was  to  design,  build,  and  test  a  CLCP  to 
verify  the  concept. 

There  were  also  a  number  of  more  detailed  objectives  such  ass 

1)  Identify  critical  issues  related  to  TEMPs,  incorporating  CLCPs 

2)  To  compare  TEMPs  incorporating  CLCP  with  those  incorporating 
EMPs,  using  permanent  magnets,  electromagnets,  or  Helmholtz 
magnets 

3)  Identify  the  temperature  limiting  items  of  CLCP— TEMPs 

4)  Determine  the  comparable  survivability  of  CLCP-TEMPs 

5)  Identify  the  busbar  and  duct  materials  to  use  in  CLCP-TEMPs 
above  1300K 

+6)  Determine  if  CLCP  masses  are  more  than  those  of  cored  EMPs  of 
the  same  size 

7)  Determine  if  ZrBs  as  a  conductor  material  is  a  potential 
alternative  to  Mo 

TASKS 

The  original  statement  of  work  contained  five  tasks.  Tasks  3  and  4 
were  deleted  and  Task  6  was  added  by  amendment  of  the  contract  when  the 
results  of  the  design  and  analyses  indicated  that  testing  was  not 
warranted.  Nevertheless,  some  work  was  done  on  the  deleted  tasks  prior 
to  their  deletion,  and  it  is  reported  below. 

The  statement  of  work  for  each  task  is  presented  below  immediately 
followed  by  a  description  of  the  work  done,  a  statement  of  results,  a 
discussion  of  its  significance,  and  any  decisions  made  in  consequence. 

An  M2  thermoelec tromagnetic  pump  (TEMP)  specification  was  initially 
selected  as  a  pump  that  would  be  used  on  the  smallest  nuclear  space 
power  system  (NSPS)  of  interest  to  SDI .  This  specification  is  one  of  a 
family  of  TEMP  specifications  developed  in  a  previous  TEMP  study* and 
defines  one  of  three  4  psi  —  100  gpm  TEMPs  that  would  be  needed  in  a  64 
kWe  NSPS  using  Li  as  coolant,  Nb— IZr  containment,  Na  heat— pipe  radiators 
and  SiSe  thermoelectric  generators  (TE6s).  At  this  size,  any  superior — 
ity  of  the  CLCP  TEMP  will  be  most  evident  since  experience  has  shown 
that  integration  of  a  thermoelectric  generator  and  an  electromagnetic 
pump  (EMP)  into  TEMP  becomes  harder  for  larger  NSPSs. 
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Task  1. 
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thermoelectric 


Helmholtz  Magnet  TEMP  Revisited 

In  a  previous  work‘>->,  a  TEMP  using  Helmholtz  rings  (coreless)  to 
produce  the  magnetic  field  was  studied.  The  unsatisfactory  results 
obtained  led  to  the  linear  coreless  magnet  concept.  It  was  decided  to 
revisit  the  Helmholtz  TEMP  (HM-TEMP)  to  establish  a  basis  for  comparison 
of  the  relative  merits  of  the  two  coreless  concepts.  Additional 
analysis  was  required  since  the  previous  work  on  the  HM-TEMP  was 
incomplete.  Also,  analytical  expressions  have  been  developed  to 
evaluate  the  effects  of  circular  conductors  and/or  coils  whose  cross- 
large  enough  to  be  significant  in  determining  the  magnetic 
fields  produced,  and  their  distributions.  Such  expressions  have  not 
been  developed,  or  are  unknown,  for  the  parallel  counter-current  (PCC) 
magnets  to  be  used  in  the  CLCP. 

The  effects  of  size  for  the  parallel-counter-current  (PCC)'  flux 
generators  or  linear  magnets  (LMs)  is  evaluated  in  Task  6  below.  The 
inclusion  of  the  large  conductor  effects  in  the  Helmholtz  model  thus 
provides  an  easy  way  of  evaluating  these  effects  for  comparison  with  the 
same  effects  in  the  PCC  magnets.  Finally,  the  Helmholtz  analytical,  or 
design,  model  with  its  circular  geometry  also  provides  a  basis  for 
assessing  the  relative  merits  of  Mo  and  of  Zr  diboride  as  the  electrical 
conductors.  The  HM-TEMP  analytical,  or  design  models  (HMM2EMP2  and 
HMM2TEB2),  are  presented  in  Appendices  A  and  B. 

1.2  Hass  Optimization  of  TEMP  Using  the  HMM2EMP2  Model 

The  HMM2EMP2  model  was  developed  for  a  4  psi  —  lOO  gpm  EMP.  It 
also  estimates  the  total  mass  of  the  TEMP  assuming  a  g/W(e)  derived  from 
previous  analyses  or  from  an  earlier  analysis  of  the  TEG  for  this 

The  mass  of  the  TEG  is  estimated  by  multiplying  the 
W(e)  needed  by  the  EMP  by  the  TEG  g/W(e).  A  multivariate  optimization 
was  carried  out  where  the  criterion  was  the  lowest  possible  mass 
estimated  (EMASS)  for  the  TEMP,  which  includes  both  the  EMP  and  the  TEG. 
The  optimization  variables  were  H,  the  height  of  the  pump  duct,  LT,  the 
length  of  the  pumping  section,  B,  the  axial-magnetic-flux  density  in  the 
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pump  duct,  and  ETA,  the  length  of  one  side  of  the  cross— section  of  the 
Helmholtz  conductor/coil.  A  typical  optimization  curve  is  shown  in 
Figure  Tl.l,  where  the  estimated  TEMP  mass  is  plotted  against  the  square 
cross-section  dimension  of  the  coil.  The  results  of  the  optimization  is 
an  estimate  of  the  HM-TEMP  mass  of  217  kg  of  which  the  mass  of  the  TEG 
is  180  kg  or  83%  of  the  total.  That  leaves  37  kg  for  the  EMP  of  which 
the  Helmholtz  coils  are  30  kg.  It  is  obvious  that  the  estimated  TEMP 
mass  is  only  as  good  as  the  specific  mass  (45  g/W)  that  was  assumed  for 
this  iteration  of  the  EMP/TEG  analysis  cycle.  When  the  EMP  is 
electrically,  thermally,  and  hydraul ical ly  matched  to  the  TEG,  EMP  mass 
usually  increases.  The  g/W(e)  of  the  EMP  derived  from  this  analysis, 
along  with  the  power  and  voltage  requirements  of  the  EMP,  must  be 
matched  in  the  optimization  of  the  TEMP  using  the  HMM2TEG2  Model. 

In  fact,  it  was  found  that  TEG  could  not  be  matched  to  the  HMM2EMP2 
in  the  highly  integrated  configuration  that  was  applicable  with 
permanent  magnet-  (PM)  and  electro  magnet-  (EM)  TEMPs.  In  order  to  get 
enough  heat  transfer  surface  to  supply  the  heat  rate  needed  to  produce 
the  required  power,  the  length  and  height  of  the  pumping  section  had  to 
be  increased  significantly.  This,  in  turn,  degraded  the  efficiency  of 
the  EMP  and  required  still  more  power  to  be  generated  by  the  TEG.  The 
system  was  obviously  divergent  with  no  practical  solution  attainable 
without  changing  the  configuration.  The  solution  is  to  make  the  TEG  a 
separate  component  with  its  own  heat  exchanger.  However,  the  power  and 
mass  of  the  TEG  could  be  estimated  from  experience  with  the  specific 
masses  of  TEGs  that  have  been  designed  for  space  power  applications. 
Table  Tl.l  is  a  comparison  of  three  types  of  TEMPs  to  provide  some 
perspective  of  the  merits  of  the  HM-TEMP  and  of  the  LM-TEMPs  to  be 
developed  in  the  remainder  of  this  report.  The  data  for  the  permanent- 
and  electro—  magnet  TEMPs  (PMM2  and  EMM2)  were  taken  from*^*,  while  that 
for  the  Helmholtz  magnet  TEMP  (HMM2)  was  developed  as  described  above. 

It  is  obvious  that  TEG  becomes  increasingly  the  dominant  mass  of 
the  TEMPs  as  the  power  required  to  produce  the  magnetic  field  increases. 
It  is  also  known  that  the  TEG's  waste  heat  radiator  is  the  dominant  mass 
of  the  TEG,  If  the  coreless  EMP  and  TEG  can  be  operated  at  higher 
temperatures  than  the  TEMPs  based  on  ferromagnetic  materials,  then  they 
may  become  competitive  in  mass  per  delivered  pumping  power  as  well.  The 
Curie  point,  the  temperature  at  which  materials  lose  their  ferromagnetic 
properties,  is  about  1175K  for  most  of  the  high— temperature  magnet 
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materials,  such  as  Alnico  V  and  Hyperco  20.  If  a  coreless  TEMP  could  be 
operated  at  200K  higher  temperature,  i.e.,  ~1375K,  then  the  radiator 
could  be  reduced  in  mass  by  a  factor  ( 1375/1175)-4,  or  1.875.  Assuming 
the  radiator  mass  to  be  507.  of  the  180  kg  of  the  HMM2  TEMP,  or  190  kg, 
90/1.875  would  be  its  lowered  radiator  mass,  or  48  kg,  and  the  HMM2  TEMP 
would  still  be  169  kg.  A  substantial  improvement,  but  still  not 
competitive  unless  there  are  other  indirect  benefits  of  the  higher 
temperature  capability.  The  case  for  the  linear  magnet  (LM)  will  be 
made  later  in  this  report. 
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Table  Tl.l 

Comparison  of  Three  Types  of  TEMPs 


PMM2 

EMM2 

HMM2 

TEMP  Mass 

(kg) 

27.3 

70.2 

217 

TEG  Mass 

(kg) 

14.4 

43.4 

ISO 

TEG/TEMP  Mass 

(■/.) 

52.7 

61.8 

S3 

Hydraulic  Power 

Delivered  (We) 

174 

174 

174 

EMP  Efficiency 

(7.) 

33.9 

19.5 

10.2 

Total  TEMP  F'ower 

(We) 

513 

1286 

4045 

Terminal  Voltage 

of 

TEG/ EMP 

(mV) 

106 

280 

317 

Optimum  Magnetic 

Flux 

Density 

(mT) 

7B 

77 

40 

Specific  Mass  of 

TEG 

(g/We) 

28 

33,8 

45 

Specific  Mass  of 

EMP 

(g/We) 

25.1 

20.8 

9.15 

1.3  The  Core less  Linear  Conduction  Pump  (CL CPI 

An  analytical  model  was  developed  for  the  Coreless  Linear 
Conduction  Pump  (CLCP),  which  is  presented  in  Appendix  C  and  is  named 
LMM2EMP2.  LMM2  designates  a  linear  magnet  and  an  M2  pump  specification, 
i.e.,  a  pump  delivering  100  gpm  at  a  head  of  4  psi .  The  EMP2  means 
simply  the  electromagnetic  pump;  analytical  model,  version  2.  When 
coupled  with  the  LMM2TE62,  a  complete  CLCP-TEMP  analytical  model  is 
created. 

The  principal  items  in  the  development  of  the  CLCP  are  the 
parallel,  counter-current  conductors  that  constitute  the  magnet  which 
develops  the  magnetic  field  to  interact  with  the  armature  current  in  the 
pump  duct  to  produce  the  pumping  action.  This  item  is  called  the  linear 
magnet,  LM,  throughout  this  document.  Throughout  the  analysis,  SI  units 
are  used  exclusively.  Only  the  pump  specification  (gpm  and  psi)  is 
stated  in  non-SI  units.  However,  in  the  models,  the  specifications  were 
converted  to  SI  units  for  all  subsequent  use. 

The  use  of  vector  analysis  was  considered,  as  is  frequently  done  in 
magnetic  field  analysis.  However,  the  size  and  directions  of  the 
vectors  must  be  accounted  for  separately  in  the  computer  modeling. 
Thus,  vector  analysis  did  not  seem  to  provide  any  advantage  in  the 
computer  model  over  computing  the  magnitude  in  the  usual  manner  and 
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using  a  geometrical  analysis  to  provide  the  directional  aspects. 
Therefore,  the  latter  methodology  seemed  more  natural  for  computer  use, 
and  thus,  it  has  been  used  here. 


1.4  Optimum  Magnetic  Flux  Density  in  EM  Conduction  Pumps 
The  pressure  delivered  by  the  pump  to  the  external  hydraulic 
circuit  is  given  by  the  following  equation; 


PD=C1*B*I-C2*B'^2*Q-C3*B'^2*Q-C4»Q^'2  Eqn.  1 

where  PO  =  Pressure  delivered  to  the  external  hydraulic  circuit 
(Pa) 

B  =  Magnetic  field  density  (Tesla)  normal  to  the  fluid  and 
the  armature  current  flow  directions 

I  =  Electric  armature  current  (A) 

Q  =  Volumetric  flow  rate  of  the  fluid  pumped  (m''3/s) 

Cx  =  Coefficents  relating  to  fixed  parameters  (temperatures, 
materials,  and  pumping  section  dimensions)  defined  by 
the  detailed  terms  of  the  equation  and  subscripts  x=l 
through  4 — identifying  the  coefficient  of  each  term. 

The  terms  on  the  right  side  of  Eqn.  1  are,  respectively; 

Term  1  -  the  total  pressure  rise  developed  in  the  pump 
Term  2  —  loss  of  pressure  from  the  braking  effect  developed 
by  the  counter— current  across  the  duct,  resulting 
from  the  counter  or  back  EMF  generated  by  the  fluid 
flowing  across  the  magnetic  field 
Term  3  —  pressure  loss  from  eddy  current  and  fringe  end 
effects  that  develop  near  the  entrance  and  exit  to 
the  pumping  section 

Term  4  —  Internal  pressure  loss  from  friction  and  from 
contraction  and  expansion  in  the  nozzle  and  diffuser 
at  the  ends  of  the  pumping  section. 

Taking  the  derivative  of  Eqn.  1  with  respect  to  B  and  setting  it 
equal  to  zero,  the  optimum  value  of  B  is  found  to  be; 

B(opt)  =  I/Q*Cl/(2!)t(C2+C3)  )  Eqn.  2 

Thus,  the  optimum  magnetic  flux  density  is  seen  to  be  proportional 
to  the  ratio  of  the  armature  current  to  the  flow  rate  where  the 
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proportional i ty  coefficient  is  the  constant  ratio  of  Cl  to  2j|c(C2+C3). 
Unfortunately,  this  optimum  cannot  be  achieved  in  TEMPs  which  are 
optimised  for  the  minimum  mass  per  unit  of  pumping  power.  Minimum  mass 
is  reached  by  reducing  the  power  that  must  be  provided  by  the  TEG,  since 
the  specific  power  (g/We)  to  be  supplied  by  the  TEG  is  greater  than  the 
specific  power  (g/We)  demanded  by  the  EMP.  While  the  power  needed  by  a 
permanent  magnet  (PM)  TEMP  is  only  the  armature  load,  in  electromagnet 
(EM)  TEMPs,  HM,  and  LM  TEMPs,  the  magnet  must  be  powered  in  addition  to 
the  armature  power.  Consequently,  the  optimum  magnetic  flux  density  and 
maximum  EMP  efficiencies  are  not  realized. 

Partitioning  of  Current  for  Armature  and  Linear  Maonet  Use 

When  considering  CLCP ,  the  amount  of  current  available  will  be 
limited.  This  raises  the  question  of  how  the  current  should  be 
partitioned  between  the  armature  current,  lA,  and  the  field  current,  IF, 
in  the  pump  as  a  whole  and  everywhere  along  the  pumping  section  length. 
The  pressure  rise  developed  by  the  pump  is: 

PO  =  B»IA*W/(W*H)  =  B/H»IA 

Now  if  we  let  lA  =  I*f  where  I  is  the  total  pump  current  and  f  is 
the  fraction  of  I  that  is  armature  current,  then 

PO  =  B/H»f»I  but  B  =  MuO»(l-f )»I/(pi»d) 

where  (1-f)  =  the  fraction  of  I  that  is  IF,  the  field  current. 

Then , 

PO  =  MuO/(PI*d»H)*I»f*(i-f )  =  CFl*(f-f''2)  Eqn.  3 

Differentiating  Eqn.  3,  setting  the  result  equal  to  zero  and  solving 
for  f  produces  the  result: 

d(PO)/df  =  CF*(l-2»f)  =  O 


therefore:  f  =  1/2 

Thus,  as  far  as  developing  the  highest  possible  head  in  the  pump 
from  a  given  supply  of  current,  I A  =  IF.  Of  course,  adding  the  other 
terms  to  the  equation  for  the  delivered  head  and  flow  rate  specified, 
Eqn.  1,  the  result  is  modified  as  shown  here 
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f 


(CF1+2»(CF2+CF3) )/(2»(CFl+CF2+CF3> ) 


Eqn .  4 


where  CFl,  CF2,  and  CF3  are  the  coefficientB  of  the  first  three  terms  on 
the  right  side  of  Eqn.  1,  respectively.  The  more  complete  equation,  4, 
for  the  optimum  partition  of  the  current  is  greater  than  1/2  but 
probably  not  very  much.  This  being  the  case,  then  previous  relationship 
for  B(opt),  Eqn.  2,  which  provides  maximum  pumping  efficiency,  cannot 
be  simultaneously  attained  by  this  arrangement.  It  should  be  noted  that 
an  electrically  parallel  arrangement  between  the  armature  and  the  field 
circuits  is  inherent  in  the  counter  current  arrangement. 

In  the  CLCP  concept,  part  of  the  field  current  would  later  (further 
down  the  pumping  section)  be  converted  to  armature  current  until  it  was 
all  consumed.  With  the  partitioning  fraction  of  the  total  current  close 
to  1/2,  half  of  the  field  current  would  go  to  armature  current  and  the 
other  half  to  residual  magnetic  field  current  for  each  subsequent 
partition.  In  a  preliminary  analysis  of  bow  the  current  would  be 

distributed,  the  CLCP  was  divided  into  5  serial ly— connected  pumping 
cells  to  be  supplied  by  common  armature  and  field  current  busbars  that 

also  consisted  of  a  linear  magnet.  When  power  was  supplied  to  the 

entrance  end  of  five  pump  cells,  and  the  current  partitioned 

half-and-half  to  the  armature  and  field  currents,  that  left  half  the 
current  for  the  second  cell,  l/4th  for  the  3rd  cell,  l/8th  for  the  4th 
cell,  and  l/16th  for  the  5th  cell.  But  the  pumping  forces,  or  pressure 
developed,  are  proportional  to  the  product  of  the  armature  and  field 
currents,  or,  since  IA=IF  in  each  cell,  the  square  of  the  current  to 
each  cell.  Thus,  the  following  proportional  pumping  powers  and  ratios 
would  be  produced  in  the  five  cells: 

Proportional  Pumping 
Pumping  Power 

Power  Ratio 


1st 

cell 

(1/2)= 

1/4 

1 

2nd 

cell 

(1/4)= 

=: 

1/16 

1/4 

3rd 

cell 

(1/S)= 

= 

1/64 

1/16 

4th 

cel  1 

(1/16)= 

= 

1/256 

1/64 

5th 

cel  1 

(1/32)= 

= 

1/1024 

1/256 

in  a  decreasing  geometric  progression. 
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Thus  the  optimum  partitioning  indicated  that  nearly  all  the  pumping 
would  be  provided  by  the  first  one  or  two  cells,  but  the  mass  of  the 
other  cells  would  be  undiminished,  or  nearly  so.  It  is  obviously  more 
effective  to  eliminate  the  last  3  or  4  cells  under  this  partitioning . 
The  optimum  partition  was  abandoned  and  smaller  partitioning  fractions 
considered.  But  if  the  partition  fraction  were  made  still  smaller,  then 
the  progression  would  be  less  steep  and  more  cells  would  be  needed.  But 
ultimately,  the  end— most  cells  would  be  relatively  very  unproductive  and 
would  be  eliminated.  After  cell  truncation,  the  last  cell  would  pump  a 
significant  fraction  of  the  first  cell.  Why,  then,  integrate  the 
currents  in  the  first  place?  The  whole  idea  of  using  the  current 
twice  first  to  provide  the  magnetic  field  and  then  to  become  armature 
current  did  not  seem  to  make  much  sense,  or  provide  much  gain,  if  any, 
over  separate  armature  and  field  circuits.  In  addition,  it  was  found 
thdt  as  voltage  between  the  two  counter— current  busbars  dropped 

theoretically  to  zero  at  the  pump  exit,  the  voltage  to  drive  the  current 

'tbe  duct  was  competing  with  the  back  EMF  generated  by  the  fluid 
flowing  through  the  magnetic  field.  For  this  reason,  some  end  cells  had 
to  be  cut  off  but  not  the  busbars  so  that  the  voltage  at  the  last  cell 
was  sufficient  to  overcome  the  back  EMF.  For  these  reasons,  and  the 

elimination  of  the  complexities  in  controlling  and  partitioning  the 

currents,  the  integrated  armature  and  field  current  concept  was 
eliminated  in  favor  of  separate  armature  and  field  circuits.  What 
remained  was  a  simpler  concept  with  the  basic  essence  of  the  CLCP; 
namely,  an  EMP  with  a  linear  magnet. 

1.4.1  The  Infinite  Linear  Magnet 

is  a  magnetic  field  associated  with  every  electric  current 
such  that,  according  to  the  right  hand  rule,  there  is  a  circular 
magnetic  field  force  encircling  a  conductor  as  the  direction  taken  by 
the  fingers  of  the  right  hand  grasping  the  conductor  with  its  thumb 
pointing  in  the  direction  of  conventional  current  flow.  The  field  is  a 
vector  quantity  which  is  perpendicular  to  both  the  current  and  a  normal 
radial  vector  from  the  conductor  to  a  point  of  interest  (see  Figure  Tl. 
2). 
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Figure  T1.2.  Ampere's  Law  as  Applied 
to  a  Long  Straight  Conductor 

The  magnetic  field,  density  B,  associated  with  a  current-carrying 
conductor  of  negligible  cross-section  and  infinite  length  is  given  by 
Ampere's  Law**»®*  as 


dH  =  I*dL»5in(phi)/(4»pi»d^2) 


where  H  =  magnetic  field  intensity  (A/m) 

I  =  current  (A) 

dL  =  differential  length  of  conductor  (m) 
phi  =  angle  made  with  dl  when  viewed  from  the  point 
interest  in  the  field 

and  d  =  distance  from  dL  to  the  point  of  interest  (m) 

If  we  integrate  the  above  equation  from  minus  infinity  to 
infinity,  the  result  is  that 


p-rr/?. 

H  =  I/(4*pi»d)»  \  COS (phi)d(phi) ) 


of 


+ 
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Since  the  integral  is  exactly  2, 


H  =  I/2*pi*d 


(A/m) 


Then  the  magnetic  flux  density  is  B  =  muD»mur*H 

where  B  =  magnetic  flux  density  (T) 

muo  =  absolute  permeability  (V-s/A-m)  or  (T-m/A) 
mur  =  relative  permeability 
muo  =  4*pi*10''-7. 

Since  mur  =  1  for  vacuum  and  air  and  for  all  but  ferromagnetic 


materials. 


B  =  muo»I/(2*pi*d)  (T) 


Eqn .  5 


For  two  infinite  parallel  counter-current  conductors,  separated  by 
a  distance  2#d,  the  flux  density  at  the  mid-separation  point,  d  is  2*B 
or 

B(mid-point)  =  muD»I/(pi»d) 

and  if  d  is  the  separation  distance  between  conductors 

E(mid-point)  =  2*muo»I/ ( pi*d ) 

1.4.2  The  Finite  Linear  Magnet 

For  a  finite  (in  length)  straight  current  conductor  (see  Figure 
Tl.o),  the  magnetic  field  intensity  at  a  point,  P,  separated  from  the 
conductor  by  d  meters  is  given  by  the  Biot-Savart  Law*'^*, 

H  =  I/(4*pi»d)»Csin{bl)+sin(b2) 3 


where  bl  = 


and  b2 


—  angle  between  a  line  from  P  to  one  end  of  the  conductor 
^od  the  normal  from  the  conductor  to  P,  or 

=  atan  (Ll/d) 

angle  between  a  line  from  P  to  the  other  end  of  the 
conductor  and  the  normal  from  P  to  the  conductor,  or 

—  atan(L2/d)  where  LI  and  L2  are  the  lengths  of  conductor 
in  each  direction  from  the  point  of  interest 


15 


Figure  T1.3.  Biot— Savart  Law  Diagram 

For  a  pair  of  finite  counter — current  conductors  of  the  same  length 
and  ends  not  displaced  from  one  another,  the  magnetic  intensity  at  the 
mid-point  between  the  conductors  will  be  just  twice  that  of  a  single 
conductor.  If  w  (for  width)  is  the  length  of  the  line  between 
conductors,  then  at  any  point,  P,  between  conductors,  the  magnetic 
intensity  will  be  the  sum  of  the  contributions  from  each  of  the  two 
conductors-  At  any  point,  P,  between  the  conductors,  the  distance  from 
the  first  conductor  is,  say  fw»w,  where  fw  is  simply  a  fraction  of  w, 
then  the  distance  from  the  second  conductor  will  be  (1— fw)»w.  Then 

H  =  Cl/(fw*w)+i/((i-fw)*w)3*I/4*pi»Csin(bl)+sin(b2)3 

The  magnetic  flux  density,  B,  is  then  simply  muo»H-  The  above 
equation  is  the  basis  for  the  linear  magnet  (LM)  and  was  written  into 
the  LMM2EMP2  model  program. 

Using  the  LMM2EMP2  LM  model,  a  typical  magnetic  field  distribution 
between  the  magnetic  field  busbars  was  calculated  as  shown  in  Figure 
T1.4  at  the  mid-length  of  the  magnet.  It  is  seen  that  the  field  is 
strongest  near  the  conductors  and  is  minimum  at  the  mid— point  between 
the  busbars.  In  Figure  T1.5,  the  magnetic  field  distribution  along  the 
pump  duct  between  the  busbars  is  shown.  The  busbars  coincide  in  length 
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with  the  pumping  section.  The  field  is  seen  to  be  strong  in  the  pumping 
section,  which  extends  from  0-2  to  0.8  of  the  total  length.  Upstream 
and  downstream  of  the  pumping  region,  the  field  is  seen  to  fall  off 
rapidly. 
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Figure  T1.4.  Magnetic  Flux  Density  Distribution  Across  the  Pump  Duct 


In  the  pumping  section,  only  the  axial  magnetic  flux  density  is 
effective  in  developing  useful  pumping  action.  The  magnetic  flux  on  the 
line  between  the  conductors  is  theoretically  all  axial  flux.  However, 
parts  of  the  pump  duct  lie  above  and  below  the  theoretical  line,  and  the 
busbars  are  not  point  current  conductors.  Thus,  it  is  necessary  to 
consider  second  order  effects  which  may  be  important.  The  magnetic 
intensity  above  the  line  between  the  conductors  will  be  considered  in 
the  next  paragraph.  The  effects  of  busbar  size  will  be  considered  in 
Task  6. 

Above  and  below  the  line  connecting  the  centers  of  the  busbars  (the 
busbar  line) ,  the  points  along  the  perpendicular  bisector  of  the  busbar 
line  are  by  definition  always  equidistant  from  the  busbars.  Thus  the 
magnetic  intensity  is  twice  that  of  one  busbar.  However,  since  the 
magnetic  vectors  are  normal  to  the  radial  vectors  from  the  busbars  to 
the  point  of  interest,  P,  they  are  not  axial  but  have  equal  opposing 
non- axial  components  that  cancel  each  other.  Their  axial  component  is, 
however,  twice  that  of  one  busbar.  The  components  are  shown  in  Figure 
T1.6  to  be  the  total  flux  times  the  cos{atan(2*h/w) ) ,  and  the  distcince 
from  each  busbar  is  d=(h''2+(w/2)''2)'‘(l/2) .  The  distribution  of  the 
total  and  axial  magnetic  flux  density  as  calculated  by  the  LMM2EMP2 
model  is  shown  in  Figure  T1.7. 


Figure  T1.6.  Vector  Diagram  of  Bbx  On  The  Bisector 
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Figure  T1.7,  Axial  Magnetic  Flux  Density  vs  Height 


When  the  point  of  interest  is  not  on  perpendicular  bisector  of  the 
busbar  line,  then  contributions  of  each  of  the  busbars  to  the  magnetic 
intensity  roust  be  calculated  separately  and  vectors  summed,  as  shown  in 
Figure  T1.8.  Thus,  in  general, 

d=(h''2+(fw!tcw)~2)^(l/2)  and  d'  =  (h'‘2+((l-fw)*w)''2)''(l/2) 

HI  =  I/4*pi*d*[sin(bl)+sin(b2)  3*cos(atan(h/(fw!»:w)  )  )  Eqn.  10 

H2  =  I/4*pi»d'*Csin(bl)+sin(b2)3*cos(atan(h/((l-fw)*w)))  Eqn.  11 

and  Bax  =  rouo^H  =  muo*(Hl  +  H2) 

®AAX 


P 

B 


Figure  T1.8.  General  Vector  Diagram  of  Bax  Between  Busbars 

1.4.4  Linear  Magnet  fLM)  TEMP 

An  analytical  model  of  the  LM-TEMP  was  developed  for  comparing  this 
concept  with  permanent  magnet,  electromagnet,  and  Helmholtz  magnet 
TEMPs.  The  LM-TEMP  model  is  essentially  the  same  model  as  the  HM-TEMP 
except  for  the  substitution  of  the  LM  model  for  the  HM  model.  The 
LM-TEMP  model  is  shown  in  Appendix  C  and  is  named  LMM2EMP2.TK. 

The  basic  configuration  of  the  LM-TEMP  is  shown  in  Figure  T1.9  and 
is  similar  to  those  used  for  evaluating  the  cored-magnet  and  HM-TEMPs. 
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The  EMP  consists  of  a  double  pass  arrangement  where  fluid  is  being 
pumped  in  opposite  directions,  as,  for  example,  to  and  from  a  heat 
source.  There  are  two  LMs,  one  just  above  and  outside  the  upper  pump 
duct,  and  the  other  directly  just  below  the  lower  pump  duct. 


HEAT  PIPE  -  FIN 


The  counter-current  busbars  in  each  magnet  carry  current  in  the 
same  direction  on  the  same  side  of  the  stacked  pump  ducts,  but  with 
cross-bars  at  opposite  ends  of  the  main  busbars.  Field  current  enters 
the  top  magnet  at  one  end  of  a  busbar  and  the  bottom  magnet  at  the 
opposite  end  of  the  busbar  on  other  side  of  the  duct.  Opposing  flow  of 
current  in  the  two  ducts  provides  self-compensation  of  the  magnetic  end 
that  would  occur  if  no  compensation  were  provided. 

The  LM  model  used  for  the  evaluation  uses  the  average  method 
described  above  in  Section  6.3.3.  While  this  model  is  not  as  accurate 
as  the  "Differential  Busbar  Model",  it  is  good  enough  for  preliminary 
comparisons,  and  much  easier  to  use.  In  fact,  the  difference  in  axial 
flux  calculated  by  the  two  methods  is  less,  the  smaller  the  conductors 
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are.  The  conductors  in  the  optimized  TEMP  calculation  are  <36  cm==  in 
cross-sectional  area,  and  most  of  them  are  less  than  20  cm=®.  The 
differences  in  the  two  methods  are  thus  <17  and  11%,  respectively, 
derived  from  Figure  T6.7. 

Some  of  the  results  of  the  calculations  are  shown  in  Table  Tl.l. 
The  results  are  similar  to  those  obtained  from  the  HM.  The  power 
required  by  coreless  pumps  of  a  given  pumping  capacity  is  much  greater 
than  required  by  the  cored  pumps,  i.e.,  the  permanent  and  electromagnet 
pumps.  Large  currents  are  required  to  produce  very  modest  magnetic  flux 
densities.  Whether  the  required  current  is  provided  as  a  single 
Amp-turn  or  many  Amp-turns,  the  total  current  must  be  the  same.  Thus, 
the  power  required  can  only  be  reduced  by  the  use  of  lower  resistivity 
materials  and  smaller  structures  (which  can  only  be  achieved  by  the  use 
of  lower  resistivity  materials).  When  more  than  one  amp-turn  is  used, 
electrical  insulation  between  turns  is  required.  Some  dilution  of  the 
current  density  by  the  insulation  is  unavoidable,  resulting  in  the  need 
for  more  amp-turns  per  unit  of  magnetic  flux  at  the  point  or  region  of 
interest.  The  use  of  insulation  at  very  high  temperatures  introduces 
additional  problems  which  would  have  to  be  addressed. 

The  magnet  power  is  critical.  Every  Watt  of  electrical  power 
required  must  be  supplied  by  the  thermoelectric  generator.  It  has  been 
found  that  in  the  coreless  magnet  TEMPs,  the  magnet  power  dominates  the 
power  required  by  the  EMP.  This,  in  turn,  means  that  the  TEG  mass  to 
supply  the  power  will  be  predominately  that  needed  by  the  LM. 

The  larger  the  TEG  becomes,  the  greater  the  thermal  energy  it  must 
be  provided  from  the  liquid  stream  being  pumped.  In  attempting  to 
integrate  the  TEG  with  the  EMP,  as  was  done  with  the  cored  TEMPs,  it  was 
found  that  heat  exchanger  requirements  distorted  the  EMP,  thereby 
reducing  its  efficiency.  The  loss  of  efficiency  is  evidenced  by  both 
armature  current  and  magnet  current  increases.  The  result  is  an  even 
larger  TEG  is  required,  etc.  That  is,  the  TEMP  design  diverges,  blows 
up,  or  there  is  no  solution.  A  configuration  change  is  required  which 
separates  the  design  into  a  thermoelectric  generator  component  with  its 
own  heat  exchanger  and  an  electromagnetic  pump  without  a  heat  exchanger. 
It  is  no  longer  a  TEMP. 
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Table  T1.2 

Comparison  of  TEMPs  Based  on  Coreless 
Helmholtz  and  Linear  Magnets 


PARAMETERS 

HM 

LM 

Axial  Magnetic  Flux 

(mT) 

40 

24 

Hass  of  Magnets 

(  kg ) 

30.8 

19.6 

Mass  of  EMP 

(kg) 

33.5 

21.3 

Mass  of  TEG 

(kg) 

227 

199 

Estimated  TEMP  Mass 

(kg) 

262 

220 

Power  Output  of  Pump 

(Wh) 

174 

174 

Hydraulic/Electric  Effic 

iency  (%) 

9.0 

7.7 

Power  Input  to  EMP 

(We) 

1935 

2165 

Power  Input  to  Magnets 

(We) 

3159 

2261 

Power  Output  of  TEG 

(We) 

5093 

4426 

Terminal  Voltage  of  EMP 

(mV) 

123 

148 

Terminal  Voltage  of  Magnets  (mV) 

251 

570 

Number  of  Pump  Cells  or 

Passes 

2 

2 

Number  of  Magnets 

2 

2 

Specific  Mass  of  EMP 

(g/We) 

6.97 

4.8 

Specific  Mass  of  TEG 

(g/We) 

44.5 

45 

Armature  Current 

(kA) 

15.7 

15.3 

Field  Current 

(kA) 

6.3 

1.9 

1.4.5  The  Feasibility  of  Core less  TEMPs 

Both  the  careless  HM  and  LM  TEMPs  are  not  feasible  except  perhaps 
in  very  small  pumps  that  are  of  no  interest  to  SDI.  While  coreless  EMPs 
are  feasible  as  separate  components,  they  may  even  be  practical  and 
useful  in  SDI  NSPSs.  When  they  are  integrated  into  a  single  component  (a 
TEMP)  where  they  must  provide  heat-exchanger  surface  and  perhaps 
current,  voltage,  and  flow  compromises  with  the  TEG,  they  have  been 
found  not  to  be  feasible.  The  hoped-for  advantages  of  the  coreless 
TEMPs  (elimination  of  the  Curie  point  loss  of  ferromagnetism  and  its 
derivative  benefits)  may  yet  be  availabe  to  SDI  with  a  pumping  system 
employing  coreless  EMPs,  heat  exchangers,  and  TEGs  as  separate 
components  working  together  to  provide  all  the  functions  of  TEMPs.  The 
investigation  of  such  systems  appears  to  have  merit. 
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No  analysis  of  larger  coreless  LM  TEMPs  for  larger  NSPSs  was 
undertaken  as  originally  planned  since  magnet  power  requirements  will  be 
more  than  proportionately  larger  because  the  volumes  to  be  magnetised 
larger,  and  the  magnets  themselves  will  be  further  removed  from  the 
location  where  the  field  must  act.  This  effect  was  observed  in  the 
previous  study*  >•>  even  with  the  PM  and  EM  TEMPs.  Thus,  the  larger  the 
coreless  TEMP  required,  the  further  from  feasibility  it  will  be. 

Task  2.  Based  on  the  results  of  Task  1 ,  a  Corel ess  Linear  Conduction 
Pump  will  be  fabricated  in  such  a  manner  that  it  can  be  tested 
in  existing  test  facilities.  A  test  plan  will  be  generated . 

2.1  Design  of  the  Test  Pump. 

Initially  it  was  confidently  thought  that  the  test  pump  could 
be  designed  to  fit  the  available  existing  loop  with  minor 
modifications.  This  did  not  turn  out  to  be  the  case. 

An  analytical  model  of  the  VCLCP  was  developed  to  design  the 
test  device  as  shown  in  Appendix  E.  The  minimum  pump  performance 
that  could  be  acccomodated  by  the  loop  was  thought  to  be  5  psi  of 
delivered  head  at  5  gallons  per  minute.  Attempts  to  design  the 
test  pump  resulted  in  very  high  magnet  currents  at  very  low 
voltages.  Test  pump  performance  was  reduced  until  a  pump  of  .2  psi 
and  1  gpm  was  reached  requiring  currents  of  12  to  30  kA  at  voltages 
^1  V.  This  would  have  resulted  in  severe  instrumentation  and  heat 
removal  problems  requiring  significant  cost  and  time  to  solve. 
Later ,  it  became  obvious  that  such  a  test  would  not  provide  any 
credible  data  for  the  CLCP.  Construction  and  testing  of  the  VCLCP 
was  deleleted  from  the  program. 

The  electrical  network  schematic  for  an  analytical  model  of 
the  verification  CLCP  is  shown  in  Figure  T2.1.  The  electrical 
network  was  developed  to  control  the  partitioning  of  the  armature 
and  field  currents  along  the  length  of  the  pumping  section.  The 
network  model  is  shown  in  Appendix  F.  This  model  was  subsequently 
incorporated  into  the  overall  VCLCP  design  model.  After  the 
decision  to  completely  separate  armature  and  field  currents,  the 
model  was  not  needed. 
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Cu  Armature 


Figure  T2,l,  Preliminary  VCLCP  Electrical  Network 


was  selected  as 


NaK-78,  the  near  eutectic  alloy  of  Na  and  K, 
the  liquid  metal  to  be  pumped  by  the  VCLCP.  Its  physical  properties 
were  taken  from  the  literature*®-^*  and  the  data  as  a  function  of 
temperature  over  the  range  of  interest  was  fitted  by  polynomial 
equations  that  could  be  readily  used  in  the  design  model.  Plots  of 
the  data,  the  fit  of  the  polynomial  equations  and  the  polynomial 
coefficients  are  given  in  a  series  of  figures  and  computer  output 
that  comprises  Appendix  B.  The  resulting  polynomial  equations  were 
incorporated  into  the  VCLCP  design  and  analysis  model.  A  similar 

set  of  equations  was  developed  for  NaK-56  alloy  before  the  NaK-78 

alloy  was  selected. 

The  design  of  the  VCLCP  never  progressed  to  the  point  where 
stress  analysis  was  appropriate  and  a  LM  cooling  system  design  was 
required . 

2.2  Fabrication  of  the  VCLCP 

The  most  critical  item  in  the  construction  of  the  VCLCP  was 
believed  (based  on  the  history  of  fabrication  dc  EMPs)  to  be  the 
attachment  of  the  armature  current  eletrodes  to  the  pumping  section 
wall.  For  the  test  pump  to  operate  at  relatively  low  temperature, 
304L  stainless  steel  and  Cu  were  selected  for  the  pumping  section 
and  eletrodes  respectively.  This  suggested  that  commercial 

fabrication  of  these  materials  for  kitchenware  might  be  applied. 
An  order  was  placed  with  Texas  Instruments'  (TI)  Foil  Division  for 
sheets  of  SS  rolled  and  bonded  together  with  Cu  to  form  a  finished 
foil  of  SS  clad  on  Cu  with  layers  of  20  and  99  mils  respectively. 
This  was  the  thickest  layer  of  Cu  that  TI  felt  they  could 

fabricate.  Additional  Cu  would  have  been  electroplated  onto  the 
electrodes  as  required.  Figure  T2.2  is  a  photomicrograph  of  the 
bonded  region  of  the  foil  and  was  used  to  verify  the  thicknesses. 
Figure  T2.3  is  a  photograph  of  a  fabrication  test  specimen  made 
from  the  foil  by  first  removing  Cu  from  selected  areas  by  precision 
machining  and  lapping,  bending  the  SS  to  the  pumping  section 
cross-section  configuration  and  then  seam  welding.  No  further 
fabrication  testing  was  conducted  after  the  decision  was  made  to 
forego  building  and  testing  of  the  pump. 

Some  electrical  insulation  fabrication  tests  were  started  but 
immediately  abandoned  upon  the  decision  to  forego  construction. 
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Task  3.  Existing  test  loops  will  be  modified  and  instrumented  to 
accommodate  the  Coreless  Pump .  and  to  collect  the  necessary 
test  data  (current  voltage,  flowrates,  pressure.  and  magnetic 
f lux) . 

3.1  The  NaK  Cart  Loop  at  ETEC 

The  NaK  Cart  Loop,  Figure  T3.1,  at  Energy  Technology  Engineering 
Center  (ETEC)  was  selected  for  testing  the  VCLCP.  However,  the  results 
of  Tasks  1  and  2  did  not  support  testing.  Modifications  to  the 
facility,  instrumentation,  power  supply  and  heat  sink  could  have  been 
made,  but  their  cost  was  beyond  the  scope  of  this  program.  This  task 
was  deleted  by  contract  modification.  However  there  are  a  few 
activities  that  were  accomplished  that  are  reported  below. 
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Figure  T3,l.  ETEC  NaK  Cart  Loop  Schematic 


An  analysis  was  made  of  the  pressure  losses  to  be  expected  in  the 
NaK  Cart  Loop  using  the  schematic  sketch  provided  by  ETEC.  The 
equivalent  velocity  head  loss  of  all  of  the  components  in  the  flow  path 
were  determined  as  shown  in  Table  T3.1.  The  losses  in  the  NaK  resevoir 
tank  and  the  oxygen  filter  were  shown  to  be  uncertain  or  excessive 
relative  to  the  low  delivery  pressure  to  be  developed  by  the  test  pump. 
The  solutions  agreed  to  by  ETEC  and  Hi-Z  were  modifications  to  the  loop 
piping  to  (1)  eliminate  the  static  head  rise  in  the  resevoir  by 
returning  the  NaK  to  the  bottom  of  the  tank  instead  of  the  top  and  (2) 
installing  valves  and  piping  to  by-pass  the  filter  during  testing. 


Table  T3.1 

Velocity  Head  Losses 


Velocity  Head 

Quantity  Element  Coefficient.  K 


2 

45®  El-3/4" 

0.45 

i 

Tee  (3/4"-l/2"-3/4" )  -  Side  exit 

5.0 

2 

Valves  1/2”  (fully  open) 

1.15 

1 

Flowmeter 

0. 

1 

EH  Pump  (nozzle  &  diffuser) 

0.08 

1 

Reservoir  (expansion  &  contraction) 

1.5 

1 

Cold  Trap  (mechanical  filter) 

lO 

1 

90°  El  (1/2”  X  3/4") 

0.94 

1 

90°  El  3/4” 

0.9 

1 

Tee  (1/2”)  (straight  through) 

1.8 

1 

3/4"  Tubing  (54”  total  length) 

2.52 

1 

1/2"  Tubing  (241”  total  length) 

16.87 

41.21 

Loop  pressure  drop  is  estimated  to  be  -^7  x  10“=  Pa  at  2  gpm  flowrate 
Static  head  of  reservoir  entrance  is  ^10“^  Pa 
(There  is  a  free  surface  in  the  reservoir  where 
cover  gas  pressure  may  be  applied) 

Static  head  (half— full  reservoir)  2l5  >'  10®  Pa 
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ETEC  indicated  that  the  existing  EM  flowmeter  had  never  been 
calibrated  for  such  low  flow  as  would  be  occuring  in  the  VCLCP  tests. 
Mapping  of  the  magnetic  flux  distribution  and  calibration  of  the  flow 
meter  would  have  been  required. 

The  existing  pressure  gage  (pitot-tube  type)  in  the  NaK  Cart  Loop 
had  a  range  that  was  not  compatible  with  the  pressure  head  to  be 
developed  by  the  VCLCP.  A  remotely  operated  differential  liquid— metal 
manonmeter  would  have  been  required. 

3.2  Magnetic  Flux  Measurement 

Measurement  of  the  magnetic  flux  and  its  distribution  can  be 
difficult.  Measurement  on  the  installed  pump  would  be  non-existent  to 
severely  limited  and  of  limited  usefulness.  Measurements  outside  of  the 
loop  can  only  confirm  that  the  stray  flux  is  low  and  very  local.  The 
areas  of  interest  are  in  the  pumping  section  and  near  the  ends  of  the 
pumping  section.  Measurements  at  these  locations  have  to  be  outside  of 
the  loop  and  at  reduced  power.  Two  methods  are  possible:  one  uses  ac 
power  to  the  magnet  with  induction  coils  that  measure  the  flux  induced 
voltage.  The  sensivity  is  low  at  the  low  frequency  of  the  power  supply. 
The  directional  aspects  are  difficult  to  measure  without  an  elaborate 
three  dimensional  positioning  system. 

A  second  method  uses  a  Hall  Effect  (HE)  sensor.  Such  sensors  are 
commercially  available.  Complete  HE  flux  meters  are  also  available 
commercially  but  can  be  expensive.  Some  preliminary  tests  were  made  at 
Hi-Z  using  a  HE  sensor  sample  given  by  the  Siemens  Co.  A  simple 
instrument  was  made  using  a  standard  circuit.  The  small  size  (<3  mm)  of 
the  sensor  was  a  considerable  advantage  over  the  size  of  an  effective 
induction  coil.  (It  could  easily  be  fit  into  the  pumping  section  duct). 
However,  the  HE  device  was  little  better  in  directional  aspects  than  the 
induction  coils.  Either  ac  or  dc  magnet  power  can  be  used  with  HE 
method.  Some  measurements  were  made  of  the  magnetic  flux  and  its 
vertical  distribution  between  two  parallel  counter — current  Cu  busbars. 
Using  a  welding  machine  power  supply  both  ac  and  dc  measurements  were 
made  successfully  and  correlated  well  with  the  analytical  model.  How¬ 
ever,  the  power  supply  was  limited  in  capacity  to  '^350A.  The  HE  method 
was  selected  but  no  further  work  on  it  was  conducted. 


Power  Supply 

A  power  supply  capable  of  supplying  <24  kA  current  at  <1  Vdc  was 
indicated  by  the  design  analysis.  ETEC  could  supply  50  kA,  at  10  Vdc  by 
a  rectified  30  source  without  a  filter.  The  details  had  not  been  worked 
out  at  the  time  the  test  was  cancelled,  but  long  heavy  Cu  cables  were 
foreseen  in  which  the  voltage  would  be  dropped  to  the  user's  requirement 
and  parasitic  power  in  the  cable  would  be  dissipated  by  forced  air 
cooling,  if  necessary.  To  get  to  IV,  about  90%  of  power  would  be 
wasted,  or  about  450  of  500  kW.  The  remaining  power,  25-50  kW,  would  be 
consumed  in  the  VCLCP.  That  could  not  be  accomplished  without  a  local 
water  heat  sink.  That  power  supply  was  to  power  the  VCLCP,  which  would 
produce  an  output  of  less  than  1/2  W  of  hydraulic  power?  Incredible! 

is  not  what  was  envisioned  in  our  proposal  which  said  "In  parallel 
with  analyses,  a  low-cost  proof-of-principle  coreless  pump  will  be 
designed  and  tested  in  a  small,  low-cost  test  loop". 

will.  ^  performed,  and  the  data  will  be  analyzed  to 
verify  the  design. 

This  task  was  deleted  by  program  contract  modification.  See 
Task  3  above. 

B.  jTBpart  will  written,  covering  the  design,  and  test 

results,  arid  will  specifically  address  the  feasibility  of  the 
coreless  pump  iri  thermoelectric  magnetic  pumps.  Based  on  the 
results,  a  Phase  1±  proposal  may  be  developed  which  would 
include  tjie  deve  1  opmen t  of  a  prototype  thermoelectric  magnetic 
pump. 

This  report  is  the  product  of  Task  5. 

The  use  ojp  ZrBg  as  magnetic —fie Id  conductors  and  busbars  will 
b^  assumed  and  the  affect  of  its  use  on  CLCP  f easibi 1 i tv  wi 1 1 
be  anajyzed  and  compared  with  the  normal  use  of  Mo  ^  1200  K 
or:  higher  for  these  purposes.  The  effects  of  large 

current  conductors  on  magnetic  f 1 ux  generation 
will  be  analyzed . 


Task  4. 


Task  5- 


Task  6. 


6.1  Comparison  of  ZrBz  &  Mo  as  High  Temperature  Conductors  in 
LM  Magnets 

6.1.1  High  Temperature  Conductor  Materials 

At  the  high  temperatures  of  space  power  systems,  the  choice  of 
electrical  and  thermal  conductors  is  limited.  Without  going  to  such 
exotic  concepts  as  canned  liquid  metals.  Mo  is  a  frequent  choice  and  is 
the  likely  choice  for  the  TEMPs.  However,  Hi— Z  Technology,  Inc.  (Hi— Z) 
is  currently  in  the  SBIR  Phase  II  development  of  ZrBa,  zirconium 
diboride,  which  is  comparable  in  electrical  resistivity  to  Mo  at  lOOOK 
and  can  be  used  up  to  2000K.  The  program  is  funded  by  SDID  and  managed 
by  The  Defense  Nuclear  Agency  (DNA).  Zirconium  diboride's  advantage  is 
its  density  which  is  only  6  Mg/m'''3  compared  to  10.1  of  Mo.  Based  on  the 
obvious  potential  for  lower  mass  conductors  using  zirconium  diboride, 
its  use  in  coreless  TEMPs  was  analyzed  to  determine  if  lower  pump  masses 
would  result  and  if  it  would  contribute  to  core less  pump  feasibility. 
The  following  analysis  was  conducted. 

This  analysis  was  carried  out  using  the  HMM2EMP2  analytical  model 
of  the  Helmholtz  electromagnetic  pump  which  would  deliver  100  gpm  of  Li 
coolant  at  a  head  4  psi  (See  Appendix  A).  The  masses  of  the  Helmholtz 
conductor  rings  made  of  Mo  and  ZrBs  were  compared  while  all  remaining 
pertinent  variables  were  adjusted  so  as  to  optimize  the  TEMP  design  with 
respect  to  the  total  TEMP  mass,  i.e.,  including  TEG  mass  as  described  in 
Section  1.2  above.  The  results  of  the  optimization  analyses  are  shown 
in  Table  T6.1. 

The  analytical  result  reflects  the  difference  in  their 
conductivities  per  unit  mass.  The  conductivity  per  unit  mass  of  ZrBa 
and  Mo  are  4.35  x  10®  and  3.54  x  10®  cm*/g.0hm,  respectively.  Thus,  the 
conductivity  per  unit  mass  of  ZrB^  is  23%  greater  than  that  of  Mo  at 
i200K.  The  results  show  the  estimated  mass  of  the  HM-TEMP  optimized  for 
the  use  of  ZrB^  is  22%  less  than  that  when  optimized  for  Mo  use. 
Similar  differences  appear  in  the  power  output  and  mass  of  the  TEG.  The 
superiority  of  ZrB*  with  respect  to  the  Mo  is  clear.  At  higher 
temperatures,  the  situation  would  seem  to  favor  Mo  because  of  its  lower 
coefficient  of  resistivity  compared  to  ZrBa,  as  shown  in  the  equations 
below.  However,  there  appears  to  be  some  confusion  in  the  accuracy  of 
these  equations  or  the  data  sources  from  which  they  were  derived.  This 
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Table  T6.1 

Comparison  of  HM-TE!iPs 
Using  Mo  versus  ZrBs  Conductors 


Conductors 


PARAMETERS 

Mo 

ZrB: 

Axial  Magnetic  Flux 

(mT) 

40 

24 

Mass  of  Magnets 

(kg) 

30.8 

12.0 

Mass  of  EMP 

(kg) 

33.5 

14.3 

Mass  of  TEB 

(kg) 

227 

190 

Estimated  TEMP  Mass 

(kg) 

262 

204 

Power  Output  of  Pump 

(Wh) 

174 

174 

Hydraul ic /Elec trie  Ef f icienc 

y  (7.) 

9.0 

6.8 

Power  Input  to  EMP 

(We) 

1935 

2570 

Power  Input  to  Magnets 

(We) 

3159 

1690 

Power  Output  of  TEG 

(We) 

5093 

4260 

Terminal  Voltage  of  EMP 

(mV) 

123 

393 

Terminal  Voltage  of  Magnets 

(mV) 

251 

244 

Number  of  Pump  Cells  or  Passes 

2 

2 

Number  of  Magnets 

2 

2 

Specific  Mass  of  EMP 

(g/We) 

6.97 

3.4 

Specific  Mass  of  TEG 

(g/We) 

44.5 

45 

Armature  Current 

(kA) 

15.7 

6.5 

Field  Current 

(kA) 

6.3 

3.5 

is  being  pursued  under  the  referenced  ZrB=  development  program  noted 
above  in  this  section. 


RhoZrB  —  *“5.710+0.03573|cT  micro— Ohm— cm ^ 
Rhoho  =  “0.183+0.0321:1:1  micrD-0hm-cm<®  > 


6.1.2  The  Electric  Resistivity  of  Molybdenum  Versus  Temperature 
The  electrical  resistivity  of  molybdenum  (Mo)  was  seen  to  be 
critically  important  in  the  design  comparison  with  ZrBa,  zirconium 


diboride,  as  the  conductors.  Resistivity  data  from  various  sources  lead 


to  some  confusion.  An  additional  attempt  has  been  made  using  well- 
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documented  Mo  data(7).  The  data  are  from  one  of  the  oldest  and  largest 
commercial  suppliers,  and  was  reviewed  and  selected  to  be  used  in  a 
highly  regarded  materials  handbook  of  recent  publication.  It  is 
regarded  to  be  as  well-founded  as  one  can  expect  at  this  writing. 
Several  points  were  taken  from  the  published  graph  of  the  full 
temperature  range  presented.  These  data  were  fit  with  a  fourth  order 
polynomial  to  derive  the  following  set  of  coefficients  for  the  equation. 
The  original  data  points  were  then  plotted  as  shown  in  Figure  T6 . 1 ,  with 
the  curve  (solid  line)  calculated  from  the  equation.  The  fit  appears  to 
be  excellent. 

RHOMo= A+B*T+C*T  **  2+D*T  3+E*T  “  4  (  micr oOhm-  cm  ) 

where  A  =  -8.144267 

B  =  +0.21863752 
C  =  -6.88*10^-10 
D  =  +2.3868*10''-8 
E  =  -5.28*10^-12 
and  T  =  Temperture  K) 

RH0ZrB2=-5 .71+0 . 0357*T  (microOhm-cm) 
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Figure  T6,l.  Verification  of  Mo  Resistivity  Equati 


Zirconium  diboride  is  actively  under  development  as  a  conductor 
material,  including  aerospace  applications*®*.  Its  electrical  resis¬ 
tivity,  among  other  properties,  is  somewhat  in  a  state  of  flux. 
However ,  the  resistivity  curve  given  above  is  regarded  as  a  reasonable 
working  model  until  the  work  in  progress  is  completed  in  1992. 

A  comparison  of  the  electrial  resistivities  of  Mo  and  ZrB^  over  a 
temperature  range  of  300  to  1500K  is  shown  in  Figure  T6.2.  At  room 
temperature,  the  resistivities  are  essentially  equal  but  Mo,  at  least 
initially,  has  a  smaller  slope.  Mo's  resistivity  is  35%  lower  at  1200K, 
and  the  curves  are  still  diverging. 
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Figure  T6*2.  Resistivity-Temperature  Curves  of  Mo  and  ZrB 


6.2  Effect  of  Conductor  Size  on  Magnetic  Field  Calculation  in  HM 

In  calculating  the  magnetic  field  produced  by  an  electric  current 
flowing  in  conductors,  it  is  generally  assumed  that  the  cross-sectional 
area  of  the  conductors  is  negligible  and  thus  two  amp-turns  will  produce 
twice  the  field  strength  of  one  amp-turn.  However,  when  the  conductors 
or  the  number  of  amp-turns  gets  very  large,  the  total  conductor  can  grow 
to  a  size  where  its  dimensions  can  no  longer  be  neglected.  For  circular 
conductors,  such  as  the  Helmholtz  conductors,  analytical  expressions  for 
the  magnetic  flux  have  been  developed  to  take  into  account  the  size  of 
the  conductors* .  Such  an  expression  is  contained  in  HMM2EMP2  for 
conductors  of  square  cross-sections.  However,  such  an  expression  for 
PCC  conductors  is  unknown  or,  at  least,  not  readily  available.  It  was 
thus  of  interest  to  use  the  Helmholtz  coil  model  to  determine  if  it 
would  be  an  important  effect  in  generating  magnetic  fields  for  coreless 
TEMPS. 

The  analysis  was  conducted  using  the  HMM2EMP2  program  (See  Appendix 
A)  and  decoupling  the  magnetic  flux  section  from  the  rest  of  the  program 
by  replacing  B,  the  axial  magnetic  flux  density  generated  in  the  pump 
duct,  by  BT,  the  test  magnetic  flux  density,  which  is  unrelated  to  the 
pumping  equations.  This  is  done  by  "quoting-out",  that  is,  by  placing  a 
quotation  mark  at  the  beginning  of  any  equation  that  causes  it  to  be 
ignored.  So  the  Bs  were  quoted— out  and  BT  inserted  for  the  coil  flux. 

The  coil  current  was  divided  by  the  coil  cross-sectional  area  to 
yield  a  current  density  that  could  then  be  held  constant.  Keeping  the 
current  density  constant,  increasing  the  cross-sectional  area,  is  equi¬ 
valent  to  raising  the  current  or  the  number  of  amp-turns  of  the 
conductors  and  vice-versa.  If  there  were  no  effect  of  the  conductor 
size,  the  flux  would  increase  in  proportion  to  the  the  cross-sectional 
area.  Otherwise,  the  increase  would  be  less  than  proportional.  Thus  by 
comparing  the  two  cases,  the  effectiveness  or  efficiency  of  the 
conductors  can  be  established.  The  result  of  such  an  analysis  is  shown 
in  Figure  T6.3.  The  magnetic  flux  is  shown  to  rise  in  proportion  to  the 
cross-sectional  area  along  the  upper  curve  when  the  size  of  the 
conductors  is  ignored,  and  to  fall  along  the  lower  curve  when  size  and 
shape  are  considered.  The  efficiency  of  the  conductors  is  seen  to 
decrease  as  the  size  grows.  These  results  were  obtained  when  the 
effective  radius  of  B  cm.  The  effect  of  size  becomes  less  as  the  radius 
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B  US  COIL  SIZE  (CD=CONSTANT) 


increases,  but  becomes  very  sensitive  to  coil  size  as  the  radius 
decreases.  It  is  clear  that  the  size  of  the  conductors  carrying  the 
magnetic-field-producing  current  cannot  be  neglected  in  coreless  TEMPs. 
The  effects  of  the  conductors  size  is  considered  next  for  the  LM. 

6.3  The  Effect  of  Conductor  Size  on  Maonetic  Field  in  Lli 

Analysis  in  previous  paragraphs  of  the  effect  of  the  size  of  the 
cross-sectional  area  on  the  efficiency  of  magnetic  field  production 
showed  that  dimensions  of  the  conductors  may  be  important  in  the  design 
of  coreless  electromagnetic  pumps.  While  analytical  expressions  are 
available  for  such  effects  in  circular  coils,  there  is  no  readily  avail¬ 
able  data  for  counter-current  conductors  such  as  those  of  interest  here. 
Thus,  the  following  method  was  developed  using  superposition  and 
adapting  analytical  techniques  available  in  the  literature. 

The  approach  is  to  subdivide  square  busbar  cross-sections  into  25 
smaller  squares  so  that  a  5  by  5  array  of  square  conductors  or  minibus- 
bars  make  up  the  total  busbar  of  each  of  two  matched  counter— current 
busbars  used  to  generate  the  magnetic  field  for  an  electromagnetic  pump. 
The  contribution  of  each  minibusbar-pair  to  the  magnetic  flux  density  at 
a  given  location  will  be  determined.  The  total  magnetic  intensity,  or 
flux  density,  may  be  obtained  by  adding  the  individual  contributions 
asserting  that  superposition  applies  for  this  purpose.  It  is  assumed 
that  the  cross-sectional  current  density  will  be  constant  for  each 
minibusbar  per  case,  but  may  be  different  from  case  to  case. 

The  equations  6  and  7  developed  in  Section  1.4.2  above  were  applied 
here.  Only  instead  of  obtaining  the  axial  magnetic  flux  density  at 
various  distances,  P,  from  a  busbar  pair,  we  now  calculated  the  axial 
magnetic  flux  density  at  a  point  contributed  to  by  busbar  pairs  located 
at  several  distances  from  the  point,  P.  The  equations  are  the  same,  only 
our  orientation  in  applying  them  changed. 

The  cross-sectional  area  of  each  of  a  pair  of  parallel,  large 
counter — current  busbars  comprising  a  LM  are  divided  into  a  five-by-five 
matrix  of  square  bars  (minibusbars)  as  shown  in  Figure  T6.4.  Equipoten— 
tial  electric  field  lines  are  circles  passing  through  the  busbar 
locations  like  those  shown  as  dotted  lines  in  Figure  T6.5.  The  magnetic 
field  vector  at  any  point  P  on  one  of  those  circles  will  be  directed 
toward  the  center  of  its  circle.  The  centers  of  these  equipotential 
circles  lie  on  the  perpendicular  bisector  of  the  busbar  lines  and  grow 
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larger  in  radius  as  they  fall  more  to  one  side  or  the  other  of  the 
busbar  line.  Lines  of  magnetic  flux  must  cross  the  equipotential  lines 
normally.  The  magnetic  flux  lines  are  also  circles  whose  centers  may  be 
found  as  follows.  First,  draw  a  radial  vector  from  the  center  of  the 
equipotential  circle  to  any  point,  P,  on  its  circumference.  Then  a  line 
drawn  tangent  to  the  circle  at  P  (and  perpendicular  to  the  radial 
vector)  will  intersect  the  extension  of  the  busbar  line  at  the  center  of 
the  magnetic-f lux-line  circle.  From  this  geometrical  information,  the 
direction  aspects  of  the  magnetic  flux  may  be  developed. 

6.3.1  P  On  or  Off  the  Minibusbar  Line 

Now  consider  the  busbar  pair  labeled  11-11  in  Figure  T6.4  and  a 
point  of  interest,  P,  at  the  mid-point  between  the  minibusbar  pair 
31-31.  Equations  6  and  7  developed  in  Section  1.4.2,  The  Finite  Linear 
Magnet,  for  the  general  case  of  P  locations  are  used  to  calculate 
contributions  of  the  minibusbar  pairs  to  the  axial  magnetic  flux  at  P, 
the  point  of  interest,  and  then  to  sum  the  contributions  for  the  total 
axial  magnetic  flux. 

Equations  for  the  25  minibusbar  contributions  were  incorporated 
into  the  MABFLDDT.TK  computer  program.  Appendix  D.  Equations  for  one  of 
the  five  vertical  minibusbar  pairs  as  a  group  in  which  the  distance  from 
the  busbar  to  point  of  interest  were  related  as  fractions  of  the  overall 
busbar  height,  eta.  Within  this  column,  the  busbar  line  width,  w,  was 
constant  but  was  written  as  wo  +  2*(nf +1/2) Iceta/S  where  eta  is  the 
length  of  one  side  of  the  overall  square  busbar  cross-section,  wo  is  the 
busbar  line  width  (center — to— center)  of  the  innermost  minibusbar  pairs, 
and  nf  is  the  column  number  (i.e.,  0-4).  Thus  there  are  five  sets  of 
indexed  equations,  one  for  each  of  five  busbar  pairs  in  a  column. 

It  was  then  necessary  to  generate  four  more  sets  of  the  column  sets 
properly  indexed,  so  there  is  one  column  set  for  each  of  the  five 
columns,  and  5  minibusbars  per  column  (layer).  The  effects  of  finite 
busbars  were  then  evaluated  below. 

6.2.2  Comparison  Between  Point  and  Finite  Busbar  Calculations 

A  comparison  of  the  flux  calculated  assuming  point  conductors  with 
busbars  of  significant  dimensions  was  conducted  similar  to  the  one  done 
for  the  Helmholtz  coils,  as  shown  in  Section  6.3  above.  A  typical 
result  is  shown  in  Figure  T6.6.  All  of  the  current  in  the  point  conduc¬ 
tors  was  assumed  to  flow  in  straight,  parallel  counter-current 
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POINT  US  FINITE  CROSS-SECTION  BUSBARS 
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conductors  of  negligible  cross-section  separated  by  a  constant  busbar 
line  of  4.2  cm  while  the  current  was  increased  as  shown.  The  calculated 
flux  at  the  mid-point  between  the  busbars  then  increases  in  a  straight 
line,  as  shown  by  the  upper  curve.  The  lower  curve  shows  the  axial 
magnetic  flux  density  at  the  same  point  as  calculated  by  summing  the 
contributions  of  the  twenty-five  minibusbars  pairs.  In  the  latter 
calculation  the  total  current,  which  is  the  same  as  that  used  for  the 
point  busbars,  was  distributed  among  the  minibusbars  as  described  above. 
The  current  density  of  130  A/cm'^2  was  used  and,  as  the  cross-section  was 
increased,  the  total  current  increased.  The  equations  take  into  account 
the  changing  distances  and  vector  relationships  as  documented  above. 
The  lower  curve  shows  the  dramatic  reduction  that  can  occur  as  a  result 
of  the  busbar  size.  Unfortunately,  the  lower  curve  represents  the  real 
world  and  becomes  a  significant  factor  in  the  potential  performance  of 
the  Coreless  Linear  Conduction  Pump  concept. 

6.3.3  Average  Busbar  Line  Length  Versus  Differential  Busbar 
Calculation 

Another  common  method  of  calculating  the  magnetic  flux  is  to  assume 
all  the  current  is  flowing  in  point  conductors  located  at  the  center  of 
the  cross-section  of  finite  busbars.  The  busbar  line  length  is  then  a 
mean,  or  average,  value  of  some  kind.  The  vertical  distribution  of  the 
current  is  ignored.  Figure  T6.7  shows  a  comparison  of  the  common  method 
with  the  differential  method  of  MAGFLDDT.TK.  As  the  cross-  section  area 
of  the  busbars  increase,  the  busbar  line  length  in  the  common  method 
increases  as  well  as  in  the  differential  method,  as  discribed  above.  In 
Figure  T6.7,  the  upper  curve  is  axial  flux  density  calculated  by  the 
differential  method  as  used  in  MAGFLDDT.TK,  and  the  lower  curve  is  the 
result  using  the  common  method.  At  relatively  small  cross-sections,  the 
difference  is  not  large;  but  at  larger  sizes,  the  difference  is  not 
negligible  being  *^21  7.  at  81  cm'^2  as  shown  here.  Apparently  in  the 
Helmholtz  coil  analytical  expression  there  an  effective  radius 
(analagous  to  half  the  busbar  line  length  as  used  here)  is  used,  which 
is  different  from  the  mean  radius. 

6.3.4  Laver  or  Co 1 umn  Effectiveness 

Another  valuable  result  of  the  differential  analysis  using 
MAGFLDDT.TK  is  shown  in  Figure  T6.8.  The  contribution  of  each  of  the 
five  columns  or  layers  from  the  inside  pairs  to  the  outside  pairs  is 
shown  both  in  terms  of  the  axial  magnetic  flux  density  produced  and  the 
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Figure  T6.7.  B  of  LM  by  Average  and  Differential  Methods 
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Figure  T6»8«  Contributions  of  Minibusbar  Columns  to 


percentage  of  the  flux  produced  by  each  column  or  layer.  The  upper 
curve  shows  the  percentage  contributed  by  each  layer  or  column  at  the 
points  (+)  indicating  their  busbar  line  lengths.  The  inner  column 
produces  more  than  three  times  as  much  as  the  outer  column.  Since  each 
column  contains  the  same  amount  of  material  and  mass  (weight),  it  will 
be  more  efficient  to  remove  the  two  outer  columns  and  use  a  rectangular 
*-*^°ss~5ec t ional  busbar  rather  than  a  square,  if  possible.  This  is 
easily  accommodated  in  MAGFLDDT.TK  by  going  to  a  3  X  5  matrix  by 
removing  the  unwanted  equations.  This  is  done  in  the  rule  sheet  of 
MAGFLDDT.TK  by  simply  "quoting"  them  out  (See  Appendix  D).  This  would 
allow  them  to  simply  be  reinstated  if  it  was  later  desired  to  use  the 
square  matrix  again. 

Effectiveness  of  Vertical  Distribution  of  Mini busbars 

The  effectiveness  of  the  vertical  distribution  or  rows  of  the  mini¬ 
busbars  that  make  up  the  overall  busbar  magnet  is  shown  in  Figure  T6.4. 
It  was  generated  by  first  assuming  a  point  of  interest  midway  between 
the  busbar  pairs  and  at  the  top  busbar  pair  level.  The  contribution  of 
each  of  the  busbar  levels  to  the  total  magnetic  flux  at  the  point  of 
interest  was  then  calculated.  The  point  of  interest  was  then  shifted 
down  to  the  second  level,  and  the  calculations  of  all  the  contributions 
from  the  other  busbar  levels  repeated.  And  so  on  until  five  points  of 
interest  at  all  five  busbar  levels  were  calculated. 

The  results  are  shown  in  Figure  T6.9  and  in  more  detail  in  Table 
T6.2.  It  is  seen  in  Figure  T6.9  that  the  axial  magnetic  flux  varies 
only  slightly  over  the  range.  The  same  result  can  be  seen  in  Table  T6.2 
for  all  the  other  columns  in  the  matrix.  The  variation  is  %  or  less. 

This  result  appears  to  justify  the  neglect  of  the  vertical  portion 
of  the  busbar  in  the  common  practice  approach  to  calculating  the 
magnetic  flux  discussed  above  in  Section  T6.3.4.  However,  when  the 
point  of  interest  exceeds  the  height  of  the  busbars,  the  vertical 
distribution  of  the  minibusbars  can  no  longer  be  ignored. 


49 


Figure  T6.9,  Contributions  of  Minibusbar  Rows  to 


Table  T6 . 2 . 
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RESULTS 

The  analysis  established  that  LM-TEMPs  and  HM-TEMPs  (i.e.,  TEMPs 
incorporating  coreless  pumps)  are  not  feasible  for  NSPS  of  64  kWe  and 
larger.  Thus,  the  first  major  objective  was  not  achieved.  The  second 
major  objective  of  building  and  testing  a  verification  pump,  and  most  of 
the  detailed  objectives,  depend  on  the  establishment  of  feasibility  and 
thus  became  moot.  Consequently,  Tasks  3  and  4  concern  with  the  test 
loop  modifications  and  the  verification  test  program  were  deleted  and 
Task  6  pertaining  to  the  use  of  ZrBz  conductors  and  calculation  of 
magnetic  flux  produced  from  large  busbar  and  their  effect  on  feasibility 
was  added.  The  use  of  ZrBa  conductors  did  result  in  better  performance, 
but  not  enough  better  to  make  the  coreless  magnet  TEMPs  feasible. 

The  loss  of  feasibility  of  the  TEMPs  results  from  the  large  power 
demand  of  the  coreless  magnets  which  must  be  supplied  by  the  TEG.  There 
is  not  enough  heat  transfer  surface  in  the  optimized  EMPs  to  supply  the 
TEG.  Heat-fluxes  required  are  too  high  by  two  orders  of  magnitude.  The 
specific  power  of  the  TEGs  was  assumed  to  be  45  g/We,  an  optimistic 
number  based  on  TEGs  built  for  use  in  space  applications.  (This  is 
equivalent  to  a  specific  weight  of  10  We/lb.) 
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The  heat  exchanger  problems  grow  worse  when  the  compromises 
required  to  match  EMP  and  TEG  electrical,  thermal,  and  hydraulic 
circuits.  The  EMPs  demand  even  more  power;  the  TEGs  get  larger  and 
require  more  heat  flow  thus  moving  the  TEMPs  further  away  from 
feasibility. 

The  normal,  or  average,  busbar  analytical  method  was  found  to  be 
adequate  for  this  feasibility  study,  but  the  "differential  method"  or 
something  better,  would  be  needed  for  detailed  design  and  performance. 
The  sizes  of  the  conductors  found  in  the  LM  magnets  were  not  so  large 
the  errors  of  <207.  in  the  flux  were  never  exceeded. 

The  LM-TEMPs  were  found  to  be  better  than  HM-TEMPs.  The  LM  mass 
was,  itself,  also  less  than  the  HM  mass. 

The  ZrBa  conductors  produced  lower  mass  in  the  otpimized  LM-TEMPs 
compared  to  the  HM— TEMPs,  but  not  low  enough  to  make  these  TEMPs 
feasible.  The  electrical  resistivity  of  ZrBz  used  in  the  analyses  was 
based  on  state-of-the-art  data,  which  is  quite  different  from  handbook 
values.  Still  lower  resistivity  of  ZrBs  is  expected  to  result  from 
ongoing  support  by  SDI  in  another  program*®*. 

Double-pass,  self-compensated  LM-TEMPs  were  found  to  be  about  a 
factor  of  4  less  massive  than  single-pass  uncompensated  LM-TEMPs. 

The  lack  of  feasibility  of  the  CLCP-TEMP  was  sufficient  reason  to 
delete  the  verification  pump  test  work  from  the  tasks  defined.  However, 
both  the  cost  and  time  to  conduct  the  program  would  have  been  grossly 
exceeded  because  of  unforeseen  pump  construction  and  instrumentation  and 
power  supply  problems  which  developed. 

ESTIMATE  OF  FEASIBILITY 

The  concept  of  CLCP  is  technically  feasible,  while  the  TEMPs  based 
on  the  CLCPs  are  not  feasible. 

ESTIMATE  OF  ECONOMIC  FEASIBILITY 

The  economic  feasibility  of  the  CLCP  is  a  moot  consideration  since 
the  CLCP-based  TEMPs  are  not  feasible. 
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APPENDIX  A 


HELMHOLTZ  MAGNET  EMP  DESIGN  AND  ANALYSIS  MODEL 

HMM2EMP2.TK 


VLCOP 

251.12174 

mV 

VT 

374.22061 

mV 

EEH 

8.9928881 

% 

POH 

174.00631 

W 

VOLND 

388.9114 

cm'' 3 

VOLNDI 

128.13048 

cm ''3 

VOLNDO 

260.78092 

cm"' 3 

S 

R2 

SELECTW 

306.20446 

1 

cm^2 

mm 

LNDI 

67.991482 

mm 

LNDO 

138.38144 

mm 

LND 

206.37293 

mm 

1200 

TL 

K 

TWO 

727 

C 

MND 

.26252216 

kg 

ML 

.50946919 

kg 

MT 

.10897914 

kg 

MB 

3.7705728 

kg 

MASS 

35.487131 

kg 

1 

TI 

mm 

MTEG 

226.67485 

kg 

EMASS 

262.16198 

kg 

44.5 

SPMTEG 

SRATIO 

VAR 

192.7467 

g/w 

136.01373 

DPTVAR 

kg 

SPMEMP 

6.9667072 

g/w 

295 

TRM 

K 

TWK 

1000 

K 

DENRTW 

8.5734271 

Mg/m^3 

DENW 

8.4286 

Mg/m'' 3 

RHOW 

44.292631 

uOhm-cm 

RHOL 

40.531 

uOhm-cm 

MU 

.2729 

cp 

DENL 

.4618 

Mg/m'' 3 

DENRTL 

.5272515 

Mg/m'' 3 

RHOB 

31.99 

uOhm-cm 

DENSB 

5.89152 

Mg/m^3 

7.3 

DENM 

Mg/m'' 3 

RHOMo 

20.138 

uOhm-cm 

DENMo 

10.1234 

Mg/m'' 3 

KMo 

1.1157822 

W/cm-K 

DENZrB 

5.89152 

Mg/m^3 

RHOZrB 

31.99 

uOhm-cm 

SCO 

Z 

15.5 

cm 

cm 

Z1 

5.65 

cm 

Z2 

9.85 

cm 

RM 

8.5 

cm 

RE 

B. 7401961 

cm 

Terminal  voltage  of  parallel  coils 
Voltage  of  EMP  +  HM 
Electrical  to  hydraulic  efficiency 
Output  power  (hydraulic) 

Volume  of  nozzle  and  diffuser 
Volume  of  nozzle 
Volume  of  diffuser 

Surface  area  of  nozzle  and  diffuser 

Equivalent  radius  of  throat  area 

1=(D-W)  &  0=(D-H)  EQNS  used  for  LNDs 

Length  of  nozzle 

Length  of  diffuser 

Length  of  nozzle  and  of  diffuser 

Temperature  of  fluid 

Temperature  of  the  walls 

Mass  of  nozzles  and  diffusers 

Mass  of  fluid 

Mass  of  pump  throat 

Mass  of  busbar 

Total  mass  of  EMP 

Thickness  of  therm.  &  elect  insulation 
Estimated  Mass  of  TEG 
Estimate  total  TEMP  mass 
Specific  mass  of  the  TEG 

EMASS 

Specific  mass  of  EMP  as  calc,  here 

MATERIALS  PROPERTIES 

Room  temperature 

Temperaure  of  wall  material 

RT  density  of  wall  material 

Density  of  wall  material 

Resistivity  of  throat  wall  material 

Resistivity  of  pumped  fluid 

Viscosity  of  fluid 

Density  of  fluid 

RT  density  of  fluid 

Resistivity  of  busbar 

Density  of  busbar 

Density  of  PM  material 

Resistivity  of  Mo 

Density  of  Mo 

Thermal  conductivity  of  Mo 

Density  of  zirconium  diboride 
Resistivity  of  zirconium  diboride 

HMl  HELMHOLTZ  MAGNET  DESIGN 
06/12/91  1730 

Separation  magnet  coils  center  lines 
Axial  displacement  to  pt  of  interest 
Z  to  center  of  nearest  pump  cell 
Z  to  center  of  2nd  pump  cell 
Mean  coil  radius 
Equivalent  radius  of  coils 
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C 

- 

Cl 

1.5469374 

- 

C2 

.88732955 

- 

Hm 

Oe 

HOCl 

267.80434 

Oe 

H0C2 

132.19566 

Oe 

HOC 

400 

Oe 

LG 

6289. 

,5481  ICO 

A 

L 

40 

B 

mT 

BT 

mT 

L 

7 

ETA 

cm 

SR 

cm 

1 

MMU 

- 

RHOCO 

31.99 

uOhm 

L 

DENCO 

5.89152 

Mg/m 

RCO 

3.9927E-5 

Ohm 

LCO 

61.157075 

cm 

ACO 

49 

cm''2 

CD 

128.35812 

A/cm 

RCOP 

3.9927E-5 

Ohm 

RCOS 

Ohm 

PCOS 

W 

PCOP 

3158.8845 

W 

VLCOS 

Volt 

2 

NCELL 

- 

RI 

5 

cm 

RO 

12 

cm 

VCO 

2616.9467 

cm^3 

L 

MCO 

30.835587 

kg 

Rl 

cm 

R 

0 

cm 

R12 

k 

- 

X 

deg 

Y1 

- 

Y2 

- 

0 

fR 

- 

RFLl 

- 

RFL2 

RFLUX 

RFLUX2 

RFLUX3 

RFLUX4 

RFLUX5 

MMB 

Separation  factor,  coils  to  duct  CL 
Separation  factor,  coil  1  to  duct  CL 
Separation  factor,  coil  2  to  duct  CL 
Magnetic  field  intensity 
Magnetic  field  intensity  from  coil  1 
Magnetic  field  intensity  from  coil  2 
Magnetic  field  intensity  both  coils 
Coil  current 
Magnetic  flux  density 

Radial  width  of  turn  cross-section 
Slant  radius 

Ratio  of  B/Hm  in  the  coil  flux  field 

Resistivity  of  coil  material 

Density  of  the  coil  material 

Resistance  of  one  coil 

Length  of  each  coil  and  lead 

Cross-sectional  area  of  the  coils 

Current  density  in  coils 

Resistance  of  parallel  connected  coil 

Resistance  of  series  connected  coils 

Power  input  to  series  coils 

Power  input  to  parallel  coils 

Terminal  voltage  of  series  coils 

Single  or  Double  Pass  Pump 

Inside  coil  radius 

Outside  coil  radius 

Volume  of  coil 

Mass  of  the  coils 

Slant  radius  to  variable  location 

Variable  radius  position 

Slant  radius  opposite  from  R1 

Modulus  of  the  elliptical  integrals 

Angle  associated  with  k 

Ellitical  integral  of  1st  kind 

Ellitical  integral  of  2nd  kind 

Frational  value  of  RE 

1ST  TERM  OF  RFLUX 

2ND  TERM  OF  RFLUX 

Radial  magnetic  flux  normalized  (0,C) 
Radial  magnetic  flux  normalized  (0,C) 
Radial  magnetic  flux  normalized  (0,C) 
Radial  magnetic  flux  normalized  (0,C) 
Radial  magnetic  flux  normalized  (0,C) 
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RULE  SHEET  MMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMM 
I  S  RuleDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 
■  "HELMHOLTZ  MAGNET  EMP  DESIGN  07-01-91  1600 


* 


* 


DPDV=M1*B*IT 

DPMB=M2>*<B^2*Q 

DPMF=M3*B^2*Q 

DPHY=(DPT+DPNAD) 

TDPDL=DPDV-DPMB-DPMF-DPHY 

CEMF=B*W*Q/AF 

IC=CEMF/(RB+RL) 

IL=IT*RB/(RB+RL) 

AF=W*H 

RL=RHOL5klO^-8W(H*LT) 

RB=RF*RW/(RF+RW) 

RF=FX5|«2 . 6  ^KRHOL*  1 0 -  8 /H 
RW=RHOW>K  1 0  "  -  8>K  ( W+  2*TW )  /  ( TW*  2*LT ) 

QE=15850*Q 

M1=RB/(RB+RL)*1/H 

M2=1/(RB+RL)*1/H"2 

M3=1/(H^2*RF) 

M4=l 

NRE=DENL*V*4*RH/(MO*10^-6) 

RH=W*H/(2*(W+H)) 

V=Q/AF 

DPCV=DENL*VMAX^2/2 

f=.043/NRE".2  "For  5000<NRE<200000 

HYHD=DENL*V^2/2 

DPT=f  *LT/(  45XRH )  *HYHD 

DPNAD=(CN+(1-E)*(1-AR''2)  )*HYHD 

AR=4WH/(PI()*D''2) 

RTW=RHOW*10 ''-8*2*TW/ ( H*LT  ) 

PIE=VL*IT 
PIT=PIE+PCOP 
BMIN=IT/Q*Ml/( 2* (M2+M3 ) ) 

"B=BMIN 

EEH=POH/PIE*100 
MT=NCELL*2*  ( H+W )  *TW>«LT*DENRTW 
ML=NCELL*  ( HWLT+VOLND )  J^tDENRTL 
" ML= 2* ( H*W*LT+4* VOLND ) *DENRTL 

S=PI()>t'(H^2+W^2)^(l/2)*(LNDI/COS(PI()/12)+LNDO/COS(PI()/24)) 

"S=(PI()*D+(W+H))*(LNDI+LNDO)/2 

"  S=PI  ()><<{  D/2+R2 )  *  ( LND^2+  (D/2-R2 )  ^  2 )  ^  ( 1/2 ) 

SELECTW=STEP(H,W) 

( 1  -SELECTW )  *LNDI  =  ( D-H )  /  ( 2>t«TAN  (PI()/12))*(  1  -SELECTW ) 

( 1 -SELECTW )  *LNDO=  ( D-H )  /  ( 2*TAN  ( PI  ( )  /24 ) )  >•<  ( 1 -SELECTW ) 
SELECTW*LNDI=(D-W)/(25KTAN(PI(  )/12)  )*SELECTW 
SELECTW*LNDO= ( D-W ) / ( 2*TAN ( PI () /24 ) ) ^SELECTW 
LND=LNDI+LNDO 
"LND=(D/2-R2)/TAN(PI( )/12) 

MND=NCELL*S*TW*DENRTW 
"MND=4*S^TW*DENRTW 
VOLNDI=(PI( )*D^2/4+(W*H) )*LNDI/2 
VOLNDO=(PI()*D"2/4+(W>KH)  )3XLNDO/2 
VOLND=VOLNDI+VOLNDO 

"VOLND=PI(  )*LND/3*(  (D/2)"2+D/2’i<R2+R2"2) 

"R2=(W^2+H^2)^(l/2)/2 
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"MATERIALS  PROPERTIES 


*  TWC=TL-273  ■ 

*  TWK=TL 

*  DENW=8. 6304-1. 895E-4*TWK-1.23E-8*TWK^2  I 

*  DENRTW=8. 6304-1. 895E-4*TRM-1.23E-8*TRM''2  I 

*  RHOW=16. 337+4. 224E-2*TWC-4.922E-6*TWC^2-3.941E-10*TWC^3  "For  Nb-Zr  in  uOhm-cm 


*  RHOL=13. 735+2. 9256*10"-2>KTL-2.46*10^-6*TL"2 

*  DENL=.  5593-1 . 133E-4*TL+1 . 58E-8=*:TL^2 

*  DENRTL= .  5593-1 . 133E-4*TRM+1 . 58E-8>i<TRM^2 

*  MU=1 . 8739-4 . 667E-3*TL+4 . 933E-6*TL^2-1 . 867E-9*TL"3 

*  RH0Mo=(-0. 183+2. 0321E-2*TL) 

"RHOMo=-3 . 174+ . 03171*TL 

*  DENMo=10 . 2698-1 . 056E-4>KTL-4 . 08E-8>»:TL^2 

*  KMo=(0. 3602-1. 142E-4>KTL+2.050E-8*TL^2)*4. 1868 


"For  Li  in  micro-Ohm-cm 
"For  Li  in  Mg/m''3 
"For  Li  in  Mg/m''3 
"For  Li  in  centipoise 
"For  Mo  in  uOhm-m 
"For  Mo  in  uOhm-cm 
"For  Mo  in  Mg/m '‘3 
"For  Mo  in  W/cm-K 


"RHOL=36 . 768-1 . 6212E-2*TL+7 . 1741E-5*TL-'2 
"DENL= .  9390-2 . 426E-4>KTL 
"DENRMLr . 9390-2 . 426E-4*TRML 

"M0=1. 1914-2. 7431E-3*TL+2.5463E-6*TL^2-8.492E-10*TL'' 


"For  NaK  in  micro-Ohra-cm 
"For  NaK  in  Mg/m^3 
"For  NaK  in  Mg/m 3 
3  NaK  in  Mg/m" 3 


*■  DENZrB=6.08*.969 
*  RH0ZrB=-3.71+.0357*TL 


"For  ZrB2  in  Mg/cm"3 
"For  ZrB2  in  uOhm-cm 


"SINGLE-PASS  PUMP  SECTION 


"DPDL=DPDLE  "SP 
"VLL=IT^RB/(RB+RL)>KRL  "SP 
"VLW=IT*RB/(RB+RL)*RTW  "SP 
"VL=VLL+VLW+CEMF+VBS  "SP 
"POH=DPDL5*cQ  "SP 
"LB=H/2+TI+3/2>t<TB+W  "SP 
"MB=<  LB-TB/2+H/2 )  *TB>t=LT*DENSB  "SP 


"DOUBLE-PASS  PUMP  SECTION 


*  DPDL=DPDLE/NCELL 

*  VLL=NCELL*IT*RB/(RB+RL)*RL 
VLW=NCELL*IT*RB/(RB+RL)5KRTW 

*  VL=VLL+VLW+VC+VBS 

*  VT=VL+VLCOP 

*  VC=NCELL*CEMF 

*  POH=NCELL*DPDL*Q 
"LB=TB+TI+H 
"AB=LT*TB 

*  LB=LT/2  "DP 

*  AB=2WTB  "DP 

^  RHOB=RHOZrB 

*  RBS=RHOB*10"-8>i<LB/AB 

*  VBS=IT*RBS 

*  DENSB=DENZrB 

*  MB=2*LB>t^AB*DENSB 

*  MASS=MCO+MB+MT+ML+MND 

*  EMASS=MTEG+MASS 

*  MTEG=SPMTEG*PIT 

*  SRATIO=100>^EMASS/DPTVAR 


"DP 

"DP 

"DP 

"DP 

"DP 

"DP 

"DP  Backside  busbar 
"DP  Backside  busbar 
Backside  connections  to  series 
Backside  connections  to  series 


"DP 


only 

only 

el ectromagnet 
electromagnet 
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.  *  SPMEMP=MASS/PIT 
■  "HMl  -  HELMHOLTZ  MAGNET  DESIGN 


RM=(LT+ETA)/2 

R1^2=(RE+R)''2+Z^2 

R12^2=(RE-R)"2+Z"2 

k^2=l-(R12/Rl)^2 

X=180/PI()>KASIN(k) 

Y1=RELATE1(X) 

Y2=RELATE2(X) 

Z 1 =H/2+TI +TW+ETA/2 
Z2=l .  5>t<H+2*TI+3*TW+ETA/2 
RE=RM* ( 1+ETA^  2/ ( 2  4*RM^  2 ) ) 

C1=RE/Z1 

C2=RE/Z2 

SCO=Zl+Z2 

•■Hin=2*PI(  )*1000*ICO*RE^2/(10*SR^3) 

H0C1=2*PI(  )*ICO/10>t<C1^3/(RE*(C1^2+l)^1.5)=t:10"-2 
H0C2=2*PI(  )’t=ICO/10*C2^3/(RE*(C2"2+l)^1.5)>K10"-2 
•’SR=(RE^2+Z^2)''0.5 
"B=2*MMU*Hin 

H0C=H0C1+H0C2  «  „  „ 

B=MMU*H0C/10"4  "10e/45*'pi ( )5t:10"-3  Oe/(A/in)*4*pi(  )’*'10"-7  T-in/A*10"3  mT/T  =.l 

'’BT=MMU*H0C/10"4  "10e/45t:pi(  )*10"-3  Oe/(A/m)*4*pi ( )*10''-7  T-m/A*10''3  mT/T 

RHOCO=RHOZrB 

RCO=RHOCO* 10"- 8*LCO/ACO 

RI=LT/2 

RO=RI+ETA 

VCO=PI( )*(RO"2-RI"2)*ETA 

DENCO=DENZrB 

MCO=NCELL*VCO>t:DENCO 

ACO=ETA"2 

CD=ICO/ACO 

"RC0S=2*RC0 

RCOP=RCO 

" VLCOS =1000*1 CO*RCOS 

VLCOP=ICO*RCOP 

PCOS=ICO"2*RCOS 

PCOP=2*ICO"2*RCOP 

"MCO=ACO*LCO*2*DENCO 

LCO=2*PI( )*RM+SCO/2 

■■Z=RE/C 

R=fR*RE 

•'  gr  2^Z 

RFL1=(C"2+1)"1.5/(PI( )*C"2*(C"2*(l+fR)"2+l)".5) 

RFL2=2*C"2*(l-fR)*Y2/(C"2*(l-fR)"2+l)+Yl-Y2 

'‘RFLUX=RFL1*RFL2 

•■RFLUX2=RFL1*RFL2 

"RFLUX3=RFL1*RFL2 

••RFLUX4=RFL1*RFL2 

RFLUX5=RFL1*RFL2 
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HELMHOLTZ  MAGNET  TEG 
HMM2TEG2 . TK 


HMMI-lMI'lI’lIiMI'MIMI'lMMMMMI'l  VARIABLE  SHEET 

St  InputDDDD  NameDDD  OutputDDD  UnitDDDDD  CommentDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 

HMM2TEG2  -  THERMOELECTRIC  GEN.  DESIGN 
SI-GE  +  GAP  MATERIALS  26kB 

5/13/88  2140 

Tcc-72,L-81,  DTl-175 
HEAT  EXCHANGER  DESIGN 


HR 

m 

DELTATl 

V 

NRE 

DH 

942.25393  AF 
PER 

28.553149  H 
W 
LT 

HYHD 
FF 
Q 
QE 
QHX 
HX 

AFHX 
PEN 
NUN 
PRN 

WDELTAT 
TWI 
TWO 
HFLX 
FDELTAT 
tDW 
CONDW 


15.945295  kW 
2.8009404  kg/s 
1.4158034  K 
6.6958047  m/s 
361831 

30 . 615945  mm 

mm^2 

123.1063  mm 
mm 

33  mm 

mm 

9951.9616  Pa 
.00355678  - 
6.3091483  Litre/s 
gpm 

15.945295  kW 
6.2589593  W/cm^2-K 
131.34449  cm^2 
5948.7223  - 
31.150438  - 
.01644061  - 
8.6714692  K 
1208.6715  K 
1200  K 

121.40057  W/cm''2 
19.396287  K 
mm 

W/cm-K 


230 

100 

.5 

.7 


Heat  transport  rate  by  fluid 
Fluid  mass  flowrate 
Temperature  drop  inlet  to  oulet 
Velocity  of  fluid  flow 
Reynolds  number 
Hydraulic  diameter 

Cross-sectional  area  in  flow  stream 

Perimeter  of  flow  duct 

Height  of  flow  duct 

Width  of  flow  duct 

Duct  length 

Hydraulic  head 

Friction  factor 

Volumetric  flowrate 

Volumetric  flowrate 

Heat  transfer  rate  fluid  to  wall 

Heat  transfer  coefficient 

Heat  transfer  area 

Peclet  number 

Nusselt  number 

Prandlt  number 

Temperature  drop  across  duct  wall 
Temperature  of  inside  of  wall 
Temperature  of  outside  of  wall 
Heat  flux  through  wall 
Temperature  drop  across  fluid  film 
Thickness  of  duct  wall 
Thermal  conductivity  of  wall  material 


DTLQ 

29.483559 

K 

1228.6  TFL 

K 

COND 

.61515654 

W/cm-K 

MU 

.25152306 

cp 

DENL 

.44394906 

Mg/m^3 

Cp 

.96056449 

cal/g-K 

RHOMo 

11.581969 

uOhm-cm 

DENMo 

10.194986 

Mg/m^3 

CONDMo 

1.2602791 

W/cm-K 

5.67E-12 

SIG 

W/cm^2- 

256 

Te 

K 

.00053045 

S 

V/K 

.05799219 

K 

W/cm-K 

.00529696  RHO  Ohm-cm 


MATERIALS  PROPERTIES 

Temperature  of  fluid 

Temperature  of  fluid 

Thermal  conductivity  of  fluid 

Viscosity  of  fluid 

Density  of  fluid 

Heat  capacity  of  fluid 

Resistivity  of  Mo 

Density  of  Mo 

Thermal  conductivity  of  Mo 

THERMOELECTRIC  CONVERTER  DESIGN 

Boltzmann  constant 

Effective  temperature  near  Earth 

Couple  Seebeck  Coef . 

Couple  thermal  conductivity 
Couple  resistivity 
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O  O  O  •  O  O  M 


0 


RHOcc 

Rhc 

Rcc 

Qte 

Qbp 

Qs 

Qrj 

Qst 

Qr 


0  Ts 

0  The 

1200  Th 

Tm 

645  Tc 

L  Tee 

Tei 

0  Tep 

Thp 
ET 
EA 
Mo 
IG 

10  I 

.15017817  L 
A 

1  NS 

L  NP 

NT 
Vc 
Vee 

L  VL 

P 

L  PO 

Aet 
AR 
Khi 
Khe 

0035  Kli 

Kce 
Kei 
Kep 
Khpj 

.7  Khp 

Ae 
Ahe 
Aj 

0  Lhi 

.6  b 

Lee 

.02  Lli 

0  Lei 

0  Lep 

.  1  tl 

0  Lj 

0  Lhpw 

.3  F 

.  1  Fbp 


ME 


1.1582E-5  Ohm-em 
Ohm 

6.1146E-5  Ohm 
19.398138  W 
1.3656913  W 
20.763829  W 
19.080411  W 
15945.295  W 
14652.538  W 
K 
K 
K 

922.5  K 
K 

578.95621  K 
K 
K 

557.53268  K 
.08741286 
.08136901 
1.3583095 
23.567306  A/em 
A 
em 

.06372309  em''2 

767.93618 
767.93618 
.16956467  V 
.00061146  V 
.16956467  V 
1.6895322  W 
1297.4529  W 
131.34449  em^2 


W/cm-K 
W/em-K 
W/cm-K 
1.2602791  W/cm-K 
W/cm-K 
W/cm-K 
W/cm-K 
W/cm-K 

.17103568  cm^2 
.13754356  cm^2 
cm"  2 


cm 

.72855315  cm 


.05585247  kg 


Resistivity  of  cold 
Resistance  of  hot  connector 
Resistance  of  cold  connector 
Thermoelectric  source  heat  rate 
Bypass  heat  in  leg  insulation 
Couple  Source  Heat  Rate 
Couple  reject  heat  rate 
Total  input  heat 
Heat  to  be  rejected  by  radiator 
Heat  Source  Temperature 
Temperature  of  hot  connector 
Hot  Junction  Temperature 
Mean  temperature  across  couples 
Cold  junction  temperature 
Temperature  of  cold  connector 
Temperature  *’  "  insulator 

Temperature  "  compression  pad 
Temperature  of  heat  pipe 
Theoretical  conversion  efficiency 
Actual 

Max.  Efficiency  load  matching  factor 
Current  gradient 

Electric  current  in  Thermoelements 
TE  element  length 
TE  "  cross-sectional  area 
No.  of  series-connected  couples 
No.  of  parallel  circuits 
Total  number  of  couples 
Couple  load  voltage 
Cold  connector  deltaV 
Generator  or  module  terminal  voltage 
Couple  output  power 
Module  output  power 
Total  area  covered  by  thermocouples 
Area  ratio  heat  source/thermocouples 
Conductivity  of  hot  insulator 
Conductivity  of  hot  connector 
Conductivity  of  leg  insulator 
Conductivity  "  cold  connector 
Conductivity  of  cold  insulator 
Conductivity  "  compression  pad 
Conductivity  of  heat  pipe  joint 
Conductivity  "  "  "  wall 

Area  of  couple  for  heat  conduction 
Hot  connector  heat  conduction  area 
"  "  joint  area 

Thickness  of  hot  insulator 
Thickness  of  hot  and  connector 
Length  of  cold  connector  or  busbar 
Thickness  of  leg  insulator 
Thickness  ”  cold  insulator 
Length  of  compression  pad 
Thickness  of  insulation  between  ducts 
Thickness  of  joint 
Thickness  of  heat  pipe  wall 
Factor  of  1D/2D  correction 
Bypass  heat  allowance  factor 
Mass  of  thermoelements 
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MC 

.50522957  kg 

MHLC 

g 

ML  I 

.01256764  kg 

MR 

57.929759  kg 

MCB 

.39258017  kg 

TM 

58.895989  kg 

SPPO 

45.393548  g/W 

3.8 

DE 

g/cm'S 

DC 

10.194986  g/cm''3 

4 

DI 

g/cm^3 

2.5 

DLI 

g/cm"3 

2.7 

DCS 

Mg/m^  3 

dCB 

1.2  cm 

Acb 

131.34449  cm^2 

RG 

1.6256E-5  Ohm 

IM 

7679.3618  A 

VOC 

.29439975  Volt 

WM 

1332.9178  W 

VLP 

.57596744 

WP 

2.4758E-5 

IP 

.00130219 

EAP 

dL 

2.6 

cm 

fl 

.2 

- 

f2 

.4 

- 

f3 

.6 

- 

f4 

.8 

- 

f5 

1 

- 

LI 

2.6 

cm 

L2 

5.2 

cm 

L3 

7.8 

cm 

L4 

10.4 

cm 

L5 

13 

cm 

Q1 

3228.5254 

W 

Q2 

2377.0369 

W 

Q3 

1631.9696 

W 

Q4 

996.75833 

W 

Q5 

461.21832 

W 

Qrl 

851.4885 

W 

Qr2 

745.06729 

W 

Qr3 

635.21124 

W 

Qr4 

535.54001 

W 

Qr5 

461.21832 

W 

9 

el 

- 

Arl 

2152.4425 

cm^2 

Ar2 

2152.4425 

cm^2 

Ar3 

2152.4425 

cm^2 

Ar4 

2152.4425 

cm"' 2 

Ar5 

2152.4425 

cm^2 

Tf 

557.35534 

K 

T1 

534.82298 

K 

T2 

518.23328 

K 

T3 

499.25034 

K 

T4 

479.92665 

K 

'•  ■■  cold  connector  -  radiator  blk 

Not  applicable 

■'  "  leg  insulation 

Mass  of  radiator 
Mass  of  clamping  blocks 
Total  Mass  of  TEG 
Specific  Power 
Density  of  thermoelements 
"  "  connectors 

"  ■'  hot/cold  side  insulators 

"  "  leg  insulation 

”  "  clamping  block 

Thickness  of  clamping  block 
Area  of  clamping  block 
Resistance  of  the  module 
Total  module  current 
Open  circuit  voltage 
Power  dissippated  in  module 
Fraction  of  VOC  dissspated  in  module 
Fraction  of  power  dissspated  in  modul 


WASTE-HEAT  RADIATOR  DESIGN 
RECTANGULAR  RADIATOR  FIN  DESIGN 

Dif f eriential  slab  length 
fractional  length  of  slab  1 

fractional  length  of  slab  2 

fractional  length  of  slab  3 

fractional  length  of  slab  4 

fractional  length  of  slab  5  ■ 

Length  of  slab  1 

Length  of  slab  2 

Length  of  slab  3 

Length  of  slab  4 

Length  of  slab  5 

Heat  entering  slab  1 
Heat  entering  slab  2 
Heat  entering  slab  3 
Heat  entering  slab  4 
Heat  entering  slab  5 
Heat  radiated  from  slab  1 
Heat  radiated  from  slab  2 
Heat  radiated  from  slab  3 
Heat  radiated  from  slab  4 
Heat  radiated  from  slab  5 
emissivity  of  radiator  surface 
Radiation  area  of  slab  1 

Radiation  area  of  slab  2 

Radiation  area  of  slab  3 

Radiation  area  of  slab  4 

Radiation  area  of  slab  5 

Root  temperature  of  fin 
Mean  temperature  of  slab  1 
Mean  temperature  of  slab  2 
Mean  temperature  of  slab  3 
Mean  temperature  of  slab  4 
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T5 

463.83209  K 

DTI 

45.064734  K 

DT2 

33.1794  K 

DT3 

37.965875  K 

DT4 

38.647373  K 

DT5 

32.189125  K 

1.5 

Kf 

W/cm-K 

Wf 

827.86252  cm 

13 

Lf 

cm 

.3 

tf 

mm 

2.7 

DR 

g/cm‘'3 

13 

N1 

- 

5 

N2 

3 

N3 

1.8 

N4 

- 

1 

N5 

- 

FEFF 

6.3755991  - 

CT 

3228.5254  W 

RT 

5063.8776  W 

MF 

18.83225  kg 

Nf 

2.2692308  - 

SHP 

2622.3354 

cm'' 2 

1 

DHP 

cm 

LHP 

834.71527 

cm 

HPLDR 

834.71527 

- 

MP 

1.1151601 

kg 

MW 

15.602895 

kg 

MFL 

.5113554 

kg 

MHP 

17.229411 

kg 

MHPT 

39.09751 

kg 

.5 

tP 

mm 

8.5 

DENP 

g/cm''3 

LE 

6.8527559 

cm 

1 

tw 

nun 

8.5 

DENW 

g/cm''3 

AE 

48.853288 

cm'' 2 

AC 

5901.8308 

cm'' 2 

.3 

fV 

- 

.5 

DENHPL 

g/cm''3 

.3 

fPOR 

- 

HPEVFX 

299.92942 

W/cm^2 

Mean  temperature  of  slab  5 

Temperature  drop  across  slab  1 

Temperature  drop  across  slab  2 

Temperature  drop  across  slab  3 

Temperature  drop  across  slab  4 

Temperature  drop  across  slab  5 

Thermal  conductivity  of  fin  material 

Width  of  fin 

Lenght  of  fin 

Unit  thickness  of  fin 

Density  of  fin  material 

No.  of  unit  layers  in  slab  1 

No.  of  unit  layers  in  slab  2 

No.  of  unit  layers  in  slab  3 

No.  of  unit  layers  in  slab  4 

No.  of  unit  layers  in  slab  5 

Fin  effectiveness 

Radiated  heat  from  slabs 

Radiated  heat  at  fin  root  temperature 

Mass  of  fins 

No.  of  fins  in  radiator 

HEAT  PIPE  MASS 

Surface  area  of  heat  pipe  surface 
Inner  diameter  of  heat  pipes 
Length  of  pipe 

Heat  pipe  length/diameter  ratio 

Mass  of  pipe  wall 

Mass  of  pipe  wick 

Mass  of  fluid 

Mass  of  heat  pipes 

Total  mass  of  heat  pipes 

Thickness  of  pipe  wall 

Density  of  pipe  wall 

Length  of  evaporator  section 

Thickness  of  wick 

Density  of  wick  material 

Heat  pipe  evaporator  area 

Heat  pipe  condenser  area 

Excess  volume  factor 

Density  of  heat  pipe  fluid 

Porosity  factor  of  wick 

Evaporator  wall  heat  flux 
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RULE  SHEET 

S  RuleDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD: 
"THERMOELECTROMAGNETIC  PUMP  TEG  DESIGN  2/22/88  1415 

"Revised  5/17/88 

"HEAT  EXCHANGER  DESIGN 

*  HR=m*10005t«Cp5t<4. 186*DELTAT1 

*  m=DENL*AF*V*1000 

*  NRE=DENL*V*DH/MU=l'1000000 

*  DH=4*AF/PER 

*  AF=W*H 

>K  PER=2*(W+H) 

*  HYHD=DENL^1000*V^2/2 

’I'  FF=.046/NRE^.2  "For  5000<NRE<200000 

3K  Q=AF*V*1000 

*  QE=15.85>t'Q 

*  QHX=HX>t<AFHX*FDELTAT 
HX=NUN*COND/DH*100 

*  PEN=PRN>KNRE 

*  NUN=5+.025>KPEN".8 
"NUN= .  023*NRE^0 . 8>KPRN^0 . 4 

*  PRN=Cp*4. 186*MU/(100*COND) 

>K  HR=QHX 

*  HFLX=QHX/AFHX 

*  AFHX=2*LT>^H 

*  WDELTAT=HFLX*tDW/(CONDW>H100) 

*  WDELTAT=TWI-TWO 

*  HR=Qst 

*  DTLQ=WDELTAT+FDELTAT+DELTAT1 
"MATERIALS  PROPERTIES 

*  MU=.00142>K100*EXP(702.4/TFL) 

*  CONDr. 2725+2. 789E-4*TFL 

*  DENL= . 5593-1 . 133E-4*TFL+1 . 58E-8*TFL^2 

*  Cp=l .  1816-4 . 429E-4*TFL+3 . 035E-7>KTFL^2-7 . 28E-11>KTFL^3 

*  RH0Mo=-0. 183+2. 0321E-2>>:Tcc 

*  DENMo=10 . 2698-1 . 056E-4>KTcc-4 . 08E-8*Tcc^2 

*  CONDMo=(0. 3602-1. 141E-4*Tcc+2.05E-8*Tcc"2)*4. 1868 

"TEG  MODULE  DESIGN 

*  Th=TWO 

"Qs=Khi*Ac* ( Ts-Thc ) /Lhi 
"  Qs=Khc*Ahc>K  ( Thc-Th )  /b 

*  Ac=2*(A+4*Lli>t^(A^(l/2)+Lli)) 

Ahc=2*A+A^(l/2)*2*Lli 

^  ET=  ( Th-Tc )  /Th>^  ( Mo- 1 )  /  ( Mo+Tc/Th ) 

*  Z=S^2/(RHO*K) 

Mo=(l+Z5t<Tin)^(l/2) 

^  Tm=(Th+Tc)/2 

Qte= S*Th>K  I  +K’t^  ( Th-  Tc )  /  ( L/A )  - 1  ^  2  *RHO*  ( L/A )  /2  - 1  ^  2>i=Rhc 
*■  Qs=Qte+Qbp 

*  EA=P/Qs 

*  Rhc=F*RHOhc*  ( 2>t=A^  ( 1  /2 )  +L1  i+b )  /  ( A^  (1  /2 )  *b ) 

*  Qbp=Fbp*K*A5t<(  Th-Tc  )/L 
IG=I*L/A 

*  IG=S*(Th-Tc)/(RHO*(Mo+l)) 

>i<  Vc=S*(Th-Tc)*(Mo/(l+Mo)) 


"For  Li  in  cp 
"For  Li  in  W/cm-K 
"For  Li  in  Mg/m'' 3 
"For  Li  in  cal/g-K 
"For  Mo  in  Mg/m^3 
"For  Mo  in  Mg/m" 3 
"For  Mo  in  W/cm-K 
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■  * 

P=  ( S*  ( Th-Tc )  /  ( 1+Mo ) )  ^2mo/  ( RHO*L/A )  - 1  2*Rcc 

1 

Qr j=Qs*(l-EA)+I"2*Rcc 

■  * 

RHOc  c =RH0Mo2<^  1 0  "  -  6 

:4c 

Rcc=RHOcc’*‘Lcc/  ( LT*b ) 

■ 

Lcc=b+tI+H 

1 

■■Rcc=F*RHOcc*  ( 2*A^  ( 1  /2 )  +Lli+b )  /  ( A^  ( 1/2 )  *b ) 

5K 

Vcc=I*Rcc 

■ 

■'Rcc=Rhc/RHOhc*RHOcc 

B  * 

Kcc=CONDMo 

*  3K 

Qr j =Kcc*Ahc* ( Tc-Tcc ) /b 

* 

Qr=Qrj*NT 

1 

" Qr  j  =Kcp*Ahc* ( Tci -Tcp ) /Lcp 

B 

Qr=Khp*AE* ( Tcc-Thp ) /tP 

Qr=Khp*AC* ( Thp-Tf ) /tP 

1 

" Qr =Khp*Ahp* ( Thp-Tf ) /Lhpw 

1  * 

VL=NS*Vc 

* 

NT=PO/P 

■  * 

NP=NT/NS 

1  * 

Qst=Qs*NT 

Act=Ac*NT 

* 

AR=Act/AFHX/10000 

1  * 

RG=NS*RHO*L/A/NP 

■  * 

IM=NS*S* ( Th-Tc ) / ( RG* ( 1+Mo ) ) 

5»C 

VOC=NS*S*( Th-Tc) 

■  * 

WM=(VOC/2)"2/RG 

1  * 

VLP=VL/VOC 

* 

WP=W/WM 

■  * 

IP=I/IM 

1 

"MASSES 

®  :4c 

ME=  ( DE^K  ( 2*A*L )  /lO  00 )  *NT 

"TE  Elements 

"MC=NT*DC*2*  ( 2*A"  ( 1/2 )  +L1  i )  *  ( A"'  ( 1/2 )  *b ) 

"Connectors 

1  * 

DC=DENMo 

1 

"MC=NT*DC5»<  ( 2*A"  ( 1/2 )  +Lli+tI )  >•'{  A"  ( 1/2 )  *b)/1000  "Pump 

connector  -  radiator 

blk 

MC=  ( 2*H+3>*=tI+4*tDW)*b*LT’>:DC 

■  * 

MLI=(2*DLI*L*(4*Lli>*‘{A^(l/2)+Lli))  )*NT/1000 

"Leg  Insulation 

MCB=DCB/10 3*  ( Acb*dCB*10 ''4-Nf *PI  ( )  *DHP'' 2/4*LE ) 

"Clamping 

Blocks 

:4c 

dCB=1.2^DHP 

»  * 

Acb=2*H*LT5t:AR 

1 

" TM= NT* ( ME+MC+MHC I +ML I +MCB ) +MR+MRHP 

"Total  Mass  of  TEG 

B  * 

TM=ME+MC+MLI+MR+MCB 

Total  Mass  ox 

* 

SPPO=TM/PO 

1 

"WASTE  HEAT  RADIATOR  DESIGN 

"RECTANGULAR  RADIATOR  FIN  DESIGN 

1  :4c 

Nf=2*(LT*10''2+Lf/2)/(2*Lf ) 

"  :|c 

HPEVFX=Qr/AE 

:4c 

HPLDR=LHP/DHP 

n  * 

LE=(2*H*10^2*AR)*1.2 

:4c 

dL=Lf/5 

fl=Ll/Lf 

f2=L2/Lf 

:K 

f3=L3/Lf 

f4=L4/Lf 

1* 

f5=l 

*  54c 

Ll=dL 

1 

L2=2*dL 

1 
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*  L3=3*dL 
*■  L4  =  4*dL 

*  L5=Lf 

*  Ql=Qrl+Q2 

*  Q2=Qr2+Q3 
*■  Q3=Qr3+Q4 

*  Q4=Qr4+Q5 

*  Q5=Qr5 

*  Qrl=SIG*el*Arl*(T1^4-Te'4) 

5k  Qr2=SIG*el5kAr2>K(T2''4-Te"4) 

*  Qr3=SIG5kel5kAr3*(T3^4-Te^4) 

*  Qr4=SIG5kel*Ar4*(T4''4-Te"4) 

*  Qr5=SIG>Kel*Ar5>k(T5^4-Te^4) 

>i<  Arl=Wf>KdL 

5k  Ar2=Wf5kdL 
5k  Ar3=Wf5kdL 

*  Ar4=Wf5kdL 
5k  Ar5=Wf5kdL 

5k  Tl=(2*Tf-DTl)/2 
5k  T2=(25kTl-DT2)/2 
5k  T3=(2:kT2-DT3)/2 
5k  T4={25kT3-DT4)/2 
5k  T5=(25kT4-DT5)/2 
5k  DTl=Ql5kdL/(Kf5kWf5kN25ktf ) 

:k  DT2=Q2*dL/(Kf5kWf5kN25ktf ) 

5k  DT3=Q35kdL/(Kf5kWf5kN35ktf ) 

:k  DT4=Q45kdL/(Kf5kWf5kN45ktf ) 

5k  DT5=Q55kdL/(Kf5kWf5kN55ktf ) 

5k  FEFF=105kCT/RT 
5k  CT=:Qrl+Qr2+Qr3+Qr4+Qr5 
5k  RT=SIG5kel5kWf5kLf5k(Tf"4-Te^4) 

5k  MF=Wf 5kdL5ktf 5k (N1+N2+N3+N4+N5 ) 5kDR5k25kNf /lOOO 
5k  Nf=Qr/(25kQl) 

5k  MR=MF+MHPT 

"HEAT  PIPE  MASS 

5k  SHP=PI()5kDHP5kLHP 
:k  MHPT=MHP5kNf/1000 
5k  MHP=MP+MW+MFL 

5k  MP=  ( SHP+PI  ( )  5kDHP^ 2/2 )  5ktP5kDENP 
5k  LHP=Wf+LE 

5k  MW=SHP5ktW5k(l-fPOR)5kDENW 
5k  MFL=SHP5ktW5kDENHPL5kfPOR5k(l+fV) 

5k  AE=PI(  )5kDHP5kLE5kNf 
5k  AC=PI(  )5kDHP5kWf5kNf 
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UNIT  SHEET 

FromDDDDD  ToDDDDDDD  Multiply  ByDD  Add  OffsetDDD 


m 

mm 

1000 

in''3/s 

Litre/s 

1000 

Ta^2 

mm^2 

1000000 

kW 

W 

1000 

Mg/s 

kg/s 

1000 

N 

kg-m/s"2 

cp 

Pa-s 

.001 

cal/g-K 

W-s/g-K 

4.186 

W/cm^2-K 

W/m^2-K 

10000 

W/cin''2 

W/m"2 

10000 

kg 

S 

1000 

cm 

mm 

10 

m'"2 

mm '‘2 

1000000 

m^2 

cm^2 

10000 

Mg 

kg 

1000 

hg 

kg 

10 

V 

dV 

10 

kg/W 

g/w 

1000 
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APPENDIX  C 


LINEAR  MAGNET  EMP 
LMM2EMP.TK 


VARIABLE  SHEET  MMMMMMMMI'MMMMMMMMMMMMMMMMMMMhM 
St  InputDDDD  NameDDD  OutputDDD  UnitDDDDD  CommentDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 

FILE:  LMM2EMP2.TK  -  34.5  kB 
CORELESS  LINEAR  CONDUCTION  PUMP 
DESIGN  AND  PERFORMANCE  PROGRAM 
07-24-91 

Delivered  pressure  from  cell 
Volumetric  flowrate 
Delivered  pressure  from  cell 
Trial  DPDL 
Volumetric  flowrate 

Differential  pressure  rise  developed 
Pressure  loss  from  magnetic  braking 
Pressure  loss  from  magnetic  fringe 
Pressure  loss  from  fluid  friction 
Pressure  loss  in  throat  flow 
Pressure  loss  in  nozzle  and  diffuser 
Magnetic  flux  noraml  to  Q  and  I 
Optimum  magnetic  flux 
Total  pump  current 
Electric  current  in  fluid 
Inside  height  of  throat  duct 
Inside  width  of  throat  duct 
Length  of  throat 
Velocity  of  liquid  in  throat 
Velocity  limited  by  cavitation 
First  term  coefficient 
'2  Second  term  coefficient 
'2  Third  term  coefficient 
Fourth  term  coefficient 
Throat  cross-sectional  flow  area 
Counter  current  from  CEMF 
Electrical  resistance  bypassing  fluid 
Electrical  resistance  across  W 
Resistance  of  fringe  path 
Empirical  factor  from  MSA  calcs 
Thickness  of  throat  walls 
Reynolds  Number  in  the  throat 
Hydraulic  radius  in  the  throat 
Pressure  margin  for  cavitation  limit 
Fluid  friction  factor  in  throat 
Hydraulic  velocity  head  in  throat 
Discharge  coefficient  of  nozzles 
Diffuser  recovery  coefficient 
Area  ratio  throat  to  pipe 
Diameter  of  circuit  piping 
Resistance  of  throat  wall  at  busbar 
Resistance  of  throat  walls 
Input  power  to  EMP  (electrical) 

Counter  EMF  of  fluid  moving  across  B 
Counter  EMF  of  fluid  moving  across  B 
Voltage  drop  across  liquid 
Voltage  drop  across  busbar 
Voltage  drop  across  wall 
Load  voltage  of  TEG  and  EMP 
Electrical  to  hydraulic  efficiency 
Output  power  (hydraulic) 


L 

4 

DPDLE 

psi 

L 

100 

QE 

gpm 

DPDL 

13.79 

kPa 

TDPDL 

kPa 

Q 

6.3091483 

litre 

DPDV 

16.883366 

kPa 

DPMB 

. 26159966 

kPa 

DPMF 

.13640714 

kPa 

DPHY 

2.6953594 

kPa 

DPT 

.43080041 

Pa 

DPNAD 

2.264559 

Pa 

L 

24 

B 

mT 

BMIN 

509.03751 

mT 

LG 

15.324931 

IT 

kA 

IL 

14.144309 

kA 

L 

H 

20.106383 

mm 

L 

47 

W 

mm 

L 

400 

LT 

mm 

V 

6.6763474 

m/s 

VMAX 

6.6958047 

m/s 

Ml 

45.903868 

1/m 

M2 

71985324 

1/Ohm' 

M3 

37535646 

1/Ohm 

M4 

1 

- 

L 

945 

AF 

mm'"  2 

IC 

.21915929 

kA 

RB 

3. 1715E-5 

Ohm 

RL 

2.6473E-6 

Ohm 

RF 

.0000659 

Ohm 

L 

1.125 

FX 

_ 

.5 

TW 

mm 

NRE 

556074.9 

- 

RH 

7.041059 

mm 

10000 

DPCV 

Pa 

f 

.003051 

- 

HYHD 

9941.9664 

Pa 

.05 

CN 

- 

.8 

E 

- 

.33333333 

AR 

— 

D 

60.080231 

mm 

RTW 

6.335E-8 

Ohm 

RW 

.00006114 

Ohm 

L 

PIE 

2261.4391 

W 

VC 

15.06184 

mVolt 

CEMF 

7.5309198 

mVolt 

VLL 

74.887951 

mVolt 

VBS 

55.824138 

mVolt 

VLW 

1.7920963 

mVolt 

L 

VL 

147.56602 

mVolt 

L 

EEH 

7.6944946 

% 

L 

POH 

174.00631 

W 

Cl 


VOLND 

427911.07 

mm'' 3 

VOLNDI 

.00014098 

mm^3 

VOLNDO 

.00028693 

mm  "'3 

S 

36993.776 

mm'"  2 

LNDI 

.07459222 

mm 

LNDO 

.15181576 

mm 

LND 

226.40797 

mm 

DENRML 

Mg/m '"3 

295 

TRML 

K 

LB 

2.5606383 

cm 

5 

TB 

mm 

RBS 

3.6427E-6 

Ohm 

.5 

TI 

mm 

1200 

TL 

K 

TWC 

927 

C 

MND 

.15858172 

kg 

ML 

.42491782 

kg 

MT 

.23013267 

kg 

MB 

.88162774 

kg 

MASS 

21.277136 

kg 

295 

TRM 

K 

TWK 

1200 

K 

DENRTW 

8.5734271 

Mg/m^3 

DENW 

8.385288 

Mg/m^3 

RHOW 

50.949923 

uOhm-cm 

RHOL 

45.2998 

uOhm-cm 

MU 

. 150844 

cp 

DENL 

.446092 

Mg/m^3 

DENRTL 

.5272515 

Mg/m'^  3 

RHOB 

28.451506 

uOhm-cm 

DENSE 

10.084328 

Mg/m^3 

I 

1899.052 

A 

.5 

fW 

WFB 

10.6 

cm 

fLl 

.22222222 

- 

fL2 

.77777778 

- 

LI 

16 

cm 

L2 

56 

cm 

.5 

fX 

- 

72 

X 

cm 

UP 

0 

- 

PUMP 

1 

- 

DOWN 

0 

- 

PHIl 

1.3117506 

rad 

PHI  2 

1.3117506 

PHI  3 

1.3117506 

PHI  4 

1.3117506 

HFl 

.96663486 

HF2 

.96663486 

HF3 

.96663486 

HF4 

.96663486 

Volume  of  nozzle  and  diffuser 
Volume  of  nozzle 
Volume  of  diffuser 

Surface  area  of  nozzle  and  diffuser 

Length  of  nozzle 

Length  of  diffuser 

Length  of  nozzle  and  of  diffuser 

Density  of  liquid  at  room  temperature 

Room  Temperature  of  Liquid 

Length  of  armature  busbars 

Thickness  of  connecting  busbars 

Resistance  of  armature  busbar 

Thickness  of  electrical  insulation 

Temperature  of  fluid 

Temperature  of  the  walls 

Mass  of  nozzles  and  diffusers 

Mass  of  fluid 

Mass  of  pump  throat 

Mass  of  connecting  busbars 

Total  mass  of  EMP 


MATERIALS  PROPERTIES  1 

Room  temperature 
Temperaure  of  wall  material 
RT  density  of  wall  material 
Density  of  wall  material 
Resistivity  of  throat  wall  material 
Resistivity  of  fluid 
Viscosity  of  fluid 
Density  of  fluid 

RT  density  of  fluid  ' 

Resistivity  of  busbar 
Density  of  busbar 

MAGNETIC  FIELD  CALCULATIONS  7-22-91 

Field  Current 
Fraction  of  WFB 

Center-center  span  between  field  buses 

Fractional  location  of  LI 

Fractional  location  of  L2 

Inlet  end  of  Field  busbars 

Discharge  end  of  Field  busbars 

Fraction  of  total  pump  length 

Total  pump  length 

Upstream  region  indicator 

Pumping  section  region  indicator 

Downstream  region  indicator 

Angular  integration  limit 

Angular  integration  limit 

Angular  integration  limit 

Angular  integration  limit 

Integrated  flux  fraction 

Integrated  flux  fraction 

Integrated  flux  fraction 

Integrated  flux  fraction 


C2 


^  ^  ^  ^  ^  ^ 


1.257E-6 


1 

L 

L  4.6 
L 

LG  1899.052 


L 


2 

1 


L 


45 


R1 

5.6956024 

cm 

Radial  distance  -  busbar  1  to  P 

R2 

5.6956024 

cm 

Radial  distance  -  busbar  2  to  P 

HFPRl 

1.9332697 

Total  integrated  fractions,  bar  1 

HFPR2 

1.9332697 

Total  integrated  fractions,  bar  2 

HFPl 

3.3352055 

mT 

Magnetic  Flux  from  bar  1 

HFP2 

3.3352055 

mT 

Magnetic  Flux  from  bar  2 

MUO 

T-m/A 

Absolute  permeability 

HAXVl 

.93054248 

Axial  flux  vector  for  HFPl 

HAXV2 

.93054248 

Axial  flux  vector  for  HFP2 

HFP 

12.895704 

mT 

Total  magnetic  flux 

HFPAX 

12 

mT 

Total  axial  magnetic  flux 

fH 

— 

Fraction  of  height  between  busbars 

HFB 

2.0856383 

cm 

Height  of  the  busbars 

ETA 

cm 

Cross-section  side  of  large  BBs 

ABB 

9.5939362 

cm'' 2 

Cross-sectional  area  of  large  busbars 

CD 

197.94294 

A/cm^2 

Current  Density 

IBB 

A 

Field  Current 

RHOBB 

28.451506 

uOhm-cm 

Resistivity  of  Busbars 

RHOMo 

28.451506 

uOhm-cm 

Resistivity  of  Mo 

VBB 

970.90634 

cm^3 

Voliime  of  the  LM  busbars 

DENBB 

10.084328 

Mg/m'' 3 

Density  of  Busbars 

DENMo 

10.084328 

Mg/m" 3 

Density  of  Mo 

DENZrB2 

5.89152 

Mg/m'' 3 

Density  of  ZrB2 

MLM 

19.581876 

kg 

Mass  of  the  Linear  Magnet  (busbars) 

RBB 

.00030012 

Ohm 

Resistance  of  the  busbars  (magnet) 

LBB 

101.2 

cm 

Length  of  LM  coil  turn 

NLM 

- 

Number  of  Linear  magnets 

NCELL 

- 

Number  of  pump  cells 

VLBBP 

569.93562 

mV 

LM  load  voltage  -  parallel  connected 

RBBP 

.00030012 

Ohm 

LM  resistance  -  parallel  connected 

VLBBS 

mV 

LM  load  voltage  -  series  connected 

RBBS 

Ohm 

LM  resistance  -  series  connected 

PLMP 

2164.6747 

W 

LM  power  -  parallel  connected 

PLMS 

W 

LM  power  -  series  connected 

KMo 

1.0579206 

W/cm-K 

Thermal  conductivity  of  Mo 

SELECTW 

0 

Selects  nozzle  equations  1=H>W  ,  0=W>H 

RHOZrB2 

37.03 

uOhm-cm 

Resistivity  of  ZrB2 

PIT 

4426.1138 

W 

Total  EMP  power  input 

SPMEMP 

4.8071823 

g/w 

Specific  mass  of  the  EMP 

SPMTEG 

g/w 

Specific  mass  of  TEGs 

MEMPTEG 

101.76476 

kg 

Estimated  mass  of  the  TEG  for  armature 

MLMTEG 

97.410363 

kg 

Estimated  mass  of  the  TEG  for  LM 

EMASS 

220.45226 

kg 

Estimated  mass  of  the  LM-TEMP 
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RULE  SHEET 

S  RuleDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDl 
■■CORELESS  LINEAR  CONDUCTION  EMP  -  07-24-91  -  DESIGN  AND  PERFORMANCE  ■ 

DPDV=Ml>kB>KlT  ■ 
DPMB=M2*B^2*Q  ■ 
DPMF=M3>KB''2>KQ 

DPHY=(DPT+DPNAD)>K1000  m 
DPDL=DPDV-DPMB-DPMF-DPHY  ■ 
■■TDPDL=DPDV-DPMB-DPMF-DPHY  * 
CEMF=B*W*Q/AF  ^ 
IC=CEMF/(RB+RL)  ■ 
IL=IT*RB/(RB+RL)  ■ 
AF=W*H 

RL=RHOL=<<10^-8*W/(H*LT)  ■ 
RB=RF*RW/(RF+RW)  ■ 
RF=2 . 6E-8^FX>KRHOL/H 

RW=RHOW*  10^-8*(W+2*TW)/(TW=«2*LT)  - 
QE=15850*Q  I 
M1=RB/(RB+RL)*1/H  ■ 
M2=1/(RB+RL)*1/H"2 

M3=l/(H''2>i<RF)  ■ 
M4=l  1 
NRE=DENL*V>K4*RH/  ( MU*1 0  ^  -  6 ) 

RH=W*H/(2>K(W+H))  m 
V=Q/AF  ■ 
DPCV=DENL*VMAX" 2/2 

f=.043/NRE^.2  "For  5000<NRE<200000  - 
HYHD=DENL*V^2/2  ■ 
DPT=f*LT/(4>t=RH)*HYHD  ■ 
DPNAD=  ( CN+ ( 1 -E )  *  ( 1  - AR^  2 ) )  >KHYHD 

AR=4WH/(PI(  )=^D^2)  ■ 
RTW=RHOW*10^-8*2*TW/(H*LT)  ■ 
PIE=VL5t:IT 

BMIN=Ml*IT/(2*(M2+M3)>t=Q)  ■ 
"B=BMIN  ■ 
EEH=POH/PIE*100  ■ 
MT=NCELL5t<2*(H+W)>KTW*LT>KDENRTW 

ML=NCELL*(H*W*LT+VOLND)*DENRTL  ■ 
"ML=2*(H>t<W>KLT+4>t<VOLND)>»<DENRTL  ■ 
S=PI(  )=^(H^2+W^2)"(l/2)*(LNDI/COS(PI(  )/12)+LNDO/COS(PI()/24) ) 

■■  S=  ( PI  ( ) >KD+2*  ( W+H ) ) *  ( LNDI+LNDO ) /2  ■ 
■■S=PI()*(D/2+R2)*(LND^2+(D/2-R2)"2)''(l/2)  | 
SELECTW=STEP(H,W) 

(1-SELECTW)*LNDI=(D-H)/(2’KTAN(PI(  )/12)  )*(1-SELECTW)  - 
(l-SELECTW)*LNDO=(D-H)/(2*TAN(PI(  )/24)  )>K(1-SELECTW)  ■ 
SELECTW*LNDI  =  ( D-W ) / ( 2  *TAN ( PI ( ) /1 2 ) ) ^SELECTW  ■ 
SELECTW^LNDO=  ( D-W ) /( 2*TAN ( PI  ( ) /24 ) ) >KSELECTW 

LND=LNDI+LNDO  ■ 
■■LND=(D/2-R2)/TAN(PI()/12)  ■ 
MND=NCELL*S*TW>KDENRTW 

■■MND=4*S*TW>KDENRTW  ■ 
VOLNDI= (PI { ) *D^2/4+(W*H) ) *LNDI/2  ■ 
VOLNDO=  ( PI  ()  ^D^'  2/4+  ( W>KH ) )  *LNDO/2 

VOLND=VOLNDI+VOLNDO  _ 
■■VOLND=PI(  )XLND/3*(  (D/2)^2+D/25KR2+R2^2)  I 
’■R2=(W^2+H^2)^(l/2)/2  ■ 
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MATERIALS  PROPERTIES 


TWC=TL-273 

TWK=TL 

DENW=8 . 6304-1 . 895E-4*TWK-1 . 23E-85KTWK‘'2 
DENRTW=8 . 6304-1 . 895E-4>KTRM-1 . 23E-8*TRM^2 

RHOW=16  337+4. 224E-2*TWC-4.922E-6*TWC''2-3.941E-10>t<TWC''3  "For  Nb-Zr  in  uOhm-cm 


"DENW=7 . 68+0 .  E-2*TWK+0 .  E-65t<TWK" 2 
"DENRTW=7. 68+0.E-2>HTRM+0.E-6*TRM"2 
"RHOW=69 . 6+8. 78E-2*TWC+1 . 0-6*TWC'^2 


"For  SS  in  MicroOhm-cm 
"For  SS  in  MicroOhm-cm 
"For  SS  in  MicroOhm-cm 


RHOL=13. 735+2. 9256E-2*TL-2.46E-6>t=TL^2  "For  Li  in  micro-Ohm-cm 

DENL=.  5593-1. 133E-4*TL+1.58E-8*TL^2  "For  Li  in  Mg/in^3 

DENRTL=.  5593-1. 133E-4;t;TRM+1.58E-8>*'TRM"2  "For  Li  in  Mg/m '‘3 

MU=1 . 8739-4. 667E-3*TL+4. 933E-6>f<TL''2-l.  867E-9>KTL''3  "For  Li  in  centipoise 

"RH0Mo=-0. 183+2. 0321E-2>t'TL  "For  Mo  in  microOhm-cm 

RHOMo=(-8. 144267+.  21863752*TL-6.88E-10*TL^2+2.3868E-8>t<TL^3-5. 28E-12>t:TL^4)/10 

"For  Mo  in  microOhm-cm 

DENMo=10. 2698-1. 056E-4*TL-4.08E-8*TL"2  "For  Mo  in  Mg/m^3 

KMo=(0. 3602-1. 142E-4*TL+2.050E-8*TL"2)*4. 1868  "For  Mo  in  W/cm-K 


"RHOL=36. 768-1. 6212E-2*TL+7. 1741E-5*TL^2 
"DENL=. 9390-2. 426E-4*TL 
"DENRML=. 9390-2. 426E-4*TRML 

"MU=1 . 1914-2 . 7431E-3*TL+2 . 5463E-65HTL"2-8 . 492E-105icTL" 


"For  NaK  in  micro-Ohm-cm 
"For  NaK  in  Mg/m^3 
"For  NaK  in  Mg/m'' 3 
3  "For  NaK  in  centipoise 


"For  Cu  in  nOhm-m 
"For  ZrB2  in  uOhm-cm 
"For  ZrB2  in  Mg/m'' 3 

"SINGLE-PASS  PUMP  SECTION 


"DPDL=DPDLE  "SP 
"LL=IT>«RB/(RB+RL)>t;RL  "SP 
"VLW=IT*RB/(RB+RL)*RTW  "SP 
"VL=VLL+VLW+CEMF+VBS  "SP 
"POH=DPDL*Q  "SP 
•■LB=2*(2*HFB)  "SP 
"MB=2*(TB5t:HFB)>t<LT*DENSB  "SP 


RHOB=RHOMo 

RBS=  RH0B*LB/(LT*TB)*10''-8 

VBS=IT*RBS 

PIT=PIE+PLMP 

MASS =MLM+MB+MT+ML+MND 

SPMEMP=MASS/PIT 

MEMPTEG=SPMTEG*PIE 

MLMTEG= SPMTEG^PLMP 

EMASS=MASS+MEMPTEG+MLMTEG 


RHOCU=16. 73+0. 068* (TEMPB-293) 
RH0ZrB2=-5. 81+.0357*TL 
DENZrB2=6.08*. 969 


"DOUBLE-PASS  PUMP  SECTION 
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DPDL=DPDLE/2 

"DP 

VLL=2*IT*RB/(RB+RL)*RL 

"DP 

VLW=  2*1 T*RB/ ( RB+RL ) *RTW 

"DP 

VL=VLL+VLW+VC+VBS 

"DP 

VC=2*CEMF 

"DP" 

POH=2*DPDL*Q 

"DP 

LB=TB+TI+H 

"DP 

DENSB=DENMo 

MB= ( 2*H+3*TI+4*TW) *TB*LT*DENSB 

"DP 

"LINEAR  CORELESS  MAGNET  CALCULATIONS  -  7/22/91 


UP=(l-STEP(fX,fLl)) 

PUMP=STEP(fX,  fLl  )>t<STEP(fL2,  fX) 
DOWN=(l-STEP(fL2,fX) ) 

Ll=X/2-LT/2 
L2=X/2+LT/2 
Ll=fLl*X 
L2=fL2*X 

"UPSTREAM  OF  PUMPING  SECTION 

UP>i<PHIl=ATAN(  (fLl-fX)>KX/(fW*WFB)  )>KUP 
UP>)<PHI2=ATAN(  (fL2-fX)*X/(fW>t=WFB)  )5KUP 
UP*PHI3=ATAN( (fLl-fX)*X/( (l-fW)*WFB) )*UP 
UP*PHI 4 = ATAN ( ( f  L2 - f X ) *X/ ( ( 1 - f W ) *WFB ) ) *UP 
UP>KHFPR1  =  (HF2-HF1  )*UP 
UP*HFPR2= (HF4-HF3 ) *OP 

"PUMPING  SECTION  REGION 

PUMP*PHI1=ATAN( (fX-fLl)*X/(fW*WFB) )*PUMP 
PUMP*PHI2=ATAN( (f L2-fX)*X/(fW*WFB) )*PUMP 
PUMP*PHI3=ATAN(  (fX-fLl)*X/(  (l-fW)’KWFB)  )*PUMP 
PUMP*PHI4=ATAN( (fL2-fX)»X/( (l-fW)*WFB) )*PUMP 
PUMP^HFPRl = ( HF1+HF2 ) *PUMP 
PUMP*HFPR2= ( HF3+HF4 ) *PUMP 
HF1=SIN(PHI1) 

HF2=SIN(PHI2) 

HF3=SIN(PHI3) 

HF4=SIN(PHI4) 

Rl=(  (fW^^WFB)^2+(fH*HFB)^2)^.5 
R2=(  ((l-fW)*WFB)'^2+(fH>KHFB)^2)".  5 
HFPl=MU0*I/(4>*:pi(  )*R1) 

HFP2=MU0*I/(4*pi( )*R2) 

HAXVl=fW*WFB/Rl 

HAXV2=fW*WFB/R2 

HFP-HFPl*HFPRl+HFP2’t:HFPR2 

HFPAX=HFP1*HFPR1»«HAXV1+HFP2*HFPR2>KHAXV2 

B=NLM*HFPAX 

"DOWNTSTREAM  OF  PUMPING  SECTION 

DOWN*PHIl=ATAN(  (fX-f  LI  )>KX/(fW>KWFB)  )>KDOWN 
DOWN*PHI2=ATAN(  (fX-fL2)>*:X/(fWWFB)  )  WWN 


"STEP=1  if  fX=,  >fLl  and  0  otherwise 
"CASE=1  for  fLl<fX<fL2 
"STEP=1  if  fX<fl2 
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D0WN*PHI3=ATAN( (fX-fLl)*X/{ ( l-fW)*WFB) )*DOWN 
D0WN>i<PHI4=ATAN(  (fX-fL2)*X/(  (l-fW)*WFE)  )*DOWN 
D0WN*HFPR1  =  ( HFl  -HF2 )  Jt^DOWN 
DOWNOTPR2= (HF3-HF4 )^DOWN 


I  =  IBB 

RHOBB=RHOMo 
CD= IBB/ABB 

WFB=W+ETA+2*  ( TW+TB+2>KTI ) 
HFB=TI/2+H+TI 
ABB=ETA*HFB 
DENBB=DENMo 

RBB=RHOBB* 1 0 ^ - S^LBB/ABB 

LBB=  2*(LT+WFB) 

VBB=LBB*ABB 

MLM=NLM>KVBB*DENBB 

VLBBP=IBB*RBBP 

RBBP=RBB 

’'VLBBS=IBB*RBBS 

■'RBBS=2*RBB 

PLMP = NLM* I BB ^  2 *RBBP 

"PLMS=IBB^2*RBBS 


MMMMMMMMMMMMMMMMMMMM  UNIT  SHEET 


FromDDDDD 

ToDDDDDDD 

Multiply  ByDD  Add  OffsetDDD 

T 

mT 

1000 

psi 

Pa 

6895 

m 

mm 

1000 

m^2 

mm^2 

1000000 

in''3 

inm^3 

1E9 

m^3/s 

litre/s 

1000 

T“in''2/ s 

mVolt 

1000 

litre/s 

gpm 

15.85 

T-A-m 

N 

N/m^2 

Pa 

Volt 

Ohm-A 

Mg/m-s 

cp 

1000000 

Mg/m-s^2 

Pa 

1000 

T-m^2/s 

Ohm-A 

Pa-m"3/s 

W 

kPa 

T-A-m/m'"2 

Mg 

kg 

1000 

Ohm -A 

mVolt 

1000 

kA 

A 

1000 

Volt 

mVolt 

1000 

in''2/s^2 

Pa -m^ 3 /Mg 

1000 

m''2 

cm"' 2 

10000 

m 

cm 

100 

m^2/s^2 

Pa-m^3/Mg 

1000 

Volt-A 

W 

cm 

mm 

10 

V 

mV 

1000 

kW 

Volt-A 

1000 

kPa 

Pa 

1000 

m 

cm 

100 

m^3 

Gm^3 

1000000 

T 

mT 

1000 

A/cm'"  2 

A/m^2 

10000 

Mg/W 

g/w 

1000000 
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APPENDIX  D 


LINEAR  MAGNET  MODEL  FOR  CONDDCTORS 
OF  LARGE  CROSS-SECTION  SIZE 
MAGFLDDT.TK 


VARIABLE  SHEET 

St  InputDDDD  NameDDD  OutputDDD  UnitDDDDD  ComraentDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 

FILE:  MAGFLDDT.TK  52  kB 
MAGNETIC  FIELD  CALCULATIONS 
7-22-91 


I 

3250 

A 

.5 

fW 

- 

4.2 

WFB 

m 

.2 

fLl 

- 

.8 

fL2 

- 

.048 

LI 

m 

L2 

.192 

m 

.5 

fX 

- 

.24 

X 

m 

UP 

0 

- 

PUMP 

1 

“ 

DOWN 

0 

PHIl 

.03427229 

rad 

PHI  2 

.03427229 

PHI  3 

.03427229 

PHI  4 

.03427229 

HFl 

.03426558 

HF2 

.03426558 

HF3 

.03426558 

HF4 

.03426558 

HFPRl 

.06853116 

HFPR2 

.06853116 

HFPl 

. 1548066 

mT 

HFP2 

.1548066 

mT 

1.257E-6 

MUO 

T-m/A 

HAXVl 

1 

HAXV2 

1 

HFP 

.02121815 

mT 

HFPAX 

.02121815 

mT 

0 

fH 

- 

HFB 

.05 

m 

ho 

3 

10 

w 

3 

h 

.5 

fw 

ThetaA 

.5404195 

ThetaB 

.5404195 

Beta 

2.0607537 

Ro 

5.6666667 

FOM 

.3 

Bfl 

.14705882 

Bf2 

.14705882 

BB 

.29411765 

dd 

5.8309519 

ddd 

5.8309519 

nnl 

.85749293 

nn2 

.85749293 

Field  Current 
Fraction  of  WFB 

Center-center  span  between  field  buses 

Fractional  location  of  LI 

Fractional  location  of  L2 

Inlet  end  of  Field  busbars 

Discharge  end  of  Field  busbars 

Fraction  of  total  pump  length 

Total  pump  length 

Upstream  region  indicator 

Pumping  section  region  indicator 

Downstream  region  indicator 

Angular  integration  limit 

Angular  integration  limit 

Angular  integration  limit 

Angular  integration  limit 

Integrated  flux  fraction 

Integrated  flux  fraction 

Integrated  flux  fraction 

Integrated  flux  fraction 

Total  integrated  fractions,  bar  2 

Total  integrated  fractions,  bar  1 

Magnetic  Flux  from  bar  1 

Magnetic  Flux  from  bar  2 

Absolute  permeability 

Axial  flux  vector  for  HFPl 

Axial  flux  vector  for  HFP2 

Total  magnetic  flux 

Total  axial  magnetic  flux 

Fraction  of  height  between  busbars 
Height  of  the  busbars 

EQUIPOTENTIAL  CIRCLE  CALCULATIONS 

Height  of  Po  above  busbar  line  at  w/2 
Length  of  busbar  line 
Height  of  P  above  busbar  line  at  fw*w 
Fraction  of  the  busbar  line  length 
Angle-radial  vector  A  to  busbar  line 
Angle-radial  vector  B  to  busbar  line 
Angle  between  radial  vectors  A  &  B 
Radius  of  the  equipotential  circle 
This  is  the  ratio  of  h/w 
Off  bisector  factor  for  busbar  1 
Off  bisector  factor  for  busbar  2 
Sum  of  off-bisector  factors  above 
distance  from  busbar  1 
distance  from  busbar  2 
Cosine  factor  for  busbar  1 
Cosine  factor  for  busbar  2 
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L  25 
L  4 
L  130 
L 


IPB  130 

ETA  5 

ABB 

Wo 

CD 

IBB  3250 

RHOBB 

RHOMo 

VBB 

LT 

DENBB 

DENMo 

MBB 

RBB 

LBB 

NCELL 

VLBBP 

RBBP 

VLBBS 

RBBS 

PBBP 

PBBS 


A  Current  per  busbar  pairs 

cm  Cross-section  side  of  large  BBs 

cm''2  Cross-sectional  area  of  large  busbars 

cm  W  of  inner  most  busbar  pairs 

A/cm" 2  Current  density 

A 


.5  fWl 


Fraction  of  WFB 


fHll  2 

fH12  1 

fH13  0 

fH14  1 

fH15  2 


cm 

Fraction 

of 

cm 

Fraction 

of 

cm 

Fraction 

of 

cm 

Fraction 

of 

cm 

Fraction 

of 

height  between  busbars 
height  between  busbars 
height  between  busbars 
height  between  busbars 
height  between  busbars 


2 

1 

0 

1 

2 


NFl 

NF2 

NF3 

NF4 

NF5 


Height  displacement  unit  from  busbar 
Height  displacement  unit  from  busbar 
Height  displacement  unit  from  busbar 
Height  displacement  unit  from  busbar 
Height  displacement  unit  from  busbar 


L  WFBl  5 

L  WFB2  7 

L  WFB3  9 

L  WFB4  11 

L  WFB5  13 


cm 

Center-center 

cm 

Center-center 

cm 

Center-center 

cm 

Center-center 

cm 

Center-center 

span  between  field  buse;^ 
span  between  field  buse* 
span  between  field  buse" 
span  between  field  buses 
span  between  field  buse* 


PHIll  1.2365983 
PHI21  1.2365983 
PHI31  1.2365983 
PHI41  1.2365983 
HFll  .94467366 
HF21  .94467366 
HF31  .94467366 
HF41  .94467366 
HFPRll  1.8893473 
HFPR21  1.8893473 


Angular  integration  limit 
Angular  integration  limit 
Angular  integration  limit 
Angular  integration  limit 
Integrated  flux  fraction 
Integrated  flux  fraction 
Integrated  flux  fraction 
Integrated  flux  fraction 
Total  integrated  fractions,  bar  1 
Total  integrated  fractions,  bar  2 


HFPlll  .51973456  mT 
HFP211  .51973456  mT 
HAXVlll  .99920096 
HAXV211  .99920096 


Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
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HFPAXll 

1.962349 

mT 

HFPl 12 

.52004619 

mT 

HFP212 

.52004619 

mT 

HAXV112 

.99980006 

HAXV212 

.99980006 

HFPAX12 

1.9647028 

mT 

HFPl 13 

.52015019 

mT 

HFP213 

.52015019 

mT 

HAXV113 

1 

HAXV213 

1 

HFPAX13 

1.9654887 

mT 

HFPl 14 

.52004619 

mT 

HFP214 

.52004619 

mT 

HAXV114 

.99980006 

HAXV214 

.99980006 

HFPAX14 

1.9647028 

mT 

HFPl 15 

.51973456 

mT 

HFP215 

.51973456 

mT 

HAXV115 

.99920096 

HAXV215 

.99920096 

HFPAX15 

1.962349 

mT 

HFPAXl 

9.8195923 

mT 

PHI  12 

1.1183214 

rad 

PHI22 

1.1183214 

PHI  32 

1.1183214 

PHI  4  2 

1.1183214 

HF12 

.89936785 

HF22 

.89936785 

HF32 

.89936785 

HF42 

.89936785 

HFPR12 

1.7987357 

HFPR22 

1.7987357 

HFPl 21 

.37138429 

mT 

HFP221 

.37138429 

mT 

HAXV121 

.99959209 

HAXV221 

.99959209 

HFPAX21 

1.3354994 

mT 

HFP122 

.37149794 

mT 

HFP222 

.37149794 

mT 

HAXV122 

.99989797 

HAXV222 

.99989797 

HFPAX22 

1.3363169 

mT 

HFPl 2 3 

.37153585 

mT 

HFP223 

.37153585 

mT 

HAXV123 

1 

HAXV223 

1 

HFPAX23 

1.3365896 

mT 

Total  axial  magnetic  flux 

Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 

Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 

Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 

Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 

Total  axial  magnetic  flux  of  layer  1 

LAYER  2 

Angular  integration  limit 
Angular  integration  limit 
Angular  integration  limit 
Angular  integration  limit 
Integrated  flux  fraction 
Integrated  flux  fraction 
Integrated  flux  fraction 
Integrated  flux  fraction 
Total  integrated  fractions,  bar  1 
Total  integrated  fractions,  bar  2 

Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 

Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 

Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 
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HFP124 

HFP224 

HAXV124 

HAXV224 

HFPAX24 

HFP125 

HFP225 

HAXV125 

HAXV225 

HFPAX25 

HFPAX2 


PHI  13 

PHI23 

PHI  33 

PHI  4  3 

HF13 

HF23 

HF33 

HF43 

HFPR13 

HFPR23 

HFP131 

HFP231 

HAXV131 

HAXV231 

HFPAX31 

HFP132 

HFP232 

HAXV132 

HAXV232 

HFPAX32 

HFP133 

HFP233 

HAXV133 

HAXV233 

HFPAX33 

HFP134 

HFP234 

HAXV134 

HAXV234 

HFPAX34 

HFP135 

HFP235 

HAXV135 

HAXV235 

HFPAX35 


.37149794  inT 
.37149794  mT 
.99989797 
.99989797 
1.3363169  niT 


Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 


.37138429  mT 
.37138429  mT 
.99959209 
.99959209 
1.3354994  mT 


Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 


6.6802221  mT 


Total  axial  magnetic  flux  of  layer  2 


LAYER  3 


1.012197  rad 

1.012197 

1.012197 

1.012197 

.8479983 

.8479983 

.8479983 

.8479983 

1.6959966 

1.6959966 


Angular  integration  limit 
Angular  integration  limit 
Angular  integration  limit 
Angular  integration  limit 
Integrated  flux  fraction 
Integrated  flux  fraction 
Integrated  flux  fraction 
Integrated  flux  fraction 
Total  integrated  fractions,  bar  1 
Total  integrated  fractions,  bar  2 


.288901  mT 
.288901  mT 
.99975318 
.99975318 
.97970836  mT 


Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 


.28895449  mT 
.28895449  mT 
.99993828 
.99993828 
.98007117  mT 


Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 


.28897233  mT 
.28897233  mT 
1 
1 

1.1202196  mT 


Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 


.28895449  mT 
.28895449  mT 
.99993828 
.99993828 
.98007117  mT 


Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 


.288901  mT 
.288901  mT 
.99975318 
.99975318 
.97970836  mT 


Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 


HFPAX3 


PHI14 

PHI24 

PHI  34 

PHI44 

HF14 

HF24 

HF34 

HF44 

HFPR14 

HFPR24 

HFP141 

HFP241 

HAXV141 

HAXV241 

HFPAX41 

HFP142 

HFP242 

HAXV142 

HAXV242 

HFPAX42 

HFP143 

HFP243 

HAXV143 

HAXV243 

HFPAX43 

HFP144 

HFP244 

HAXV144 

HAXV244 

HFPAX44 

HFP145 

HFP245 

HAXV145 

HAXV245 

HFPAX45 

HFPAX4 


PHI  15 
PHI  25 
PHI  3  5 
PHI  4  5 
HF15 
HF25 
HF35 
HF45 
HFPR15 


5.0397787 

mT 

Total  axial  magnetic  flux  of  layer  3 

LAYER  4 

.91846543 

rad 

Angular  integration  limit 

.91846543 

Angular  integration  limit 

.91846543 

Angular  integration  limit 

.91846543 

Angular  integration  limit 

.79467101 

Integrated  flux  fraction 

.79467101 

Integrated  flux  fraction 

.79467101 

Integrated  flux  fraction 

.79467101 

Integrated  flux  fraction 

1.589342 

Total  integrated  fractions,  bar  1 

1.589342 

Total  integrated  fractions,  bar  2 

.23639283 

niT 

Magnetic  Flux  from  bar  1 

.23639283 

mT 

Magnetic  Flux  from  bar  2 

.99983475 

Axial  flux  vector  for  HFPl 

.99983475 

Axial  flux  vector  for  HFP2 

.75129395 

mT 

Total  axial  magnetic  flux 

.23642213 

mT 

Magnetic  Flux  from  bar  1 

.23642213 

mT 

Magnetic  Flux  from  bar  2 

.99995868 

Axial  flux  vector  for  HFPl 

.99995868 

Axial  flux  vector  for  HFP2 

.75148021 

mT 

Total  axial  magnetic  flux 

.2364319 

mT 

Magnetic  Flux  from  bar  1 

.2364319 

mT 

Magnetic  Flux  from  bar  2 

1 

Axial  flux  vector  for  HFPl 

1 

Axial  flux  vector  for  HFP2 

.75154231 

mT 

Total  axial  magnetic  flux 

.23642213 

mT 

Magnetic  Flux  from  bar  1 

.23642213 

mT 

Magnetic  Flux  from  bar  2 

.99995868 

Axial  flux  vector  for  HFPl 

.99995868 

Axial  flux  vector  for  HFP2 

.75148021 

mT 

Total  axial  magnetic  flux 

.23639283 

mT 

Magnetic  Flux  from  bar  1 

.23639283 

mT 

Magnetic  Flux  from  bar  2 

.99983475 

Axial  flux  vector  for  HFPl 

.99983475 

Axial  flux  vector  for  HFP2 

.75129395 

mT 

Total  axial  magnetic  flux 

3.7570906 

mT 

Total  axial  magnetic  flux  of  layer  4 

LAYER  5 

.83644866 

rad 

Angular  integration  limit 

.83644866 

Angular  integration  limit 

.83644866 

Angular  integration  limit 

.83644866 

Angular  integration  limit 

.74226804 

Integrated  flux  fraction 

.74226804 

Integrated  flux  fraction 

.74226804 

Integrated  flux  fraction 

.74226804 

Integrated  flux  fraction 

1.4845361 

Total  integrated  fractions,  bar  1 
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HFPR25  1.4845361 

HFP151  .20003409  inT 
HFP251  .20003409  mT 
HAXV151  .99988168 
HAXV251  .99988168 
HFPAX51  .59384538  mT 

HFP152  .20005184  mT 
HFP252  .20005184  mT 
HAXV152  .99997042 
HAXV252  .99997042 
HFPAX52  .59395079  mT 

HFP153  .20005776  mT 
HFP253  .20005776  mT 
HAXV153  1 
HAXV253  1 

HFPAX53  .59398594  mT 

HFP154  .20005184  mT 
HFP254  .20005184  mT 
HAXV154  .99997042 
HAXV254  .99997042 
HFPAX54  .59395079  mT 

HFP155  .20003409  mT 
HFP255  .20003409  mT 
HAXV155  .99988168 
HAXV255  .99988168 
HFPAX55  .59384538  mT 

HFPAX5  2.9695783  mT 

L  HFPAXT  28.266262  mT 

BAXll  9.8195923  mT 
BAX12  16.499814  mT 
BAX13  21.539593  mT 
BAX14  25.296684  mT 
BAX15  28.266262  mT 

LAYICON  34.739621 
LAY2CON  23.6332 
LAY3CON  17.829661 
LAY4CON  13.291784 
LAY5CON  10.505734 


Total  integrated  fractions,  bar  2 

Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 

Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 

Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 

Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 

Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Axial  flux  vector  for  HFPl 
Axial  flux  vector  for  HFP2 
Total  axial  magnetic  flux 

Total  axial  magnetic  flux  of  layer  5 

Total  axial  magnetic  flux  -  all  layer 

Axial  magnetic  flux  -  Layer  1 
Axial  magnetic  flux  -  Layers  1  and  2 

Axial  magnetic  flux,  Layers  1  thru  3 

Axial  magnetic  flux,  Layers  1  thru  4 

Axial  magnetic  flux.  Layers  1  thru  5 
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RULE  SHEET 

IS  Ru 1 eDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 
"MAGNETIC  FLUX  CALCULATIONS  FOR  A  LINEAR  CORELESS  MAGNET 
"7/22/91 


"MAGNETIC  FIELD  CALCULATIONS 

UP=(l-STEP(fX,fLl))  "STEP=1  if  fX=,  >fLl  and  0  otherwise 

PUMP=STEP(fX,fLl)*STEP(fL2,fX)  "CASE=1  for  fLl<fX<fL2 

DOWN=(l-STEP(fL2,fX))  "STEP=1  if  fX<fl2 

Ll=fLl*X 
L2=fL2*X 

"UPSTREAM  OF  PUMPING  SECTION 

UP*PHI  1  =ATAN ( ( f  LI -f X )  *X/ ( f  W*WFB ) ) >KUP 
UP*PHI 2  =  ATAN  ( ( f  L2 -f X )  ♦X/  ( f W>tcWFB ) )  >*<UP 
UP*PHI3=ATAN(  (fLl-fX)>KX/(  (l-fW)*WFB)  )>>:UP 
UP=t'PHI4=ATAN(  (fL2-fX)*X/(  (l-fW)*WFB)  )>KUP 
UP*HFPR1 = ( HF2-HF1 ) *UP 
UP*HFPR2= (HF4-HF3 )*DP 

"PUMPING  SECTION  REGION 

PUMP>KPHI1=ATAN(  (fX-fLl)*X/(fW*WFB)  )*PUMP 
PUMP*PHI2=ATAN( (fL2-fX)*X/(fW*WFB) )*PUMP 
PUMP*PHI3=ATAN( (fX-f LI )*X/( ( l-fW)*WFB) )*PUMP 
PUMP*PHI4=ATAN( (fL2-fX)*X/( ( l-fW)*WFB) )*PUMP 
PUMP>*‘HFPR1  =  ( HFl  +HF2 )  *PUMP 
PUMP*HFPR2= ( HF3+HF4 ) *PUMP 
HF1=SIN(PHI1) 

HF2=SIN(PHI2) 

HF3=SIN(PHI3) 

HF4=SIN(PHI4) 

HFPl  =MU0  *  I  /  ( 4  2)«pi  ( )  *  ( ( f  W*WFB ) "  2  +  ( f  H*HFB ) '' 2 ) .  5 ) 

HFP2=MU0*I/(4*pi( )*( ( (l-fW)*WFB)^2+{fH*HFB)^2)^. 5) 
HAXVl=COS(ATAN(fH*HFB/(fW*WFB) ) ) 

HAXV2=COS  ( ATAN  ( f  H*HFB/  ( ( 1  -f  W )  >KWFB ) ) ) 

HFPAX=HFP 1 *HFPR1 *HAXV 1 +HFP2 *HFPR2  *HAXV2 
HFP=HFP1*HFPR1+HFP2*HFPR2 

"DOWNTSTREAM  OF  PUMPING  SECTION 

DOWN*PH  1 1 = AT  AN  ((fX-fLl)*X/(  f  W>t^WFB ) )  *DOWN 
DOWN*PHI2=ATAN( (fX-fL2)*X/(fW*WFB) )*DOWN 
DOWN*PHI3=ATAN((fX-fLl)*X/( (i-fW)*WFB) )*DOWN 
DOWN*PHI4=ATAN( (fX-fL2)*X/( (l-fW)*WFB) )*DOWN 
DOWN*HFPRl=  (HF1-HF2 )  >i'DOWN 
DOWN*HFPR2= (HF3-HF4 ) *DOWN 

"EQUIPOTENTIAL  ELECTRIC  FIELD  AND  MAGNETIC  VECTOR  CALCULATIONS 

2=««atan(ho/(w/2)  )=atan(h/(fw2t:w)  )+atan(h/(  (  l-fw)*w) ) 

ThetaA=atan(h/(fw*w) ) 

ThetaB=atan(h/( ( l-fw)*w) ) 

Beta=pi ( ) - ( ThetaA+ThetaB ) 

Ro=W'2/(8*ho)+ho/2 
FOM=h/w 
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*  Bf  l=cos(atan(h/(fw*w)  ) )/(  (h''2+ (fw*w)  ^^2  ) '"  ( 1/2) ) 

*  Bf  2  =  cos(atan(h/(  (1-fw) *V7) ) )/( (h"2+ ( ( 1-fw)  *w) ''2 ) "  ( 1/2 ) ) 

*  BB=Bfl+Bf2 

*  dd=(h''2+(fw*w)^2)^(l/2) 

*  ddd=(h^2+((l-fw)*w)^2)''(l/2) 

*  nnl=cos (ThetaA) 

*  nn2=cos (ThetaB) 

"LAEGE  CONDUCTOR  CALCULATIONS 

*  IPB= IBB/25 

*  I=IBB 

^  RHOBB=RHOMo 

*  CD= IBB/ABB 
ABB=ETA^2 

*  HFB=ETA 

*  DENBB=DENMo 

*  RBB=RHOBB>K10''-8*LBB/ABB 

*  LBB=  2*(LT+WFB) 

*  VBB=2*LBB*ABB 

*  MBB=NCELL*VBB*DENBB 

*  VLBBP=IBB*RBBP 

*  RBBP=RBB 
■■VLBBS=IBB*RBBS 
■■RBBS=2^RBB 

*  PBBP=2*IBB^2»RBBP 
■’PBBS=IBB''2*RBBS 

*  fHll=NFl*ETA/5 

*  fH12=NF2*ETA/5 

*  fH13=NF3*ETA/5 

*  fH14=NF4*ETA/5 

*  fH15=NF5*ETA/5 

■■WFB=Wo+0.2 

*  WFB1=Wo+1*ETA/5 

*  WFB2=Wo+3*ETA/5 

*  WFB3=Wo+5>t<ETA/5 

*  WFB4=Wo+7*ETA/5 

*  WFB5=Wo+9*ETA/5 

"LAYER  I 

*  PUMP*PHIll=ATAN(  (fX-fLl)*X/(fWl*WFBl)  )>^PUMP 

*  PUMP*PHI21=ATAN(  (fL2-fX)*X/(fWl>KWFBl )  )>KPUMP 

*  PUMP*PHI31=ATAN(  (fX-fLl)*X/(  (1-fWl)  WBl)  )*PUMP 

*  PUMP»<PHI41=ATAN((fL2-fX)W(  (l-fWl)*WFBl)  )5t<PUMP 

*  PUMP*HFPR11=(HF11+HF21)*PUMP 

*  PUMP*HFPR21=(HF31+HF41)*PUMP 

*  HF11=SIN(PHI11) 

*  HF21=SIN(PHI21) 

*  HF31=SIN(PHI31) 

*  HF41=SIN(PHI41) 

*  HFPl  1 1=MU0*IPB/ ( 4*pi  ( ) * ( ( fWl^WFBl )  ''2+  ( f HI l^HFB) " 2)^.5) 

HFP2 1 1  =MU0* IPB/  ( 4*pi  ( ) >K  ( ( ( 1  -f W1 )  *WFB1 )  ^ 2+  ( f  HI  1>KHFB )  ^ 2 ) .  5 ) 

*  HAXVlll=COS(ATAN(fHll>cHFB/(fWl*WFBl ) ) ) 

*  HAXV211=COS(ATAN(fHll>KHFB/(  ( 1-fWl  )*WFB1 ) ) ) 
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*  HFPAX11=HFP111*HFPR11*HAXV111+HFP211*HFPR21*HAXV211 


*  HFP112=MU0*IPB/(4*pi(  )>t^(  (fWl^WFBl )  ^2+(fH125t'HFB) ''2 )  ^  .  5) 

*  HFP212=MU05i<IPB/(45Xpi(  )*(  ( (1-fWl  )*WFBlK2+(fH12*HFB)  "2)'' .  5) 

*  HAXV 1 1 2 = COS  ( ATAN  ( f  HI  2 5*^HFB/  ( f  W 1  *WFB  1 ) ) ) 

*  HAXV212=COS(ATAN(fH12*HFB/( (l-fWl)*WFBl) ) ) 

*  HFPAXl 2  =HFP1 1 2*HFPR1 1 ^HAXV 1 1 2+HFP2 1 2*HFPR21 *HAXV  212 

>«  HFPl  1 3=MU0*IPB/ ( 4>i^pi  ( )  *  ( ( f  Wl^WFBl ) " 2+  ( f HI 3*HFB ) " 2 ) " .  5 ) 

*  HFP213=MU0*IPB/(4>KpiO*(  ( ( 1-fWl  )*WFB1 )  "2+(fH13*HFB)  "2)  .  5 ) 

*  HAXV113=COS(ATAN(fH13*HFB/(fWl*WFBl) )) 

*  HAXV213=COS(ATAN(fH13*HFB/(  (l-fWl)>t^WFBl) ) ) 

*  HFPAXl  3 =HFP1 1 3 *HFPR1 1  *HAXV  1 1 3+HFP2 1 3*HFPR2 1  J^HAXV 213 

*  HFP114=MU0*IPB/(4=t<pi(  )*(  (fWl*WFBl)^2+(fH14*HFB)^2)^.5) 

*  HFP214=MU0*IPB/(4*pi(  )*(  ( (l-fWl)5»‘WFBl)'‘2+(fH14>KHFB)"2)^.  5) 

*  HAXV114=COS(ATAN(fH14*HFB/(fWl*WFBl))) 

*  HAXV214=COS(ATAN(fH14*HFB/((l-fWl)*WFBl)) ) 

*  HFPAXl  4  =HFP  114 *HFPR1 1  ’KHAXV 114 +HFP2 1 4 *HFPR2 1  *HAX V 214 

*  HFPl  1 5=MOO*  IPB/  ( 4*pi  ( )  *  ( ( f  W1  *WFB1 )  ^ 2+  ( f  HI  5*HFB ) ''  2 ) "  .  5 ) 

*  HFP215=MU0*IPB/(4*pi(  )*(  ( ( 1-fWl  )*WFB1 ) "2+ (fHlSJt^HFB) "2) ^  .  5) 

*  HAXV115=COS{ATAN(fH15*HFB/(fWl>KWFBl) ) ) 

*  HAXV215=COS(ATAN(fH155t«HFB/(  (l-fWl)>f<WFBl) ) ) 

*  HFPAXl  5  =HFP1 1 5>t<HFPRl  1*HAXV1 1 5+HFP2 1 5*HFPR2 1*HAXV2 1 5 

*  HFPAXl =HFPAX1 1 +HFPAX1 2+HFPAXl 3+HFPAXl 4  +HFPAX1 5 
"LAYER  2 

PUMP*PHI12=ATAN( (fX-fLl)*X/(fWl*WFB2) )*PUMP 
PUMP*PHI22=ATAN( (fL2-fX)*X/(fWl*WFB2) )*PUMP 
PUMP*PHI32=ATAN( (fX-fLl)*X/( (l-fWl)*WFB2) )*PUMP 
PUMP*PHI42=ATAN( (f L2-fX)*X/( ( 1-fWl )*WFB2 ) )*PUMP 
PUMP*HFPR1 2  =  ( HFl 2  +HF2  2 ) *PUMP 
PUMP*HFPR2 2 = ( HF3 2+HF4 2 ) *PUMP 
HF12=SIN(PHI12) 

HF22=SIN(PHI22) 

HF32=SIN(PHI32) 

HF42=SIN(PHI42) 

HFP121=MU0>t‘IPB/(45t^pi(  )*(  (fWl*WFB2)''2+(fHll*HFB)"2)".5) 
HFP221=MU0*IPB/(4*pi( )*( ( ( 1-fWl )*WFB2) ^2+(fHll*HFB) ^2)" . 5) 
HAXV121=COS(ATAN(fHll^HFB/(fWl>*!WFB2) ) ) 

HAXV22 1 =COS ( ATAN ( f HI 1 *HFB/ ( ( 1 -f W1 ) *WFB2 ) ) ) 

HFPAX2 1  =HFP1 2 1  *HFPR1 2  *HAXV1 2 1 +HFP2  2 1 3»<HFPR2  2*HAXV  221 


*  HFP122=MU0*IPB/{4*pi(  )*(  (fWl*WFB2)  ^2+(fH12*HFB)''2) .  5) 

*  HFP222=MU0*IPB/(4*pi(  )=*=(  ( ( 1-fWl  )*WFB2)  ^2+(fH12*HFB)  "2)^  .  5 ) 

*  HAXV122=COS(ATAN(fH12’KHFB/(fWl*WFB2))) 

*  HAXV222=COS(ATAN(fH12*HFB/(  (l-fWl)><<WFB2)) ) 

*  HFPAX22=HFP122*HFPR125t'HAXV122+HFP222*HFPR22*HAXV222 


*  HFPl 23=MU0*IPB/ ( 4>*^pi  ( ) *  ( ( f W1=XWFB2 )  ^2+ ( f  HI S^t^HFB)  ^ 2)^.5) 

*  HFP223=MU0* IPB/  ( 4>Kpi  ( )  >K  ( ( ( 1  -f  W1 )  5»«WFB2 ) ''  2+  ( f  HI  3*HFB )  ^  2 )  ^  .  5 ) 
^  HAXV123=COS(ATAN(fH13>»<HFB/(fWl>tcWFB2) ) ) 

*  HAXV223=COS(ATAN(fH13>KHFB/(  ( 1-fWl  )5)cWFB2) ) ) 

*  HFPAX23=HFP123*HFPR12*HAXV123+HFP223*HFPR22*HAXV223 
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*  HFP124=MU0*IPB/(4*pi()*(  (fWl*WFB2)"2+(fH14>t^HFB)^2)".  5) 

*  HFP224=MU0*IPB/(4*pi(  )*(  (  (1-f  W1 ) >tcWFB2 )  ^2+ ( f H14:^cHFB)  ^2 ) .  5 ) 

*  HAXV 1 24 = COS  ( AT  AN  ( f  HI  4  *HFB/  ( f  W1  WFB2 ) ) ) 

*  HAXV224=COS(ATAN(fH14*HFB/( (1-fWl )*WFB2) ) ) 

*  HFPAX24=HFP124*HFPR12*HAXV124+HFP224*HFPR22>t:HAXV224 

*  HFP125=MU0>*^IPB/(4>Kpi(  )>»=(  (fWl5»cWFB2)^2+(fH15*HFB)^2)''.5) 

*  HFP2 25=MU0*  IPB/  ( 4>Kpi  ( )  >K  ( ( ( 1  -f  W1 )  ^WFB2 )  '^  2+  ( f  H 1  S^KHFB ) "  2 ) "  .  5 ) 

*  HAXV125=COS(ATAN(fH15*HFB/(fWl*WFB2))) 

=t<  HAXV225=COS(ATAN(fH15*HFB/(  (l-fWl)*WFB2) ) ) 

*  HFPAX25=HFP125’t<HFPR12*HAXV125+HFP225*HFPR22>tcHAXV225 

*  HFPAX2=HFPAX21+HFPAX22+HFPAX23+HFPAX24+HFPAX25 
"LAYER  3 

*  PUMP^PHI  1 3=ATAN ( ( f X-f  LI )  *X/ ( f  W1 WB3 ) )  *PUMP 

*  PUMP*PHI23=ATAN( (fL2-fX)*X/(fWl*WFB3) )*PUMP 

*  PUMP*PHI33=ATAN(  (fX-fLl)*X/(  (l-fWl)^WFB3)  )5»^PUMP 

*  PUMP^PHI43=ATAN(  (f  L2-fX)  W(  ( 1-fWl  )>KWFB3)  )*PUMP 
PUMP*HFPR13=(HF13+HF23)*PUMP 

*  PUMP5i<HFPR23=(HF33+HF43)*PUMP 

*  HF13=SIN(PHI13) 

*  HF23=SIN(PHI23) 

*  HF33=SIN{PHI33) 

*  HF43=SIN(PHI43) 

*  HFP131=MU0*IPB/(4>Kpi(  )5»c(  (fWl*WFB3)^2+(fHll*HFB)^2)^.5) 

*  HFP231=MU0*IPB/(4*pi(  )5k(  ( ( 1-fWl  )*WFB3) ''2+(fHll*HFB) ^2) ^ .  5) 

*  HAXVl 31 =COS ( ATAN ( f HI 1*HFB/ ( f W1*WFB3 ) ) ) 

*■  HAXV231=COS(ATAN(fHll>i‘HFB/(  ( 1-fWl  )*WFB3 ) ) ) 

*  HFPAX31=HFP131*HFPR135*«HAXV131+HFP231*HFPR23*HAXV231 

*  HFP132=MU0*IPB/(4*pi( )*( (fWl*WFB3)^2+(fH12*HFB)"2)^ . 5) 

*  HFP232=MOO*IPB/(4*pi(  )3tc(  ( (l-fWl)*WFB3)"2+(fH12*HFB)^2)^.5) 

*  HAXV132=COS(ATAN(fH12*HFB/(fWl>KWFB3) ) ) 

*  HAXV232=COS(ATAN(fH12OTB/((l-fWl)»^WFB3))  ) 

*  HFPAX32=HFP132*HFPR13*HAXV132+HFP232*HFPR23*HAXV232 

*  HFP133=MU0*IPB/(4*pi(  )*(  (fWl5tcWFB3)''2+(fH13*HFB)''2)^.5) 

*  HFP233=MU0*IPB/(4*pi(  )*(((l-fWl)>KWFB3)^2+(fH13*HFB)^2)^.5) 

*  HAXV133=COS(ATAN(fH13*HFB/(fWl*WFB3) ) ) 

*  HAXV233=COS(ATAN(fH13>XHFB/(  (l-fWl)*WFB3) ) ) 

*  HFPAX33=HFP133OTPR13*HAXV133+HFP223*HFPR23*HAXV233 

*  HFP134=MU0*IPB/(4*pi(  )♦(  (fWl*WFB3 ) '^2+ (fH14>t'HFB) ^2 ) ^  .  5) 

*  HFP234=MU05t:IPB/( 4*pi  ( ) 5*=  ( ( (1-fWl ) >KWFB3 )  ^2+  ( f H145t^HFB)  ^2 )  ^  .  5 ) 

*  HAXV134=COS(ATAN(fH14*HFB/(fWl*WFB3) ) ) 

*■  HAXV234=COS(ATAN(fH14>i«HFB/(  ( 1-fWl )  WB3) ) ) 

*  HFPAX34=HFP134*HFPR135i<HAXV134+HFP234*HFPR23>KHAXV234 

*  HFPl 35=MU05*'IPB/(  4*pi  ( ) * ( ( f W1  *WFB3 )  ^ 2+  ( f  HI 5*HFB )  ^ 2 )  ^  .  5 ) 

*  HFP235=MU0*IPB/ ( 4*pi  ()*(((  1-f W1 ) >KWFB3 )  ^2+  ( f HI S^J'HFB)  ^2 ) .  5 ) 

*  HAXV135=COS(ATAN(fH15*HFB/(fWl*WFB3) ) ) 

*  HAXV235=COS(ATAN(fH15*HFB/(  (l-fWl)5t<WFB3) ) ) 

*  HFPAX3  5 = HFP 135  *HFPR  1 3  5*^HAX V 1 3  5 +HFP 2  3  5 *HFPR  2  3  *H AX V  235 
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*  HFPAX3=HFPAX31+HFPAX32+HFPAX33+HFPAX34+HFPAX35 


■  "LAYER  4 

I*  PUMP*PHI14=ATAN( (fX-fLl)*X/(fWl^WFB4) )*PUMP 
>K  PUMP*PHI24=ATAN(  (fL2-fX)  W(fWl*WFB4)  )>KPUMP 

*  PUMP^i^PH  134  =  ATAN  ((fX-fLl)*X/({l-fWl)  *WFB4 ) )  >*^PUMP 

I*  PUMP*PHI44=ATAN(  (fL2-fX)*X/(  ( 1-fWl  )*WFB4 )  )>t<PUMP 
*  PUMP*HFPR14=(HF14+HF24)*PUMP 
*  PUMP*HFPR24=(HF34+HF44)’t<PUMP 
*  HF14=SIN(PHI14) 

I*  HF24=SIN(PHI24) 

*  HF34=SIN(PHI34) 

*  HF44=SIN(PHI44) 

I  *  HFP141=MU0*IPB/(4*pi(  )*(  (fWl*WFB4)''2+(fHll*HFB)"2)'‘ .  5) 

■  *  HFP241=MU0*IPB/(4>t=pi()=^(((l-fWl)>^WFB4)-2+(fHll*HFB)"2)'^.5) 

I*  HAXV141=C0S(ATAN(fHll*HFB/(fWl*WFB4) ) ) 

*  HAXV241=COS(ATAN(fHll*HFB/(  (l-fWl)>t«WFB4) ) ) 

*  HFPAX4 1  =HFP1 4 1  *HFPR1 4 ’KRAX V 1 4 1 +HFP2 4 1  *HFPR2 4 *HAXV 241 

I*  HFP142=MU0*IPB/(  4*pi  ()  *  ( ( f  W1*WFB4 )  ''2+  ( f  H12*HFB)  "2 ) "  .  5 ) 

*  HFP242=MU0*IPB/(4*pi(  )*(  ( (l-fWl)5t=WFB4)"2+(fH12>t^HFB)"2)".  5) 
HAXV142=COS  ( ATAN ( f  HI 2>KHFB/ ( f Wl>t<WFB4 ) ) ) 

■  *  HAXV242=COS(ATAN(fH12*HFB/( (l-fWl)*WFB4) ) ) 

I  ♦  HFPAX42=HFP142*HFPR14*HAXV142+HFP242*HFPR24>t^HAXV242 

*  HFP143=MU0*IPB/(4*pi(  )*(  (fWl*WFB4)''2+(fH13*HFB)"2)'‘.5) 

I*  HFP243=MU0*IPB/(4*pin*(  ( (l-fWl)*WFB4)"2+(fH13*HFB)"2)"  .  5) 
*  HAXV143=COS(ATAN(fH13*HFB/(fWl*WFB4) ) ) 

*  HAXV243=COS(ATAN(fH13*HFB/( (l-fWl)*WFB4) ) ) 

■  *  HFPAX43=HFP143*HFPR14*HAXV143+HFP243*HFPR24*HAXV243 

*  *  HFP144=MU0*IPB/(4*pi(  )*(  (fWl*WFB4)''2+(fH14*HFB)''2)".  5) 

I*  HFP244=MU0*IPB/(4*pi(  )*(  ( (l-fWl)*WFB4)"2+(fH14*HFB)"2)".  5) 
*  HAXV144 =COS ( ATAN ( f H14*HFB/ ( f W1 *WFB4 ) ) ) 

*  HAXV244=COS(ATAN(fH14*HFB/(  (l-fWl)5KWFB4) ) ) 

HFPAX44=HFP144*HFPR14*HAXV144+HFP244*HFPR24*HAXV244 

I  *  HFP145=MU0>t=IPB/(  4*pi  ()  *  ( ( fWl*WFB4  )  ^2+  ( f  H15*HFB)  "2  )  ^  .  5 ) 

*  HFP245=MU0*IPB/(4*pi(  )*(  ( (l-fWl)>KWFB4)"2+(fH15*HFB)"2)‘'.  5) 

I*  HAXV145=COS(ATAN(fH15*HFB/(fWl*WFB4) ) ) 

*  HAXV245=COS(ATAN(fH15>t<HFB/(  (l-fWl)=t<WFB4) ) ) 

*  HFPAX4 5 =HFP1 4 5*HFPR1  A*EAXV  1 4 5+HFP 245 >t:HFPR2 4 *HAXV 245 

1 HFPAX4=HFPAX41+HFPAX42+HFPAX43+HFPAX44+HFPAX45 

"LAYER  5 

I  *  PUMP*PHI15=ATAN( (fX-fLl)*X/(fWl*WFB5) )*PUMP 
*  PUMP*PH  125  =  AT  AN  ( ( f  L2  -  f  X )  *X/  ( f  W 1  *  WFB5 ) )  >KPUMP 

I*  PUMP>kPHI35=ATAN(  (fX-fLl)*X/(  (l-fWl)*WFB5)  )5t«PUMP 
*  PUMP*PHI45=ATAN(  (fL2-fX)*X/(  (l-fWl)>KWFB5)  )*PUMP 
*  PUMP>t^HFPR15=(HF15+HF25)*PUMP 

I*  PUMP>t=HFPR25=(HF35+HF45)5>:PUMP 
HF15=SIN(PHI15) 

*  HF25=SIN(PHI25) 

*  HF35=SIN(PHI35) 


Dll 


*  HF45=SIN(PHI45) 


5K  HFP151=MU0*IPB/(4>Kpi(  )>*<(  (fWl*WFB5)^2+(fHll>t<HFB)^2)^.5) 

*  HFP251=MU0*IPB/(4*pi(  )*(  ( ( 1-fWl  )*WFB5)  ^2+(fHll*HFB) ''2)  ^  .  5) 

*  HAXV151=COS(ATAN(fHll>HHFB/(fWlWB5) ) ) 

*  HAXV251=COS(ATAN(fHll5t^HFB/(  (1-fWl  )*WFB5) ) ) 

HFPAX5 1  =HFP1 5 1  *HFPR1 5*HAXV1 5 1 +HFP2  5 1  >t<HFPR2  5 *HAXV2  5 1 

*■  HFP152=MU0*IPB/(4>*:pi(  )*(  (fWl*WFB5)"2+(fH12*HFB)^2)^.  5) 

*  HFP252=MU0*IPB/( 4^pi  ()*(((  1-fWl )  ^♦'‘WFBS )  ^2  +  ( f  H12*HFB)  ^2 )  ^  .  5 ) 

*  HAXV152=COS(ATAN(fH12*HFB/(fWl>KWFB5) ) ) 

*  HAXV252=COS(ATAN(fH12*HFB/( (l-fWl)*WFB5))) 

*  HFPAX52=HFP152*HFPR15*HAXV152+HFP252*HFPR25*HAXV252 

*  HFP153=MU0*IPB/(43t'pi(  )*(  (fWl*WFB5)"2+(fH13*HFB)"2)^.5) 

*  HFP2 53=MU0* IPB/ ( 4*pi  ()*(((  1  -fWl )  *WFB5 )  ^ 2+  ( f  HI 3*HFB )  ^ 2 )  .  5 ) 

*  HAXV153=COS(ATAN(fH13>t«HFB/(fWl*WFB5) ) ) 

*  HAXV253=COS(ATAN(fH13=XHFB/(  (l-fWl)>>;WFB5) ) ) 

*  HFPAX53=HFP153>*‘HFPR15>KHAXV153+HFP253OTPR255tcHAXV253 

*  HFP154=MU0*IPB/(4*pi(  ):<<(  (fWl*WFB5)''2+(fH14>»:HFB)^2)''.  5) 

*  HFP254=MU05KlPB/(4>t«pi()*(  ((l-fWl)WB5)''2+(fH14*HFB)^2)'-.5) 

*  HAXV154=COS(ATAN(fH145KHFB/(fWl*WFB5))) 

*  HAXV254=COS(ATAN(fH14>t:HFB/(  (l-fWl)>»<WFB5) ) ) 

*  HFPAX54=HFP154*HFPR15>kHAXV154+HFP2545KHFPR25*HAXV254 

*  HFP155=MU0*IPB/(4*pi(  )*(  (fWl*WFB5)^2+(fH155KHFB)^2)^.5) 

*  HFP255=MU0*IPB/(4*pi{)>K((  (l-fWl)*WFB5)^2+(fH15*HFB)^2)''.5) 
>t:  HAXVl  55 =COS  ( AT  AN  ( f  HI  5*HFB/  ( f  W1  >*^WFB5 ) ) ) 

*  HAXV255=COS(ATAN(fH15*HFB/(  (l-fWl)>»«WFB5) )  ) 

5*^  HFPAX5  5 = HFP 1 5  5  sKHFPR  1 5  >I«HAX V 155 +HFP 255 *HFPR2  5  *HAX V  255 

*  HFPAX5=HFPAX51+HFPAX52+HFPAX53+HFPAX54+HFPAX55 

*  BAX11=HFPAX1 

*■  BAX12=BAX11+HFPAX2 

*  BAX13=BAX12+HFPAX3 

*  BAX14=BAX13+HFPAX4 

*  BAX15=BAX14+HFPAX5 

*  LAyiCONTR=HFPAXl/HFPAXT*100 

*  LAY2CONTR=HFPAX2/HFPAXT*100 

*  LAy3CONTR=HFPAX3/HFPAXT*100 

*  LAy4CONTR=HFPAX4/HFPAXT*100 

*  LAy5CONTR=HFPAX5/HFPAXT*100 

*  HFPAXT=HFPAX1+HFPAX2+HFPAX3+HFPAX4+HFPAX5 
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TABLE  SHEET 

NameDDDDDDDDDD  TitleDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 
HFP_2_8_91  HL,  L=12,  W= . 6  DISTRIBUTION  OF  MAGNETIC  FIELD 

H  ANGLE_2_8_91  HL,  L=12,  W=.6  MAGNETIC  FIELD  INTEGRATION  ANGLES 

H_SIN_2_8_91  HL,  L=12,  W=.6  MAGNETIC  FIELD  INTEGRATION  FRACTIONS 

HL_2_8_91  HL,  MAGNETIC  FIELD  ALONG  DUCT  (fW=.5) 

HW_2_8_91  HW,  MAGNETIC  FIELD  ACROSS  DUCT  (fL=.5) 

BH_2_11_91  BH,  MAGNETIC  FIELD  VS  HEIGHT  BETWEEN  BUSBARS 

CIRCCALC  EQUIPOTENTIAL  CIRCLE  CALCULATIONS 

HW_7_15_91  HW,  MAGNETIC  FIELD  ACROSS  DUCT  (fL=.21) 

LARGE_LM  AXIAL  FLUX  OF  LAYERS  VS  HEIGHT 

LARGE.LMl  AXIAL  FLUX  VS  CROSS-SECTIONAL  AREA  OF  BUSBARS 

LARGE_LM2  CONTRIBUTIONS  OF  LAYERS  TO  MAGNETIC  FLUX 


APPENDIX  E 


DESIGN  MODEL 
VCLCP 


VARIABLE  SHEET 

St  InputDDDD  NameDDD  OutputDDD  UnitDDDDD  CommentDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 

FILE:VCLCP.TK  37  kB 


.125 

DPDLE 

psi 

2 

QE 

gpm 

DPDL 

.861875 

kPa 

TDPDL 

kPa 

Q 

.12618297 

litre/s 

DPDV 

395.28887 

kPa 

DPMB 

220.43099 

kPa 

DPMF 

173.82796 

kPa 

DPHY 

.16804477 

kPa 

DPT 

.02522894 

Pa 

DPNAD 

.14281583 

Pa 

B 

3041.8864 

mT 

BOPT 

1902.7648 

mT 

IT 

1.2510426 

kA 

IL 

1.0026123 

kA 

lA 

1.2510426 

kA 

2 

IF 

kA 

6.35 

H 

mm 

12.7 

W 

mm 

35 

LT 

mm 

V 

1.5646719 

m/s 

VMAX 

m/s 

Ml 

129.61002 

1/m 

M2 

188793020 

l/0hm-m^2 

M3 

148878825 

l/0hm-m^2 

M4 

1 

- 

AF 

80.645 

mm^2 

IC 

.46015419 

kA 

RB 

108113.25 

nOhm 

RL 

23247.805 

nOhm 

RHOL 

406.83658 

nOhm-m 

RA 

26788.63 

nOhm 

RWP 

308034.09 

nOhm 

RHOW 

786.94841 

nOhm-m 

RWS 

3540.8252 

nOhm 

RFR 

166578.76 

nOhm 

RFB 

.12903451 

Ohm 

RHOFB 

24.5 

nOhm-m 

.000006 

AFB 

m''2 

RAB 

38582.677 

nOhm 

RHOAB 

24.5 

nOhm-m 

RCELL 

111273.31 

nOhm 

RTP 

45902.002 

nOhm 

RHOTP 

786.94841 

nOhm-m 

1 

FX 

_ 

.5 

TW 

mm 

NRE 

20989.073 

- 

RH 

2.1166667 

mm 

DPCV 

Pa 

.0000001 

EPSOD 

VERIFICATION  CORELESS  LINEAR 

CONDUCTION  PUMP  (VCLCP)  PROGRAM 
DESIGN  AND  PERFORMANCE  PROGRAM 
3/28/91 

Delivered  pressure  from  cell 
Volumetric  flowrate 
Delivered  pressure  from  cell 
Trial  DPDL 
Volumetric  flowrate 

Differential  pressure  rise  developed 

Pressure  loss  from  magnetic  braking 

Pressure  loss  from  magnetic  fringe 

Pressure  loss  from  fluid  friction 

Pressure  loss  in  throat  flow 

Pressure  loss  in  nozzle  and  diffuser 

Magnetic  flux  noraml  to  Q  and  I 

Optimum  magnetic  flux 

Total  pump  current 

Electric  current  in  fluid 

Armature  busbar  current 

Field  busbar  current 

Inside  height  of  throat  duct 

Inside  width  of  throat  duct 

Length  of  throat 

Velocity  of  liquid  in  throat 

Velocity  limited  by  cavitation 

First  term  coefficient 

Second  term  coefficient 

Third  term  coefficient 

Fourth  term  coefficient 

Throat  cross-sectional  flow  area 

Counter  current  from  CEMF 

Electrical  resistance  bypassing  fluid 

Electrical  resistance  across  W 

Resistivity  of  fluid 

Electrical  resistance  armature  current 
Parallel  Resistance  of  throat  walls 
Resistivity  of  duct  wall 
Series  Resistance  of  throat  wall 
Resistance  of  fringe  path 
Resistance  of  field  busbar(coil) 
Resistivity  of  field  busbar  material 
Cross-section  area  of  coil  wires 
Resistance  of  armature  busbar 
Resistivity  -  armature  busbar  material 
FBB  to  FBB  Resistance  of  pump  cell 
Resistance  of  tap  pin 
Resisitivity  of  Tap  Pin 
Empirical  factor  from  MSA  calcs 
Thickness  of  throat  walls 
Reynolds  Number  in  the  throat 
Hydraulic  radius  in  the  throat 
Pressure  margin  for  cavitation  limit 
Ratio  e/D 


El 


LAMINAR 

0 

TURBLNT 

1 

f 

.00587584 

- 

HYHD 

1038.6606 

Pa 

.05 

CN 

- 

.8 

E 

- 

.75 

AR 

- 

D 

11.700735 

mm 

PIE 

202194.35 

W 

VC 

60.446338 

mVolt 

VLL 

23.308535 

mVolt 

VWS 

3.5500749 

mVolt 

VA 

26.85861 

mVolt 

VAB 

48.268574 

mVolt 

VTP 

57.425361 

mVolt 

VFB 

161.42768 

Volt 

VCELL 

192.99888 

mVolt 

VL 

161.62067 

Volt 

EEH 

5.3787E-5 

% 

POH 

.10875394 

W 

VOLND 

mm'' 3 

VOLNDI 

mm '"3 

VOLNDO 

mm''3 

S 

mm^2 

R2 

mm 

LNDI 

mm 

LNDO 

mm 

LND 

mm 

373 

TL 

K 

373 

TDW 

K 

MND 

kg 

ML 

kg 

MT 

kg 

MB 

kg 

MASS 

kg 

TI 

mm 

295 

TRML 

K 

295 

TRM 

K 

TWK 

373 

K 

CEMF 

60.446338 

mV 

I 

2000 

A 

.5 

fW 

- 

.02 

WFB 

m 

.25 

fLl 

- 

.75 

fL2 

- 

LI 

.06 

m 

L2 

.18 

m 

.5 

fX 

- 

.24 

X 

m 

UP 

0 

- 

PUMP 

1 

- 

DOWN 

0 

- 

PHIl 

1.4056476 

rad 

PHI  2 

1.4056476 

1  =  the  flow  IS  laminar,  0  =  IS  NOT 
1  =  the  flow  IS  transient  or  turbulen 
Fluid  friction  factor  in  throat 
Hydraulic  velocity  head  in  throat 
Discharge  coefficient  of  nozzles 
Diffuser  recovery  coefficient 
Area  ratio  throat  to  pipe 
Diameter  of  circuit  piping 
Input  power  to  EMP  (electrical) 
Counter  EMF  of  fluid  moving  across  B 
Voltage  drop  across  liquid 
Voltage  drop  across  wall 
Voltage  drop  outside  of  duct 
Voltage  drop  across  armature  buses 
Voltage  drop  across  tap  pin 
Voltage  drop  across  field  buses 
Voltage  drop  across  pump  cell 
"Load  or  terminal  voltage  of  TEG  or 
Electrical  to  hydraulic  efficiency 
Output  power  (hydraulic) 

Volume  of  nozzle  and  diffuser 
Volume  of  nozzle 
Volume  of  diffuser 

Surface  area  of  nozzle  and  diffuser 
Equivalent  radius  of  throat  area 
Length  of  nozzle 
Length  of  diffuser 
Length  of  nozzle  and  of  diffuser 
Temperature  of  fluid 
Temperature  of  the  walls 
Mass  of  nozzles  and  diffusers 
Mass  of  fluid 
Mass  of  pump  throat 
Mass  of  busbar 
Total  mass  of  EMP 

Thickness  of  electrical  insulation 


Room  temperature 
Temperaure  of  wall  material 


MAGNETIC  FIELD  CALCULATIONS 
2-13-91 
Field  Current 
Fraction  of  WFB 

Center-center  span  between  field  buse 
Fractional  location  of  LI 
Fractional  location  of  L2 
Inlet  end  of  Field  busbars 
Discharge  end  of  Field  busbars 
Fraction  of  total  pump  length 
Total  pump  length 
Upstream  region  indicator 
Pumping  section  region  indicator 
Downstream  region  indicator 
Angular  integration  limit 
Angular  integration  limit 
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PHI  3 

1.4056476 

PHI4 

1.4056476 

HFl 

.98639392 

HF2 

. 98639392 

HF3 

.98639392 

HF4 

.98639392 

HFPRl 

1,9727878 

HFPR2 

1.9727878 

HFPl 

20.005776 

mT 

HFP2 

20.005776 

mT 

1.257E-6 

MUO 

T-m/A 

N 

38.536938 

turns 

.01 

d 

m 

HFP 

78.934305 

mT 

.35 

LFB 

m 

.00635 

FEW 

m 

.0635 

LAB 

m 

.00635 

TAB 

m 

.000508 

LTP 

m 

ATP 

8.709E-10 

m^2 

.00333 

DOTP 

m 

0 

DITP 

m 

0 

fH 

* 

.0127 

HFB 

m 

373 

TCU 

K 

Angular  integration  limit 
Angular  integration  limit 
Integrated  flux  fraction 
Integrated  flux  fraction 
Integrated  flux  fraction 
Integrated  flux  fraction 
Total  integrated  fractions,  bar  2 
Total  integrated  fractions,  bar  1 
Magnetic  Flux  from  bar  1 
Magnetic  Flux  from  bar  2 
Absolute  Permeability 
Number  of  coil  turns 

Midpoint  between  long  sides  of  coils 
Total  magnetic  flux 
Length  of  field  busbar 
Thickness  of  field  busbar 
Length  of  amature  busbar 
Thickness  of  armature  busbar 
Length  of  Tap  Pin 
Cross-sectional  area  of  Tap  Pin 
Outside  diameter  of  Tap  Pin 
Inside  diameter  of  Tap  Pin 
Fraction  of  height  between  busbars 
Height  of  the  busbars 
Temperature 

MATERIALS  PROPERTIES 


STAINLESS  STEEL 


DENRTW 

Mg/m'' 3 

RT  density  of  wall  material 

RHOSS 

786.94841 

nOhm-m 

Resistivity  of  300  series 

DENW 

7.68 

Mg/m^3 

Density  of  wall  material 

NAK-78  ALLOY 

PV 

Pa 

Vapor  Pressure 

K 

W/m-K 

Thermal  Conductivity 

RHO 

DEN 

Mg/m^3 

Density 

RH078 

406.83658 

nOhm-m 

Resisitivity  of  NAK-78 

MU 

.53554943 

cp 

Viscosity  of  fluid 

DENL 

.8485102 

Mg/m'' 3 

Density  of  fluid 

DENRTL 

Mg/m ^3 

RT  density  of  fluid 

DENRML 

Mg/m^3 

RT  density  of  fluid 

COPPER 

DENSE 

Mg/m'' 3 

Density  of  busbar 

RHOCU 

24.5 

nOhm-m 

Resistivity 

E3 


MMMMMMMMMMMMMMMMMMMM  RULE  SHEET 

S  Ru 1 eDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDj 
"VCLCP  DESIGN  AND  PERFORMANCE  PROGRAM  ■ 

"04-1-91 


*  QE=15850*Q 

*  DPDV=M1>kB5*«IL 

*  DPMB=M2*B^2*Q 

*  DPMF=M35»<B''2*Q 

*  DPHY=(DPT+DPNAD)5t:1000 

*  DPDL=DPDV-DPMB-DPMF-DPHY 

" TDPDL=DPDV-DPMB-DPMF-DPHY 
>x  CEMF=B*W*Q/AF 

*  AF=W*H 

*  IC=CEMF/(RB+RL) 

*  IL=IA*RB/(RB+RA) 

" I A= I T*RFB/ ( RFB+RCELL ) 

IA=IT 

*  RL=RHOL>KW/(H*LT) 

*  RH0L=RH078 

*  RA=RL+RWS 

*  RWS=RHOW*2*TW/(H*LT) 

RHOW=RHOSS 

RB=RFR*RWP/(RFR+RWP) 

*  RFR=2.6*FX*RHOL/H 

*  RWP=RHOW*(W+2*TW)/(TW=)<2*LT) 

*  RFB=RHOFB>KN/(AFB)*(2*LFB+6>KWFB) 

*  RHOFB=RHOCU 
RCELL=RTP+RAB+RA 

*  RAB=RHOAB*LAB/(H*TAB) 

*  RHOABrRHOCU 

>K  RTP=RHOTP*LTP/ATP 

*  RHOTP=RHOSS 

*  ATP=PI()*(DOTP^2-DITP^2)/4 

*  M1=RB/(RB+RL)*1/H 

*  M2=1/(RB+RL)*1/H^2 

*  M3=1/(H^2*RFR) 

*  M4  =  l 

NRE=DENL*V*4*RH/(MU*10^-6) 

*  RH=W>KH/(2*(W+H)) 

V=Q/AF 

*  DPCV=DENL*VMAX^2/2 

*  LAMINAR=STEP(2000,NRE) 

>K  TURBLNT=1 -LAMINAR 

*  LAMINAR*! =64/NRE*LAMINAR 

"(l-LAMINAR)*l/sqrt(f )=-2*log(EPSOD/3.7+2.51/(NRE*sqrt(f )))*(1-LAMINAR) 

*  (TURBLNT)*f=.043/NRE^. 2*(TURBLNT)  "For  5000<NRE<200 

*  HYHD=DENL*V^2/2 

*  DPT=f*LT/(4*RH)*HYHD 

*  DPNAD=(CN+(1-E)*(1-AR^2) )*HYHD 

*  AR=4*W*H/(PI()*D^2) 

*  PIE=VL*IA 

*  B0PT=M1*IA/(2*(M2+M3)*Q) 

••B=BOPT 

*  EEH=POH/PIE*100 
"MT=2*(H+W)*TW*LT*DENRTW 
'•ML=2*  ( H*W*LT+VOLND )  *DENRTL 
"ML=2*(H*W*LT+4*VOLND)*DENRTL 


E4 


"S=(PI( )*D+2*(W+H) )*(LNDI+LND0)/2 
■■S=PI()*(D/2+R2)*(LND^2+(D/2-R2)^2)^(l/2) 

’'LNDI=(H-D)/(2*TAN(PI()/12)) 

•' LNDO=  ( H-D )  / ( 2>l'TAN ( PI  ( )  /24 ) ) 

”LND=LNDI+LNDO 
" LND= ( D/2 -R2 ) /TAN ( PI ( ) /I 2 ) 

■’MND=2*S*TW*DENRTW 
"MND=4*S*TW*DENRTW 
"  VOLNDI  =  ( PI  ( ) 2/4+W*H )  5KLNDI/2 
" VOLNDO= ( PI { ) 2 /4  +W*H ) *LNDO/2 
" VOLND= VOLNDI+ VOLNDO 

■'VOLND=PI(  )*LND/3*(  (D/2)"2+D/2>kR2+R2^2) 

•’R2=(W''2+H^2)^(l/2)/2 

*  TWC=TL-273 

*  TWK=TL 

"MATERIALS  PROPERTIES 
"NB-ZR  ALLOY 

"DENW=8 . 6304-1 . 895E-4*TWK-1 . 23E-85t«TWK''2 
"DENRTW=8. 6304-1. 895E-4*TRM-1. 23E-8*TRM''2 

"RHOW=16. 337+4 . 224E-2*TWC-4 . 922E-6*TWC"2-3 . 941E-10*TWC"3  "raicroOhm-cm 
"316  STAINLESS  STEEL 

5(c  DENW=7.68+0.E-2*TWK+0.E-6*TWK^2  "Mg/m^3 

*  RHOSS=(47. 874+7. 1365*10''-2*TDW+3.02*10''-5>t'TDW'^2)*10''-8  "Ohm-m 

"LITHIUM 

"RHOL=13. 735+2. 9256E-2’KTL-2. 46E-6*TL''2  "micro-Ohm-cm 

"DENL= .  5593-1 . 133E-4*TL+1 . 58E-8>i<TL''2  "Mg/m '“3 

"DENRTL= . 5593-1 . 133E-4*TRM+1 . 58E-8*TRM"2  "Mg/m^3 

"MU=1 . 8739-4 . 667E-3*TL+4 . 933E-6*TL"2-1 . 867E-9*TL''3  "centipoise 

"MOLYBDENUM 

"RH0Mo=-0. 183+2. 0321E-2*TL  "microOhm-cm 

"DENMo=10 . 2698-1 . 056E-45t‘TL-4 . 08E-8>t«TL‘"2  "Mg/m''3 

"KMo= (0 . 3602-1 . 142E-4*TL+2 . 050E-8*TL"2 )*4 .1868  "W/cm-K 

"NAK-78 

"RHOL=36. 768-1. 6212E-2>i'TL+7.1741E-55t:TL"2  "micro-Ohm-cm 

*  RH078=(36.75-1.631*10''-2*TL+7.2*10''-55t<TL''2)*10"-8  "Ohm-m 

*  DENL=.  9390-2. 426E-4*TL  "Mg/m ''3 

"DENRML= .  9390-2 . 426E-4*TRML  "Mg/m ''3 

"MU=1. 1914-2. 7431E-3*TL+2. 54 63E-6*TL''2-8. 492E-10*TL''3  "centipoise 

*  MU=3. 437-1. 62*10^-2*TL+3. 1845i<10^-5>t;TL''2-2.833*10^-8>i<TL^3+9.375*10^-12*TL'^4 

"mPa-s 

"K=-6. 917+1. 872*10''-l*TL-4. 377*10^-4*TL^2+4. 858*10^-7*TL^3-2.08*10^-10*TL''4 

"W/m-K 

"PV=3. 1622*10''9/(EXP(10686/T) )  "Pa 

"COPPER 

*  RHOCU=(-.  96471+6. 454*10''-2*TCU+l*10"-5*TCU''2)*le-9  "Ohm-m 
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"SINGLE-PASS  PUMP  SECTION 

*  DPDL=DPDLE  "SP 

*  VLL=IL*RL  "SP 

VWSrIL’t'RWS  "SP 

*  VA=IL»RA  "SP 

VC=CEMF  "SP 

VCELL=VC+VLL+VWS+VAB+VTP  "SP 

VL=VCELL+VFB  "SP 

*  VTP=IA*RTP  "SP 

*  VFB=IA>KRFB  "SP 

*■  VAB=IA*RAB  "SP 

*  POH=DPDL*Q  "SP 

"LB=H/2+TI+3/2*TB+W  "SP 

"MB=  ( LB-TB/2+H/2 )  *TB*LT*r)ENSB  "SP 


RHOB=RHOCD 

"RBS=  RHOB*LB/(LT=t'TB)*10^-8 
"VBS=IT>t<RBS 
" MASS=MM+MB+MT+ML+MND 

"DOUBLE-PASS  PUMP  SECTION 


"DPDL=DPDLE/2  "DP 
"VLL=2>KIT*RB/(RB+RL)>KRL  "DP 
"VLW=2*IT*RB/(RB+RL)*RTW  "DP 
"VL=VLL+VLW+VC+VBS  "DP 
"VC=2’i^CEMF  "DP" 
”POH=2*DPDL*Q  "DP 
"LB=TB+TI+H  "DP 
"MB=(2>kH+35>:TI+4*TW)*TB*LT*DENSB  "DP 
"DENSB=DENMo 


"MAGNETIC  FLUX  CALCULATIONS  FOR  A  RECTANGULAR  COIL  ELECTROMAGNETIC  PUMP 
"2/7/91 

*  UP=(l-STEP(fX,fLl))  "STEP=1  if  fX=,  >fLl  and  0  otherwise 

*  PUMP=STEP(fX,fLl)*STEP(fL2,fX)  "CASE=1  for  fLl<fX<fL2 

*  DOWN=(l-STEP(fL2,fX))  "STEP=1  if  fX<fl2 

*  Ll=fLl*X 
L2=fL2*X 

"UPSTREAM  OF  PUMPING  SECTION 

*  UP5*^PHI1=ATAN(  (fLl-fX)>i<X/(fW*WFB)  )*UP 
^  UP^t^PHI  2 = ATAN  ( ( f  L2  -  f  X )  sf-X/  ( f  W*WFB ) ) 

*■  UP*PHI3=ATAN(  (fLl-fX)*X/(  ( l-fW)*WFB)  )*UP 

*  UPJt^PHI  4  =  ATAN  ( ( f  L2  -  f  X )  *X/  ( ( 1  -  f  W )  *WFB ) )  *UP 

*  UP>KHFPR1=(HF2-HF1)*UP 
UP5KHFPR2=(HF4-HF3)*DP 

"PUMPING  SECTION  REGION 

*  PUMP^PHI  1  =ATAN  ( ( f  X-f  LI )  *X/  ( f  W^KWFB ) )  *PUMP 

*  PUMP*PHI2=ATAN(  (f L2-fX)*X/(fW>KWFB)  )*PUMP 

*  PUMP*PHI3=ATAN(  (fX-fLl  )*X/(  ( l-fW)5HWFB)  )*PUMP 

*  PUMP*PHI4=ATAN(  (f L2-fX)*X/(  ( l-fW)*WFB)  )>t<PUMP 
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*  PUMP*HFPR1=(HF1+HF2)*PUMP 

*  PUMP*HFPR2  = ( HF3+HF4 ) *PUMP 

*  HF1=SIN(PHI1) 

*  HF2=SIN(PHI2) 

*  HF3=SIN(PHI3) 

*  HF4=SIN(PHI4) 

^  I-IF 

*  HFP1=MU0*I/(4*PI(  )5K(  (fW*WFB)"2+(fH>t=HFB)^2)^.5) 

*  HFP2=MU0*I/(4>KPI(  )*(  ( (l-fW)>KWFB)''2+(fH*HFB)"2)''.  5) 

*  HFP=HFPl*HFPRl+HFP2=t'HFPR2 

*  B=N*HFP 

"DOWNTSTREAM  OF  PUMPING  SECTION 

*  DOWN*PHIl=ATAN(  (fX-f LI )>KX/(fW*WFB)  )3KDOWN 

*  DOWN*PHI2=ATAN(  (fX-fL2)>*'X/(fW*WFB)  )*DOWN 

*  DOWN5»;PHI3=ATAN(  (fX-f LI  )>i<X/(  ( l-fW)*WFB)  )*DOWN 

*  DOWN>i^PHI4=ATAN(  (fX-fL2)>t‘X/(  (l-fW)*WFB)  )*DOWN 

*  D0WN5KHFPR1  =  (HF1-HF2)*D0WN 

*  DOWN*HFPR2=(HF3-HF4)*DOWN 


*  N=4/ ( HFPRl  +HFPR2 )  *  ( Ml / ( 2*  ( M2+M3 ) ) *PI  ( )  *d/  ( MUO^t^Q ) ) 
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MMMMMIWIMMMMMMMMMMMM  UNIT  SHEET 

FromDDDDD  ToDDDDDDD  Multiply  ByDD  Add  OffsetDDD 


T 

mT 

1000 

psi 

Pa 

6895 

N/m^2 

Pa 

Mg/in-s"2 

Pa 

1000 

Pa-m^3/s 

W 

kPa 

T-A-m/m''2 

kPa 

Pa 

1000 

in"2/s^2 

Pa-m'‘3/Mg 

1000 

kPa 

T-A-m/m^2 

cm 

mm 

10 

m"2 

cm'' 2 

10000 

m 

cm 

100 

m 

mm 

1000 

m'‘2 

mm ''2 

1000000 

m"3 

mm  ■'3 

1E9 

m''3/s 

litre/s 

1000 

litre/s 

gpm 

15.85 

T-m''2/s 

mVolt 

1000 

T-A-m 

N 

Volt 

Ohm-A 

T-m^2/s 

Ohm-A 

T-in^2/s 

Ohm-A 

Ohm-A 

mVolt 

1000 

Volt 

mVolt 

1000 

Mg/ra-s 

cp 

1000000 

Mg 

kg 

1000 

kg 

hg 

10 

Mg 

kg 

1000 

kA 

A 

1000 

kA 

hA 

10 

Volt-A 

kW 

.001 

V 

mV 

1000 

Pa-m''3/s 

W 

kW 

Volt-A 

1000 

W 

Volt-A 

Ohm-m 

nOhm-m 

1E9 

Ohm-m 

uOhra-m 

1000000 

Ohm-m 

mOhm-m 

1000 

Ohm 

uOhm 

1000000 

Ohm 

mOhm 

1000 

Ohm 

nOhm 

1E9 

FLUX  VS  VARIABLES 

1  FLUX 

10 

15 

20 

25 

30 

35 

EFF 

13.3833872 

25.1620972 

36.2249292 

45.2431519 

51.9858677 

56.7062036 

■  LOADV 

623.767694 

493.226541 

451.170387 

444.508234 

455.543081 

476.690499 

■  AMPS 

1750.5502 

1177.52522 

894.160246 

726.659284 

617.090326 

540.625369 

MASS 

79890.905 

711249.162 

712606.82 

713969.364 

715341.09 

716725.398 

_  POWER 

1091.93666 

580.786689 

403.418625 

323.006035 

281.111228 

257.710977 

I  WIDTH 

664.6 

664.6 

664.6 

664.6 

664.6 

664.6 

■  LENGTH 

11750 

11750 

11750 

11750 

11750 

11750 

HIEGHT 

221.5 

221.5 

221.5 

221.5 

221.5 

221.5 

m  MAGNET 

2748.8978 

4107.15498 

5464.81244 

6827.3565 

8199.08262 

9583.39045 

1  BUSBAR 

4246.71221 

4246.71221 

4246.71221 

4246.71221 

4246.71221 

4246.71221 

DIFFUSER 

52.9553838 

52.9553838 

52.9553838 

52.9553838 

52.9553838 

52.9553838 

■  LIQUID 

2753.07589 

2753.07589 

2753.07589 

2753.07589 

2753.07589 

2753.07589 

1  THROAT 

89.2637365 

89.2637365 

89.2637365 

89.2637365 

89.2637365 

89.2637365 

1 

FLUX  VS  VARIABLES 

■  FLUX 

40 

45 

50 

55 

60 

65 

1  EFF 

59.7845301 

61.5819715 

62.3982191 

62.4696543 

61.9794225 

61.0689519 

LOADV 

504.158274 

535.839621 

570.470467 

607.246405 

645.631161 

685.253469 

B  AMPS 

484.850414 

442.868793 

410.542506 

385.23828 

365.2006 

349.214108 

■  MASS 

718125.004 

719542.098 

720978.458 

722435.54 

723914.539 

725416.448 

■  POWER 

244.441348 

237.306646 

234.202375 

233.934561 

235.784887 

239.300179 

WIDTH 

664.6 

664.6 

664.6 

664.6 

664.6 

664.6 

■  LENGTH 

11750 

11750 

11750 

11750 

11750 

11750 

■  HIEGHT 

221.5 

221.5 

221.5 

221.5 

221.5 

221.5 

MAGNET 

10982.9966 

12400.0907 

13836.4511 

15293.5324 

16772.5319 

18274.441 

■  BUSBAR 

4246.71221 

4246.71221 

4246.71221 

4246.71221 

4246.71221 

4246.71221 

■  DIFFUSER 

52.9553838 

52.9553838 

52.9553838 

52.9553838 

52.9553838 

52.9553838 

LIQUID 

2753.07589 

2753.07589 

2753.07589 

2753.07589 

2753.07589 

2753.07589 

1  THROAT 

89.2637365 

89.2637365 

89.2637365 

89.2637365 

89.2637365 

89.2637365 

■  FLUX  VS 

VARIABLES 

*  FLUX 

70 

-  EFF 

59.8477639 

■  LOADV 

725.84814 

■amps 

336.410694 

MASS 

726942.092 

■  power 

244.183077 

■  WIDTH 

664.6 

LENGTH 

11750 

■  HIEGHT 

221.5 

■  MAGNET 

19800.0849 

•  BUSBAR 

4246.71221 

_  DIFFUSER 

52.9553838 

■  LIQUID 

2753.07589 

■  THROAT 

1 

89.2637365 

1 

1 

E9 

(s)  Screen  or  Printer:  Printer  25+/!j 

Display  options:  Screen  Printer 

TABLE:  FLUX 

Screen  or  Printer:  Printer  I 

Title:  FLUX  VS  VARIABLES  ■ 

Vertical  or  Horizontal:  Horizontal 

Row  Separator:  j 

Column  Separator:  | 

First  Element:  1 

Last  Element:  13  j 

ListDDDDD  Numeric  FormatDD  WidthDD  HeadingDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDl 


B 

10 

FLUX 

EEH 

10 

EFF 

VL 

10 

LOADV 

IT 

10 

AMPS 

MASS 

10 

MASS 

PIE 

10 

POWER 

W 

10 

WIDTH 

LT 

10 

LENGTH 

H 

10 

HIEGHT 

MM 

10 

MAGNET 

MB 

10 

BUSBAR 

MND 

10 

DIFFUSER 

ML 

10 

LIQUID 

FI  Help  F2  Cancel 

F5  Edit 

F9  Solve  /  ' 

Window  switch 
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APPENDIX  F 


RESISTANCE  NETWORK  ANALYTICAL  MODEL 
OF  CLCP 
RESCALCS.TK 


VARIABLE  SHEET 

St  InputDDDD  NameDDD  OutputDDD  UnitDDDDD  CommentDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 

RESISTANCE  CALCULATIONS  FOR  CLCP 


RFB 

5.3568E-5 

Ohm 

.0000018 

RHOCU 

Ohm- cm 

12 

LT 

cm 

1.27 

FBH 

cm 

.635 

FBW 

cm 

RAB 

.00013392 

Ohm 

1.27 

ABH 

cm 

.254 

ABW 

cm 

RSDW 

.00000048 

Ohm 

.000072 

RHOSS 

Ohm- cm 

.0508 

TW 

cm 

RPDW 

.0000375 

Ohm 

.635 

W 

RNAK 

1.6959E-6 

RHONAK 

4.0702E-5 

1.27 

H 

CM 

373 

TL 

RCEMFEQ 

.0000127 

Ohm 

CEMF 

.00635 

V 

1000 

B 

G 

500 

lA 

A 

1000 

V 

cm/s 

RTP 

.00025 

Ohm 

.0508 

LTP 

cm 

.34559261 

DOTP 

cm 

.3175 

DITP 

cm 

ATP 

cm'' 2 

DOTPE 

.13606008 

in. 

DITPE 

.125 

in. 

R7 

1 

R4 

3.5331751 

R5 

26.458333 

R3 

78.125 

R1 

111.60022 

R2 

279.00056 

R6 

520.83333 

FILE:RESCAL'CS.TK 

Resistance  of  field  busbar 
Resistivity  of  field  busbar 
length  of  busbar 
Field  bubar  height 
Field  bubar  width 

Resistance  of  armature  busbar 
Armature  busbar  height 
Armature  busbar  width 

Resistance  of  series  duct  wall 
Resistivity  if  duct  wall 
Thkn  of  series  duct  wall 

Resistance  of  2-parallel  duct  walls 
Width  of  duct 

Resistance  of  the  fluid  (NaK) 
Resistivity  of  NaK 
Inside  duct  height 


Equivalent  resistance  of  CEMF 
Counter  EMF  of  fluid  in  magnetic  field 
Magnetic  field  density 
Armature  current 
Velocity  of  fluid  in  duct 

Resistance  of  tap  pins 
Length  of  tap  pins 
Outside  diameter  of  tap  pins 
Inside  diameter  of  tap  pins 
Cross-sectional  area  of  tap  pins 


Resistance  ratio  RSDW/RSDW 
Resistance  ratio  NaK/RSDW 
Resistance  ratio  RCEMFEQ/RSDW 
Resistance  ratio  RPDW/RSDW 
Resistance  ratio  RFB/RSDW 
Resistance  ratio  RAB/RSDW 
Resistance  ratio  RTP/RSDW 
RTP  can  be  varied  to  a  ratio  of  104 
or  less. 


FI 


RULE  SHEET 

S  Ru 1 eDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDC: 
"RESISTANCE  OF  CLCP  MATERIALS 

>K  RFB=RHOCU*2*LT/(FBH5t;FBW) 

*  RAB=RHOCU*2*LT/(ABH*ABW) 

*  RSDW=2>KRHOSS*TW/(ABH*LT) 

*  RPDW=RH0SSW(2*TW*LT) 

*  RNAK=RH0NAK*W/(H*LT) 

*  RHONAK=(36. 768-1. 6212*10^-2>t:TL+7. 1741>*‘10^-5*TL^2)*10^-6 

*  CEMF=B>KV*W*10'‘-8 

*  RCEMFEQ=CEMF/IA 
"ATP=PI()>K(DOTP"2-DITP"2)/4 
"RTP=RHOSS*LTP/ATP 
"ATP=PI( )*DOTP^2/4 

*  RTP=RHOSS*LTP/ ( PI ( ) * ( DOTP^2-DITP^ 2 ) /4 ) 

*  DOTPE=DOTP/2 . 54 

*  DITPE=DITP/2.54 

*  R1=RFB/RSDW 

*  R2=RAB/RSDW 

*  R3=RPDW/RSDW 

*  R4=RNAK/RSDW 

*  R5=RCEMFEQ/RSDW 
5*;  R6=RTP/RSDW 

*  R7=RSDW/RSDW 


MMMMMMMMMMMMMMMMMMMM  TABLE  SHEET 
NameDDDDDDDDDD  TitleDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 
DIAMS  DIAMETERS  OF  THE  TAP  PINS 

DIAMS2  ALTERNATE  DIAMETERS  FOR  TAP  PINS 


DIAMETERS  OF  THE  TAP  PINS 


RES 

OD 

ID 

.00005 

.173383194 

.  125 

.0001 

.151140219 

.  125 

.00015 

.142958889 

.125 

.0002 

.138687357 

.125 

.00025 

.136060084 

.  125 

F2 


APPENDIX  G 


NaK-78  PHYSICAL  PROPERTIES  VERSOS  TEMPERATDRS 

HaKPRO.PTY 

DENSITY  (ing/B3) 

ELECTRICAL  RESIST  .. .  ^TY  (mOhm-in) 

THERMAL  CONDUCTIVITY  (W/m-k) 

VISCOSITY  (mPa-s) 


VAPOR  P.  EQN 


VAPOR  PRESSURE  (Pa) 


VARIABLE  SHEET 

St  InputDDDD  NameDDD  OutputDDD  UnitDDDDD  CommentDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 

FILE:NAK78PRO.PTY 

PHYSICAL  PROPERTIES  OF  NAK-78 


3.436541 

AMU 

-.0162033 

BMU 

3. 1836E-5 

CMU 

-2.833E-8 

DMU 

9.375E-12 

EMU 

.5333 

MU781 

cp 

.3663 

MU782 

cp 

.2856 

MU783 

cp 

.2389 

MU784 

cp 

.1964 

MU785 

cp 

373 

TNAKl 

K 

773 

TNAK2 

K 

TNAK3 

K 

TNAK4 

K 

TNAK5 

K 

MU78 

mPa-s 

373 

TNAK 

K 

-.0691676 

AK 

.00187204 

BK 

-4. 377E-6 

CK 

4.8583E-9 

DK 

-2.08E-12 

EK 

.232 

KNAK781 

W/cm-K 

.247 

KNAK782 

W/cm-K 

.256 

KNAK783 

W/cm-K 

.262 

KNAK784 

W/cm-K 

.263 

KNAK785 

W/cm-K 

KNAK78  23.2 

W/m-K 

36.747088 

ARHO 

-.016312 

BRHO 

.000072 

CRHO  • 

4.34E-20 

DRHO 

-1. 97E-23 

ERHO 

40.68 

RH0781 

f Ohm- cm 

45.14 

RH0782 

f Ohm- cm 

51.04 

RH0783 

f Ohm-cm 

58.38 

RH0784 

f Ohm- cm 

67.16 

RH0785 

fOhm-cm 

RH078  406.8 

nOhm-m 

VISCOSITY 

Polynomial  Equation  Coefficicent 
Polynomial  Equation  Coefficicent" 
Polynomial  Equation  Coefficicent 
Polynomial  Equation  Coefficicent 
Polynomial  Equation  Coefficicent 
MU78  in  1st  Polynomial  Equation 

MU78  in  2nd  Polynomial  Equation 

MU78  in  3rd  Polynomial  Equation 

MU78  in  4th  Polynomial  Equation 

MU78  in  5th  Polynomial  Equation 

Temperature  of  NaK  in  1st  Equation 
Temperature  of  NaK  in  2nd  Equation 
Temperature  of  NaK  in  3rd  Equation 
Temperature  of  NaK  in  4th  Equation 
Temperature  of  NaK  in  5th  Equation 

Viscosity  of  NaK~78  at  TNAK 
Temperature  of  NaK 

THERMAL  CONDUCTIVITY 

Polynomial  Equation  Coefficicent 
Polynomial  Equation  Coefficicent 
Polynomial  Equation  Coefficicent 
Polynomial  Equation  Coefficicent 
Polynomial  Equation  Coefficicent 
KNAK78  in  1st  Polynomial  Equation 
KNAK78  in  2nd  Polynomial  Equation 
KNAK78  in  3rd  Polynomial  Equation 
KNAK78  in  4th  Polynomial  Equation 
KNAK78  in  5th  Polynomial  Equation 

Thermal  Conductivity  at  TNAK 

ELECTRICAL  RESISTIVITY  OF  NAK-78 

Polynomial  Equation  Coefficicent 
Polynomial  Equation  Coefficicent 
Polynomial  Equation  Coefficicent 
Polynomial  Equation  Coefficicent 
Polynomial  Equation  Coefficicent 
RH078  in  1st  Polynomial  Equation 
RH078  in  2nd  Polynomial  Equation 
RH078  in  3rd  Polynomial  Equation 
RH078  in  4th  Polynomial  Equation 
RH078  in  5th  Polynomial  Equation 

Electrical  Resistivity  at  TNAK 

DENSITY 


G1 


.847 

DENI 

Mg/m''  3 

Density  in 

eqn  1 

.751 

DEN  2 

Mg/m" 3 

Density  in 

eqn  2 

. 93652 

ADEN 

Polynomial 

Coefficient 

-.00024 

BDEN 

Polynomial 

Coefficient 

DEN 

Mg/m" 3 

Density  of 

NaK-78 

.847 

DENI 

Mg/m-'  3 

Density  in 

eqn  1 

.751 

DEN  2 

Mg/m^  3 

Density  in 

eqn  2 

.93652 

ADEN 

Polynomial 

Coefficient 

-.00024 

BDEN 

Polynomial 

Coefficient 

DEN 

Mg/m^3 

Density  of 

NaK-78 

G2 


RULE  SHEET 

S  Ru 1 eDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 

"PHYSICAL  PROPERTIES  OF  NAK-78 
"VISCOSITY 

"MU7  81  =AMU+BMU*TNAK1  +CMU)XTNAK1  ^  2+DMU*TNAKl "  3+EMU*TNAKl  ^  4 
"MU782=AMU+BMU>t<TNAK2+CMU*TNAK2^2+DMU*TNAK2''3+EMU*TNAK2^4 
"MU783=AMU+BMU*TNAK3+CMU*TNAK3"2+DMU*TNAK3"3+EMU><'TNAK3'‘4 
"  MU7  84= AMU+BMU*TNAK4+CMU*TNAK4  2+DMU*TNAK4  3+EMU*TNAK4  "  4 
"  MU7  8  5=AMU+BMU>>:TNAK5+CMU*TNAK5 ''  2+DMU^TNAK5  "  3  +EMU*TNAK5  "  4 

"MD78=AMU+BMU*TNAK+CMU5i<TNAK''2+DMU*TNAK"3+EMU>t<TNAK"4 

"THERMAL  CODUCTIVITY 

" KNAK7  8 1 = AK+BK*TNAK1 +CK*TNAK1 " 2  +DK*TNAK1 ^  3+EK*TNAKl " 4 
"  KNAK7  8  2  =AK+BK=t^TNAK2+CK*  TNAK2 ''  2 +DK*TNAK2  3+EK*TNAK2 4 
"KNAK783=AK+BK>i'TNAK3+CK>ttTNAK3^2+DK*TNAK3^3+EK^TNAK3^4 
"  KNAK7  8 4  =  AK+BK*  TN  AK4 +CK*  TNAK4  "  2 +DK*TNAK4  3+EK*  TNAK4  ^  4 
"  KNAK7  8  5  =AK+BK*TNAK5 +CK>kTNAK5  "  2  +DK*TNAK5  "  3+EK*TNAK5  "  4 

*  KNAK78=(AK+BK5t=TNAK+CK*TNAK^2+DK*TNAK^3+EK>KTNAK^4)*100 
"ELECTRICAL  RESISTIVITY 

*  RH078=  ( ARHO+BRHO>tcTNAK+CRHO*TNAK''2+DRHO*TNAK"3+ERHO*TNAK''4 ) >^10 
"DENSITY 

"DEN1=ADEN+BDEN*TNAK1 

"DEN2=ADEN+BDEN*TNAK2 

*  DEN=ADEN+BDEN*TNAK 


G3 


SCOSITY  OF  NAK-78,  BLOMBURG 


NaK  TEMPE 


SISTIUITV  OF  HAK-78,  BLOMBURG 


VARIABLE  SHEET 

St  InputDDDD  NameDDD  OutputDDD  UnitDDDDD  CommentDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD! 

FILE:  VAPORP.EQN 

DISKETTE:  CLCP-TECH 

Vapor  Pressure  vs  Temperature 


P 

1000 

Pa 

Pressure 

PO 

.01 

atm 

Pressure 

714 

T 

K 

Temperature 

TP 

1.4005602 

1/K 

Inverse  Temperature 

TO 

441 

deg 

C 

Temperature 

TOP 

2.2675737 

1/C 

Inverse  Temperature 

3. 16228E9 

A 

Coefficient 

10686.297 

B 

Coefficient 

C 

Coefficient 

1000 

PI 

Pa 

Pressure  in 

1st  eqn 

PIO 

atm 

Pressure  in 

1st  eqn 

714 

T1 

K 

Temperature 

in  1st  eqn 

TIO 

441 

deg 

C 

Temperature 

in  1st  eqn 

100 

P2 

Pa 

Pressure  in 

2nd  eqn 

P20 

atm 

Pressure  in 

2nd  eqn 

617 

T2 

K 

Temperature 

in  2nd  eqn 

T20 

344 

deg 

C 

Temperature 

in  2nd  eqn 

10 

P3 

Pa 

Pressure  in 

3rd  eqn 

P30 

atm 

Pressure  in 

3rd  eqn 

546 

T3 

K 

Temperature 

in  3rd  eqn 

T30 

273 

deg 

C 

Temperature 

in  3rd  eqn 

10.151865 

ASb 

11121.563 

BSb 

659.379 

CSb 

32.570051 

ABi 

231760.16 

BBi 

Coefficients 

for  Bi  equation 

6891.4834 

CBi 

11.512364 

ACd 

6189.6255 

BCd 

Coefficients 

for  Cd  equation 

86.026332 

CCd 

9.2825295 

APb 

4614.2737 

BPb 

Coefficients 

for  Pb  equation 

-607.3846 

CPb 

10.959671 

AP 

2756.5896 

BP 

Coefficients 

for  P  equation 

.44140838 

CP 

10.959671 

AS 

3675.4528 

BS 

Coefficients 

for  S  equation 

23.921878 

CS 

98.307981 

ATI 

2538860.4 

BTl 

Coefficients 

for  T1  equation 

25762.936 

CTl 

11.699187 

ASn 

20397.825 

BSn 

Coefficients 

for  Sn  equation 

559.66503 

CSn 

G8 


34.814147  AZn 
171873.22  BZn 
4786.0722  CZn 


Coefficients  for  Zn  equation 


16.613534  AMg 

29055.626  BMg  Coefficients  for  Mg  equation 

1356.0747  CMg 

3.16228E9  ANaK78  Coefficients  for  NaK-78  equation 

10686.297  BNaK78 


G9 


RULE  SHEET 

S  RuleDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD: 


"Fluid  Properties 

"Vapor  pressure  curve  coefficients 


*  T=T0+273 

*  TP=(1/T)>I'10^3 

*  T0P=(1/T0)*10^3 

*  Tl=T10+273 

*  T2=T20+273 

*  T3=T30+273 

*  P=P0*10^5  "10^5  Pa/atm 

"P1=P10*10^5 

"P2=P20*10"5 

"P3=P30*10^5 

"P=P0*10''6 

"P1=P10*10^6 

"P2=P20*10'^6 

"P3=P30*10"6 

" LOG ( P ) =ASb-BSb/ ( T+CSb ) 

" LOG ( PI ) =ASb-BSb/ ( Tl+CSb ) 

" LOG ( P2 ) = ASb-BSb/ ( T2+CSb ) 

" LOG ( P3 ) =ASb-BSb/ ( T3+CSb ) 

" LOG ( P ) =ABi-BBi/ ( T+CBi ) 
"LOG(Pl)=ABi-BBi/(Tl+CBi) 

" LOG ( P2 ) =ABi-BBi/ ( T2+CBi ) 

" LOG ( P3 ) =ABi-BBi/ ( T3+CBi ) 

" LOG ( P ) = ACd-BCd/ ( T+CCd ) 

" LOG ( PI ) =ACd-BCd/ ( Tl+CCd ) 

" LOG ( P2 ) =ACd-BCd/ ( T2+CCd ) 

" LOG ( P3 ) =ACd-BCd/ ( T3+CCd ) 

" LOG ( P ) = APb-BPb/ ( T+CPb ) 
"LOG(Pl)=APb-BPb/(Tl+CPb) 

" LOG ( P2 ) = APb-BPb/ ( T2+CPb ) 

" LOG ( P3 ) =APb-BPb/ ( T3+CPb ) 

"LOG(P)=AP-BP/(T+CP) 

"L0G(P1)=AP-BP/(T1+CP) 

" LOG ( P2 ) = AP-BP/ ( T2+CP ) 
"LOG(P3)=AP-BP/(T3+CP) 

" LOG ( P) =AS-BS/ ( T+CS ) 

" LOG ( PI ) =AS-BS/ ( T 1 +CS ) 

" LOG ( P2 ) =AS-BS/ ( T2+CS ) 

" LOG ( P3 ) =AS-BS/ ( T3+CS ) 


"L0G(P)=AT1-BT1/(T+CT1) 
"L0G(P1)=AT1-BT1/(T1+CT1) 
"  LOG ( P2 ) =AT1-BT1/ ( T2+CT1 ) 


GIO 


"L0G(P3)=AT1-BT1/(T3+CT1) 

" LOG ( P ) =ASn-BSn/ ( T+CSn ) 

” LOG ( PI ) =ASn-BSn/ ( T 1 +CSn ) 
" LOG ( P2 ) =ASn-BSn/ ( T2+CSn ) 
" LOG ( P3 ) =ASn-BSn/ ( T3+CSn ) 

■'LOG(P)=AZn-BZn/(T+CZn) 
"LOG(Pl)=AZn-BZn/(Tl+CZn) 
" LOG ( P2 ) =AZn-BZn/ ( T2+CZn ) 
" LOG ( P3 ) =AZn-BZn/ ( T3+CZn ) 

" LOG ( P ) =AMg-BMg/ ( T+CMg ) 

" LOG ( PI ) =AMg-BMg/ ( T 1 +CMg ) 
" LOG ( P2 ) =AMg-BMg/ ( T2+CMg ) 
" LOG ( P3 ) =AMg-BMg/ ( T3+CMg ) 
P=A/EXP(B/T) 

PI  =  A/EXP  (B/Tl) 

"P2= A/EXP (B/T2) 
"P3=A/EXP(B/T3) 


"T=B/LN(A/P) 


APOR  PRESSURE  OF  NAK 


